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Tunneling  Transistors  Based  on 
Graphene  and  2-D  Crystals 

Graphene-based  tunneling  transistors  and  how  these  compare  to  2-D  transistors  made 
from  the  GaAs/AlGaAs  materials  systems  is  the  topic  of  discussion  in  this  paper. 

By  Debdeep  Jena,  Member  IEEE 


ABSTRACT  |  As  conventional  transistors  become  smaller  and 
thinner  in  the  quest  for  higher  performance,  a  number  of 
hurdles  are  encountered.  The  discovery  of  electronic-grade  2-D 
crystals  has  added  a  new  “layer”  to  the  list  of  conventional 
semiconductors  used  for  transistors.  This  paper  discusses  the 
properties  of  2-D  crystals  by  comparing  them  with  their  3-D 
counterparts.  Their  suitability  for  electronic  devices  is  dis¬ 
cussed.  In  particular,  the  use  of  graphene  and  other  2-D  crystals 
for  interband  tunneling  transistors  is  discussed  for  low-power 
logic  applications.  Since  tunneling  phenomenon  in  reduced 
dimensions  is  not  conventionally  covered  in  texts,  the  physics 
is  developed  explicitly  before  applying  it  to  transistors.  Though 
we  are  in  an  early  stage  of  learning  to  design  devices  with  2-D 
crystals,  they  have  already  been  the  motivation  behind  a  list  of 
truly  novel  ideas.  This  paper  reviews  a  number  of  such  ideas. 

KEYWORDS  |  Graphene;  semiconductors;  transistor;  tunneling 

I.  INTRODUCTION 

Semiconductors  come  in  many  crystal  forms.  Since  their 
discovery  in  the  early  20th  century,  the  semiconductors 
used  in  electronic  and  optical  devices  are  of  the  3-D  crystal 
form.  Three-dimensional  crystal  semiconductors  have  re¬ 
mained  at  the  heart  of  such  devices  from  the  earliest 
“cat’s  whisker”  detectors  [1]  to  the  latest  billion-transistor 
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silicon  complementary  metal-oxide-semiconductor 
(CMOS)  [2],  [3]  and  quantum-well  (QW)  lasers  [4].  As 
the  understanding  of  the  physics  of  electron  transport  and 
electron-photon  coupling  sharpened,  it  became  clear  that 
controlling  the  potential  energy  landscape  of  electrons 
could  lead  to  massive  boosts  in  device  functionality  and 
performance. 

The  first  level  of  direct  control  of  the  “energy-band 
diagrams”  was  by  chemical  doping,  which  involved  replac¬ 
ing  a  small  number  of  atoms  of  the  3-D  semiconductor  by 
those  with  higher  or  lower  valence.  The  next  advance 
involved  varying  the  chemical  nature  of  the  crystal  along 
specific  directions,  which  marked  the  birth  of  semiconduc¬ 
tor  heterostructures  [5].  These  advances  taught  electrons 
“new  tricks,”  and  made  possible  the  smallest  and  fastest 
electronic  switches  [6],  high-density  memories,  and  the 
most  efficient  light-emitting  diodes  (LEDs)  and  lasers  [7]. 
These  devices  form  the  bedrock  of  computation,  data 
storage,  solid-state  lighting,  and  communication  in  today’s 
information  age. 

At  this  time,  in  the  early  part  of  the  21st  century,  these 
building  blocks  based  on  traditional  device  concepts  are 
approaching  their  performance  limits.  Therefore,  new 
ideas  and  new  materials  are  necessary.  For  example, 
photonic  crystal,  metamaterial,  and  plasmonic  concepts 
are  advancing  the  area  of  optoelectronic  devices  beyond 
what  was  thought  possible  before  [8],  [9].  Strong  light- 
matter  interaction  has  been  exploited  to  demonstrate 
polariton  lasers  that  take  advantage  of  Bose-Einstein 
condensation  at  room  temperature  for  ultralow  threshold 
lasing  [10]. 

Similarly,  for  electronic  switching  devices,  a  number  of 
approaches  are  being  taken  to  address  the  future  beyond 
scaling.  Conventional  field-effect  and  bipolar  transistors 
operate  on  the  basis  of  energy  filtering  of  electrons  (or 
holes)  flowing  over  a  barrier.  The  barrier  is  electrosta¬ 
tically  controlled  with  a  voltage.  In  an  electrostatically 
well-designed  device,  all  of  the  control  voltage  is  spent  in 
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moving  the  barrier.  The  electrons  carrying  the  current 
are  spread  in  a  band  according  to  the  Fermi-Dirac 
distribution,  with  a  Boltzmann  tail  in  energy.  The  energy 
filtering  thus  leads  to  a  current  dependence  of  the  form 
I  ~  exp [qV/kfiT],  where  q  is  the  electron  charge,  V  is  the 
voltage,  T  is  the  temperature,  and  ks  is  the  Boltzmann 
constant.  When  operated  in  this  fashion,  the  current 
cannot  be  changed  any  steeper  than  S  ~  (fegT/q)  In  10  ~ 
60  mV/decade.  This  subthreshold  swing  (SS)  “limit”  is 
often  referred  to  as  the  “thermal  limit”  or  the 
“Boltzmann  limit”  (though  Boltzmann  did  not  set  this 
limit).  We  refer  to  this  condition  as  the  SS  limit  to  avoid 
confusion. 

An  electronic  switch  must  have  its  on-  and  OFF-states 
clearly  demarcated  for  performing  digital  (Boolean)  logic. 
Let  us  say  this  demarcation  is  set  to  Ion/Joff  =  104.  To 
achieve  it,  a  voltage  supply  of  at  least  4  x  60  mV  =  0.24  V 
is  necessary.  Since  the  speed  of  switching  and  the  dynamic 
and  static  power  dissipation  of  transistors  are  strong 
functions  of  the  supply  voltage,  the  SS  limit  sets  a  floor  of 
minimum  power  dissipation.  This  issue  is  described  in 
sufficient  detail  in  a  number  of  recent  articles  that  moti¬ 
vate  the  search  for  new  materials  and  ideas  for  going 
beyond  the  SS  limit  [11]— [13]. 

Now  there  is  nothing  particularly  fundamental  about 
the  SS  limit.  Devices  that  do  not  operate  on  the  traditional 
transistor  mechanism  exist  today  and  operate  below  the  SS 
limit.  An  example  is  a  nanoelectromechanical  system 
(NEMS),  which  is  the  analog  of  a  mechanical  relay.  Sub¬ 
stantial  progress  has  been  made  in  this  area  [14].  Due  to 
mechanical  moving  parts,  these  devices  are  currently  slow, 
but  are  expected  to  improve  with  scaling. 

A  number  of  relatively  new  ideas  are  being  explored  at 
this  time  for  switching  devices  beyond  the  SS  limit.  Some 
exploit  impact  ionization  to  obtain  sub-SS  limit  operation 
[15],  [16].  Other  devices  aim  to  use  correlated  electron 
effects;  for  example,  if  electrons  can  be  made  to  “pair  up” 
similar  to  Cooper  pairs  in  superconductors,  but  at  room 
temperature,  the  SS  limit  would  be  cut  in  half.  If  the 
control  voltage  could  be  internally  “stepped  up”  through 
novel  ferroelectric  gates,  sub-SS  limit  devices  can  be 
realized  [17].  Other  routes  involve  the  internal  trans¬ 
duction  of  the  voltage  into  other  state  variables  such  as 
strain,  spin,  or  electron  localization  [18].  Among  these 
strategies,  a  transistor  concept  based  on  interband  tun¬ 
neling  transport  has  emerged  as  an  attractive  candidate 
for  switching.  This  paper  will  focus  on  this  device.  The 
tunneling  field-effect  transistor  (TFET)  can  be  realized 
in  traditional  3-D  crystal  semiconductors  and  their 
heterostructures. 

However,  since  the  discovery  of  graphene  in  2004, 
device  engineers  have  a  new  class  of  materials  in  2-D 
crystals  at  their  disposal.  In  this  paper,  we  discuss  pos¬ 
sible  realizations  of  TFETs  with  2-D  crystals,  and  com¬ 
pare  them  with  3-D  crystal  counterparts.  In  the  process 
of  this  discussion,  a  number  of  novel  features  of  2-D 


crystals  will  emerge  that  distinguish  them  from  tradi¬ 
tional  3-D  crystal  semiconductors.  These  novel  features 
of  the  new  material  family  offer  a  compelling  case  for 
investigating  them  further.  To  motivate  their  suitability 
for  electronic  devices,  we  first  discuss  the  various  2-D 
crystal  materials  and  their  properties.  We  do  so  against 
the  backdrop  of  their  ubiquitous  3-D  crystal  semicon¬ 
ductor  counterparts. 

II.  TWO-DIMENSIONAL  CRYSTALS 

Two-dimensional  crystals  exploded  into  the  limelight  in 
2004  with  the  remarkable  reports  of  the  isolation  of  ato¬ 
mically  thin  graphene  [19]-[21].  What  is  often  overlooked 
is  that  the  early  reports  [22]  also  presented  evidence  of  the 
isolation  of  single-layers  of  BN — an  insulator  or  a  wide- 
bandgap  semiconductor,  M0S2 — a  traditional  semicon¬ 
ductor,  and  NbSe2 — a  superconductor  with  possible 
charge-density  wave  electronic  phases.  Single  layers  of 
the  cuprate  high-Tc  superconductors  were  also  isolated.  It 
is  interesting  to  note  that  the  voltage  “scaling”  of  silicon 
CMOS  processors  stalled  around  the  same  time,  marking 
the  move  toward  multicore  processors  [23].  One  of  the 
reasons  for  paradigm  shift  was  the  unsustainable  increase 
in  dynamic  and  OFF-state  power  dissipation  due  to  the  SS 
limit  and  high-frequency  operation.  Whether  2-D  crystals 
can  help  in  this  arena  remains  to  be  seen.  We  first  discuss  a 
few  properties  of  2-D  crystals  and  their  suitability  for 
electronic  devices. 

Fig.  1  is  a  schematic  representation  of  the  structure  of 
crystals  of  various  dimensions.  The  bottom  row  shows  the 
atomic  building  blocks.  The  first  column  shows  the  ubi¬ 
quitous  3-D  crystal  semiconductors.  The  second  column 
shows  the  emerging  family  of  2-D  crystals  and  their 
many  variants.  The  third  and  fourth  columns  indicate 
ideal  1-D  and  0-D  structures.  Atomic  chains  have  been 
investigated  for  their  transport  properties  [24],  and  a 
benzene  ring  can  be  considered  either  as  an  atomic 
“ring,”  or  even  a  basic  2-D  crystal  unit.  An  atom  is  a 
perfect  0-D  structure  in  which  electrons  are  localized  in  all 
three  dimensions.  We  note  that  the  electrons  in  an  atom 
still  move  in  3-D,  but  their  energy  spectra  are  discrete  and 
gapped;  they  do  not  form  bands  that  are  necessary  for 
transport.  It  is  in  this  sense  that  they  are  0-D.  We  focus  our 
attention  on  2-D  crystals,  and  their  differences  from  3-D 
crystal  semiconductors. 

The  building  blocks  for  3-D  semiconductors  are  typi¬ 
cally  tetrahedrally  bonded  atoms.  The  lattice  is  3-D,  and 
the  basis  typically  consists  of  two  atoms.  For  example, 
electrons  in  3-D  crystals  from  group  IV  elements  (Si,  Ge, 
etc.)  occupy  [core]  ms2mp2  orbitals,  where  m  is  the  row 
number  in  the  periodic  table,  and  [core]  represents  the 
core  electrons  that  do  not  participate  in  chemical 
bonding  directly.  Electrons  from  the  outermost  s  and  p 
orbitals  of  nearest  neighbor  atoms  pair  up  to  form  sp3 
bonds.  An  sp3  bond  is  inherently  3-D,  and  so  is  the 
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Fig.  1.  A  schematic  representation  of  “crystals”  of  the  many  spatial  dimensions  that  result  from  various  building  blocks.  The  building  blocks 
contain  atomic  bases  that  form  3-D  bonds  in  the  first  column ,  2-D  planar  bonds  in  the  second  column,  and  1-D  linear  bonds  in  the  third  column. 
The  ideal  O-D  structure  is  an  atom  in  the  fourth  column. 


resulting  semiconductor  crystal.  The  natural  crystal  is 
thus  a  bulk  3-D  semiconductor.  A  termination  such  as  a 
surface  results  in  dangling  bonds,  a  fraction  of  which 
might  reconstruct. 

The  corresponding  building  block  of  a  2-D  crystal 
consists  of  a  planar  2-D  lattice.  For  graphene  and  BN,  the 
basis  consists  of  two  atoms  attached  to  a  hexagonal  planar 
lattice.  These  chemical  bonds  in  the  two-atom  basis  for 
graphene  and  BN  are  of  the  sp2  type.  So  the  chemical 
bonds  of  their  basis  are  also  planar.  In  the  second  column 
of  Fig.  1,  the  underlying  planar  structure  of  2-D  crystals  is 
shown.  Attached  to  each  point  of  intersection  is  one 
carbon  atom  for  graphene,  alternating  B  and  N  atoms  for 
BN,  and  a  basis  of  X-M-X  for  transition  metal  dichal- 
cogenides  (TMDs).  TMD  2-D  crystals  share  the  same 
planar  lattice  geometry  of  graphene  and  BN.  But  the  basis 
of  TMD  2-D  crystals  consists  of  three  atoms  of  the  form 
MX2,  where  M  is  the  transition  metal  chemically  bonded 
to  two  chalcogenide  atoms  X.  The  chemical  bonds  in 
TMD  2-D  crystals  (e.g.,  M0S2,  WSe2,  WS2,  etc.)  involve 
s-,  p-,  and  d-orbitals,  and  the  two  M-X  bonds  stick  out  of 

DISTRIBUTION  A:  Distribution 


the  center  2-D  plane  containing  the  transition  metal  atom 
M  [25].  Thus,  unlike  its  lattice,  the  basis  of  TMD  2-D 
crystals  is  not  perfectly  planar.  Recent  reports  also  in¬ 
dicate  the  possible  existence  of  2-D  forms  of  Si  (silicene), 
Ge  (germanene),  and  possibly  AlN  and  GaN  [26]— [28], 
[102].  Single  layers  of  2-D  crystals  are  typically  less  than 
1  nm  in  thickness.  An  exotic  form  of  a  2-D  crystal  semi¬ 
conductor  may  also  exist  when  two  surfaces  of  topolog¬ 
ical  insulators  come  close  to  each  other  [29].  These 
materials  have  been  less  explored  than  the  others  dis¬ 
cussed  here. 

Unlike  a  perfect  3-D  crystal,  a  perfect  2-D  crystal  has 
no  broken/dangling  bonds  on  its  surface.  The  quasi-low- 
dimensional  structures  formed  from  3-D  crystals  such  as 
2-D  nanomembranes,  1-D  nanowires,  and  0-D  nano¬ 
crystals  are  still  volume  elements  deriving  from  3-D 
bonding,  and  necessarily  have  dangling  bonds  on  their 
surfaces.  These  broken  bonds  may  be  passivated  by  either 
dielectrics,  or  by  lattice-matched  or  strained  heterostruc¬ 
tures.  In  contrast,  the  various  dimensional  structures  de¬ 
riving  from  2-D  crystals  are  “hollow”  and  are  “all-surface.” 
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Two-dimensional  crystal  sheets  may  be  stacked  to  form 
3-D  structures  with  weak  van-der-Waal’s  interlayer 
bonding.  They  can  be  rolled  up  into  quasi-l-D  nanotubes, 
or  into  0-D  buckyballs  (Ceo)-  The  symmetry  of  a  2-D 
crystal  is  broken  at  its  edge.  Similar  to  the  surface  state 
reconstruction  or  passivation  of  the  surfaces  of  3-D  crys¬ 
tals,  the  edge  states  can  reconstruct  and  tie  up  the  dangling 
bonds.  For  special  cases,  such  as  in  buckyballs,  the  chemi¬ 
cal  bonding  is  seamless  and  there  are  no  broken  bonds. 
Indeed,  the  icosahedral  geometry  of  the  buckyball  belongs 
to  one  of  the  five  platonic  solids,  which  have  mathemat¬ 
ically  represented  “perfection”  in  shape  since  the  earliest 
times  [30],  [31]. 

For  electronic  devices  using  field  effect,  the  absence  of 
dangling  bonds  is  a  major  advantage  for  planar  2-D  crystals, 
since  electrons  trapped  in  them  serve  to  shield  electric 
field  lines  from  entering  the  bulk  of  the  corresponding  3-D 
semiconductors.  We  now  discuss  the  electronic  properties 
of  2-D  crystals  and  compare  them  to  those  of  3-D  crystal 
semiconductors. 

A.  Electronic  Properties  of  2-D  Crystals 

The  electronic  orbitals  that  form  the  family  of  2-D 
crystals  are  shown  in  Fig.  2.  Electron  states  at  the  conduc¬ 
tion  and  valence  band  edges  of  3-D  semiconductors  derive 
from  various  admixtures  of  sp3  bonds.  For  direct-gap 
semiconductors  such  as  GaAs  and  GaN,  the  conduction 
band  edge  is  mostly  s-like.  The  spherical  symmetry  of  the 
s-orbitals  imparts  electrons  in  the  conduction  band  their 
isotropic  nature.  The  electronic  states  at  the  valence 
band  edge  on  the  other  hand  are  more  p-like.  Because 
p-orbitals  are  directional,  the  hole  effective  mass  is 


anisotropic.  The  imbalance  of  the  nature  of  chemical 
bonding  in  3-D  crystals  semiconductors  thus  also  results 
in  an  asymmetry  in  the  curvature  or  the  effective  mass  of 
the  conduction  and  valence  band  states.  In  modern 
complementary  logic  devices,  symmetry  is  a  highly  de¬ 
sirable  characteristic.  The  degree  of  asymmetry  between, 
for  example,  nMOS  and  pMOS  devices  dictates  the 
geometry  and  layout  of  circuits  that  could  be  considerably 
simplified  by  symmetry. 

The  covalent  bonds  in  graphene  and  BN  are  of  the  sp2 
kind.  They  are  responsible  for  the  structural  properties  of 
the  crystal.  The  leftout  pz  orbital  sticks  out  of  the  2-D 
plane.  The  electrons  in  these  orbitals  can  hop  between 
nearest  neighbors,  leading  to  the  electronic  conductivity 
and  optical  properties  of  such  crystals.  In  graphene  and 
BN,  the  structural  properties  such  as  thermal  conductivity 
and  mechanical  stability  derive  from  the  covalent  sp2 
bonds.  But  the  electronic  and  optical  properties  derive 
from  the  delocalized  pz  orbitals.  There  is  a  wide  energy 
separation  between  the  sp2  and  pz  energy  bands.  In  this 
sense,  the  electronic  properties  of  such  2-D  crystals  have  a 
different  origin  than  their  structural  properties.  This  is  in 
contrast  to  3-D  semiconductors,  where  the  structural  and 
electronic  properties  derive  from  the  same  sp3  electronic 
band  states. 

Electrons  in  3-D  crystals  can  be  quantum-mechanically 
confined  to  move  in  2-D  and  1-D,  or  localized  in  0-D  by 
chemical  and  geometrical  constraints  in  heterostructures, 
as  shown  in  the  first  column  in  Fig.  1.  This  is  achieved  by 
taking  advantage  of  energy  band  offsets  around  the 
bandgap.  Conduction  band  offsets  AEq  confine  electrons, 
and  valence  band  offsets  AEV  confine  holes.  We  note  here 
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Fig.  2.  Energy  band  alignments  of  various  2-D  crystals  compared  to  silicon.  The  relative  energy  band  offsets  of  graphene,  BN,  and  transition-metal 
dichalcogenides  are  shown.  The  numbers  at  the  center  indicate  the  respective  bandgaps  reported  at  this  time,  but  are  subject  to  refinement 
with  further  experiments.  An  energy  scale  from  the  vacuum  level  is  also  indicated,  showing  a  work  function  (or  electron  affinity)  of  intrinsic 
zero-gap  2-D  graphene  to  be  4.5  eV.  The  conduction  and  valence  band  edge  states  of  Si,  graphene,  and  BN  are  formed  of  linear  combinations 
of  |s  >-  and  |p  -orbitals ,  whereas  those  of  the  transition-metal  dichalcogenide  2-D  crystals  involve  |d  >-orbital  states  at  the  band  edges. 

The  presence  of  6-orbital  states  near  the  Fermi  level  implies  that  some  of  these  2-D  crystals  can  exhibit  electronic  phenomena  that  require 
many-particle  effects  such  as  magnetism  and  superconductivity. 
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that  the  2-D  confinement  of  electrons  in  a  quantum  well  in 
a  3-D  crystal  leads  to  a  quasi-2-D  electron  gas  (2-DEG). 
This  means  there  are  multiple  2-D  electronic  subbands 
whose  spacing  in  energy  grows  as  the  inverse  square  of  the 
spatial  confinement.  In  sharp  contrast,  there  is  just  one 
band  for  2-D  electron  systems  in  single-layer  2-D  crystals, 
since  the  electron  wave  function  cannot  spread  sufficiently 
out  of  the  plane  in  equilibrium. 

The  energy  bandgaps  and  the  band  lineups  of  a  few  2-D 
crystals  are  shown  in  Fig.  2.  The  figure  also  indicates  the 
chemical  bonding  schemes  that  characterize  them,  along 
with  their  relative  positions  with  respect  to  the  vacuum 
energy  level  [32]  — [34] .  A  distinctive  feature  of  the  2-D 
crystals  is  that  their  energy  gap  windows  are  not  populated 
by  surface  states  in  sufficiently  crystalline  sheets,  as  is 
necessarily  the  case  for  3-D  crystals.  Thus,  the  measure¬ 
ments  of  their  band  alignments  are  relatively  simpler,  as 
described  in  [33]. 

Graphene  is  a  zero-bandgap  semiconductor,  with  the 
energy  dispersion  E(kx,ky)  =  |k| ,  where  H  is  the  re¬ 

duced  Planck’s  constant,  vF  =  108  cm/s  is  called  the  Fermi 
velocity,  and  |k|  =  +  kj  =  27r/A  is  the  electron  wave 

vector.  The  dispersion  is  taken  around  the  Dirac  points  in 
the  band  structure,  which  are  located  at  the  twofold  de¬ 
generate  K-points  in  the  k-space,  as  shown  in  Fig.  3.  These 
states  are  similar  to  the  conduction  band  edge  or  valence 
band  edge  states  of  3-D  semiconductors.  The  positive 
branch  is  the  conduction  band,  and  the  negative  branch  is 
the  valence  band.  We  note  the  perfect  symmetry  of  the 
bands,  which  is  quite  distinct  from  traditional  3-D  semi¬ 
conductors.  This  symmetry  is  special,  and  has  an  important 
bearing  on  tunneling  transistors  discussed  later.  The  ener¬ 
gy  bandgap  is  zero.  The  density  of  states  (DOS)  of  2-D 
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Fig.  3.  The  k-space  picture  of  2-D  crystals  such  as  graphene,  BN, 
and  the  transition-metal  dichalcogenide  MX,  compounds.  A  good 
understanding  of  the  k-space  picture  is  important  for  choosing  the 
right  materials  for  device  applications,  and  especially  important  for 
tunneling  transistors.  Since  the  real-space  lattice  is  hexagonal  in  the 
2-D  plane,  so  is  the  k-space  lattice.  Since  the  interlayer  separation  is 
larger  than  the  in-plane  lattice  constant,  the  hexagonal  Brillouin  zone 
is  shorter  in  the  vertical  direction.  The  important  high-symmetry 
points  are  labeled.  Graphene,  BN,  and  single-layer  MoS2  have  their 
conduction  band  edge  and  valence  band  edges  at  the  K -points,  which 
leads  to  twofold  degeneracy.  The  conduction  band  edge  of  multilayer 
MoS2  at  this  point  is  believed  to  be  along  the  r  -  K  minimum  as  shown, 
which  makes  it  an  indirect-bandgap  semiconductor,  and  imparts 
to  it  a  valley  degeneracy  of  6  by  symmetry,  similar  to  silicon. 


graphene  is  given  by  p^D(E)  =  [gsgv/27r(RvF)2]  X  |E|, 
where  gs  =  2  is  the  spin  degeneracy  and  gv  =  2  is  the 
valley  degeneracy  [35]. 

Two-dimensional  BN  has  an  energy  bandgap  of 
~6.0  eV  as  a  consequence  of  the  broken  crystal  symmetry 
in  the  basis,  but  its  band  extrema  also  occur  at  the  K- points 
in  the  Brillouin  zone.  Thus,  it  has  the  same  valley  de¬ 
generacy  as  graphene.  The  effective  mass  characterizing 
the  symmetric  conduction  and  valence  bands  of  2-D  BN  is 
m*  ~  0.6mo,  where  mo  is  the  free-electron  mass  [36].  The 
DOS  of  2-D  BN  looks  like  those  of  conventional  2-DEGs, 
Pbn  (E)  =  (gsgvm*/nH2)  x  0[E  -  Ec].  The  major  differ- 
ence  is  the  absence  of  higher  subbands  owing  to  the 
absence  of  atoms  out  of  the  plane. 

The  bandstructures  of  2-D  crystal  semiconductors  of 
the  TMD  family  are  being  evaluated  at  this  time  [37],  [38]. 
Initial  experiments  and  theoretical  models  point  out  that 
they  too  have  their  band  extrema  at  the  K-points.  The 
conduction  and  valence  bands  in  single-layer  TMDs  appear 
less  symmetric  than  graphene  and  BN,  but  much  more 
symmetric  than  traditional  3-D  semiconductor  crystals.  Ef¬ 
fective  masses  ranging  from  m*  ~  0.34mo  —  0.76mo  have 
been  calculated,  and  are  expected  to  undergo  refinement 
through  experimental  measurements  [39]. 

Two-dimensional  crystal  sheets  typically  occur  in 
nature  in  their  stacked  layered  forms.  The  band  structures 
of  the  multilayer  variants  of  graphene,  BN,  and  TMDs  are 
distinct  from  the  single-layer  counterparts.  The  bandgap  of 
a  stacked  2-D  crystal  is  smaller  than  the  single  layer  [40], 
[41].  For  example,  graphite  becomes  a  semimetal  with  a 
negative  bandgap.  Similarly,  when  2-D  crystals  are  used  to 
form  1-D  nanotubes  (Fig.  1),  quasi-l-D  subbands  appear, 
and  the  bandgaps  increase  due  to  additional  quantum 
confinement.  The  DOS  then  acquires  van  Hove  singular¬ 
ities  in  a  manner  similar  to  quasi-l-D  nanowires  or  quan¬ 
tum  wires  formed  of  3-D  semiconductor  crystals.  In  this 
paper,  we  maintain  focus  on  single-layer  2-D  crystals  and 
occasionally  mention  their  quasi-l-D  and  quasi-3-D  va¬ 
riants  when  they  appear  in  context.  The  discussion  of  the 
electronic  band  structures  of  2-D  crystals  leads  us  naturally 
to  a  point  where  we  can  gauge  their  suitability  for  electro¬ 
nic  devices.  We  start  by  discussing  their  suitability  for 
traditional  field-effect  transistors  (FETs). 

B.  Suitability  of  2-D  Crystals  for 
Traditional  Transistors 

The  operation  of  a  FET  hinges  on  electrostatics  and 
transport  of  charge  carriers.  FETs  based  on  3-D  crystal  se¬ 
miconductors  have  been  scaled  to  ~10  s  of  nanometer 
channel  lengths  in  the  quest  to  achieve  higher  perfor¬ 
mance.  As  the  source/drain  separations  have  been  scaled,  it 
has  become  necessary  to  reduce  the  channel  thickness.  This 
requirement  is  driven  by  the  need  for  a  gate  metal  to 
exercise  electrostatic  control  over  mobile  electrons  and 
holes.  If  the  gate  is  farther  away  from  the  carriers  than  the 
S/D  distance,  it  loses  control  over  them.  The  device  then 
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cannot  be  switched  on  or  off  as  effectively  as  is  needed  for 
the  transistor  to  operate  in  a  circuit.  This  necessity  is  at  the 
root  of  the  reason  for  the  move  to  silicon-on-insulator  (SOI) 
and  FinFET  type  of  topologies  [42].  The  silicon  channels 
have  thus  become  more  2-D  in  SOI  structures  [43],  and 
closer  to  1-D  in  FinFETs  and  nanowire  geometries. 

A  quantitative  statement  of  the  importance  of  electro¬ 
statics  is  obtained  from  a  solution  of  the  Poisson  equation 
for  a  FET.  For  a  FET  with  a  semiconductor  layer  of 
thickness  ts  of  dielectric  constant  es  gated  through  an  in¬ 
sulator  of  thickness  tox  and  dielectric  constant  eox,  the 
Poisson  equation  for  the  electric  potential  V  takes  the  form 
d2xv  V/l2,  where  l  =  \ftstox{es/  eox)  is  the  characteristic 
“scaling  length”  [44].  This  length  determines  the  smallest 
distances  over  which  electric  potential  may  be  dropped. 
Therefore,  for  scaling  to  the  smallest  lengths,  high-fC  insu¬ 
lators  and  ultrathin  channels  are  desirable.  This  argument, 
in  conjunction  with  the  absence  of  dangling  bonds  and  the 
associated  interface  traps  highlights  the  attractive  feature 
of  2-D  crystals  for  ultrascaled  FETs  based  on  electrostatics 
arguments  alone.  Furthermore,  2-D  crystal  insulators  such 
as  BN  can  eliminate  dangling  bonds  altogether  in  planar 
FET  geometries. 

As  the  channels  of  3-D  semiconductors  are  thinned 
down,  the  roughness  of  the  surfaces  causes  degradation  of 
the  carrier  transport  due  to  surface-roughness  scattering. 
The  root  of  this  form  of  scattering  is  the  effect  of  the 
roughness  on  the  quantization  of  energy  levels.  For  exam¬ 
ple,  in  a  SOI  structure  of  thickness  t,  the  quantization 
energy  of  subbands  varies  as  E  ~  h2 /m*t2.  Variation  of  the 
layer  thickness  by  At  leads  to  a  perturbation  of  the  subband 
edge  by  AE  ~  (2h2/m*t3)At.  Since  the  scattering  rate  is 
proportional  to  the  square  of  the  perturbation,  the  mobility 
degrades  as  fi  ~  t6,  i.e.,  roughly  as  the  sixth  power  of  the 
width  [45].  Thus,  for  very  thin  layers  of  a  3-D  semicon¬ 
ductor,  such  as  those  used  in  ultrathin  body  (UTB)  tran¬ 
sistors,  the  transport  properties  suffer  from  the  surface 
roughness.  Two-dimensional  crystals  offer  an  ideal  solu¬ 
tion  to  this  problem.  Two-dimensional  crystals  are  intrin¬ 
sically  of  an  atomically  thin  body  (ATB)  nature.  When 
sufficiently  pure,  they  do  not  have  surface  roughness.  The 
attractiveness  of  TMD  2-D  crystal  semiconductors  was 
brought  to  sharp  focus  with  the  demonstration  of  single¬ 
layer  M0S2  FETs  [46].  A  FET  with  108  on/off  ratio  at 
room  temperature  and  electron  mobility  of  ~200  cm2/Vs 
was  achieved  with  a  single  layer  of  M0S2  2-D  crystal  of 
thickness  <  1  nm.  The  SS  was  close  to  ideal,  thanks  to  the 
absence  of  broken  bonds  and  associated  interface  traps. 
Such  performance  has  never  been  measured  in  devices 
made  from  3-D  crystals  of  the  same  thickness.  Though  the 
initial  results  look  promising,  the  dynamic  range  and  re¬ 
liability  of  the  performance  metrics  will  be  assessed  care¬ 
fully  in  the  next  few  years. 

Additional  novel  features  of  charge  transport  in  2-D 
crystals  that  have  been  predicted  and  recently  observed 
include  dielectric-mediated  carrier  mobilities.  The  basic 


premise  is  that  the  Coulomb  interaction  V  ~  q/4i rer  be¬ 
tween  charged  impurities  and  mobile  channel  carriers  is 
mediated  by  the  dielectric  constant  5  of  the  space  sepa¬ 
rating  them.  In  3-D  semiconductors,  the  Coulomb  inter¬ 
action  is  dominated  by  the  bulk  dielectric  constant  of  the 
semiconductor  itself  (i.e.,  £  =  es)  since  the  charged 
impurity  and  the  charge  carrier  are  effectively  buried  in¬ 
side  and  in  close  proximity.  On  the  other  hand,  in  2-D 
crystals,  most  of  the  electric  field  lines  connecting  the 
charged  impurity  to  the  mobile  carrier  actually  lie  outside 
the  2-D  crystal  itself,  in  the  surrounding  dielectric.  This 
effectively  provides  an  external  knob  to  damp  Coulomb 
scattering  and  improve  carrier  mobilities,  since  e  ~  £ox 
for  this  interaction  [47].  Use  of  high-K  dielectrics  has 
been  observed  to  damp  scattering  and  improves  charge 
mobility  in  2-D  crystals  such  as  graphene  [48],  [49]  and 
M0S2  [46].  The  exact  mechanisms  likely  also  include 
phonons. 

At  this  time,  the  understanding  of  transport  in  2-D 
crystals  is  evolving.  It  is  clear  that  the  interactions  that 
limit  charge  transport  in  3-D  semiconductors  and  hete¬ 
rostructures  were  intrinsic  to  the  3-D  crystal  itself.  But  for 
2-D  crystals,  these  interactions  can  be  tuned  based  on 
what  we  put  around  them.  This  is  because  in  2-D  crystals 
we  have  direct  access  to  the  electrons,  their  spins,  and 
atomic  vibrations  to  an  unprecedented  degree.  As  our 
understanding  of  these  mechanisms  evolves,  the  level  of 
direct  access  to  the  physical  properties  may  well  prove  to 
be  the  defining  factor  that  differentiates  2-D  crystal  de¬ 
vices  from  their  3-D  counterparts.  This  feature  is  simul¬ 
taneously  an  advantage  and  a  challenge,  since  noise  and 
reliability  of  the  desired  nanoscale  devices  must  be  robust 
for  usability. 

C.  Possibility  of  2-D  Crystal  Heterostructures 

Heterostructures  based  on  3-D  crystals  take  advantage 
of  energy  band  offsets  that  originate  from  differences  in 
chemical  composition.  The  concept  of  quasi-electric  fields 
in  heterostructures  breaks  the  symmetry  of  electrical 
forces  acting  on  electrons  and  holes.  In  a  semiconductor  of 
constant  chemical  composition  (uniformly  doped,  or  p-n 
homojunctions),  the  electric  force  acting  on  electrons 
and  holes  is  the  same.  This  is  not  true  in  a  heterostructure 
[5].  This  broken  symmetry  is  central  to  quantum  confine¬ 
ment  and  high  oscillator  strengths  that  have  led  to  high- 
efficiency  LEDs  and  lasers.  QW  FETs  and  even  the  MOSFET 
gain  from  the  concept  of  quantum  confinement.  In  graded- 
base  heterostructure  bipolar  transistors  (HBTs),  the  broken 
symmetry  is  central  in  speeding  up  electrons  with  a  quasi¬ 
electric  field  in  the  same  region  in  space  where  there  is 
no  field  acting  on  holes  [50].  Examples  of  such  hetero¬ 
structures  based  on  3-D  crystal  semiconductors  include 
SiGe/Si,  AlGaAs/GaAs/InGaAs,  and  AlGaN/GaN/InGaN 
material  systems.  Except  in  special  cases,  most  of  such  3-D 
crystal  heterostructures  have  strain  due  to  the  lattice 
mismatch.  Strain  can  be  desirable  for  affecting  the  carrier 
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transport  or  as  the  driving  force  for  the  formation  of 
quantum  dots  by  the  Stranksi-Krastanov  mechanism  during 
epitaxy.  Strain  can  often  be  undesirable,  since  it  can  lead  to 
relaxation  and  defect  formation  beyond  certain  critical 
thicknesses. 

Heterostructures  based  on  2-D  crystals  are  at  their  in¬ 
fancy.  However,  a  number  of  interesting  features  are  likely 
to  emerge  in  them.  Initial  demonstrations  of  in-plane  2-D 
crystal  heterostructures  such  as  graphene  seamlessly  con¬ 
nected  to  BN  have  been  experimentally  observed,  and 
provide  exciting  opportunities  in  device  design  [51]. 
Hybrid  heterostructures  composed  of  2-D  crystals  such 
as  graphene  placed  on  3-D  semiconductors  such  as  silicon 
have  been  used  to  demonstrate  new  device  concepts.  One 
recent  example  is  a  graphene- Si  Schottky  diode  where 
graphene  may  be  thought  of  as  the  Schottky  “metal” 
contact.  However,  unlike  a  typical  metal,  the  Fermi  level 
of  graphene  can  be  tuned  with  a  third  gate  electrode, 
which  leads  to  a  variable  Schottky-barrier  height  [52].  This 
idea  was  used  to  demonstrate  a  variable-barrier  transistor 
(or  the  so-called  “Barristor”). 

Out-of-plane  or  vertical  heterostructures  are  also  real¬ 
ized  when  2-D  crystals  are  stacked  on  each  other.  Such 
heterostructures  do  not  suffer  from  lattice  mismatch  re¬ 
quirements,  since  there  are  no  interlayer  covalent  bonds. 
The  weak  van  der  Waal’s  interlayer  bonding  in  principle 
allows  unstrained  integration  of  2-D  crystal  layers  of  dif¬ 
ferent  material  properties.  One  may  envision  vertical 
heterostructures  of  2-D  crystal  metals,  semiconductors, 
insulators,  and  perhaps  a  wider  range  of  materials.  Due  to 
the  absence  of  broken  bonds,  the  interfaces  are  expected  to 
be  pristine  and  devoid  of  electronic  trap  states.  Interlayer 
transport  of  electrons  would  involve  tunneling.  The  rota¬ 
tional  alignment  of  the  2-D  crystal  layers  might  play  an 
important  role  in  such  heterostructures.  These  features  are 
currently  under  investigation,  and  are  certain  to  lead  to  a 
range  of  new  applications.  Initial  demonstrations  of  a 
graphene-BN-graphene  and  graphene-MoS2-graphene 
heterostructures  tunneling  transistors  have  been  recently 
reported  [53].  A  proposed  device  called  the  bilayer  pseudo¬ 
spin  FET  (BiSFET)  is  based  on  many-body  excitonic  con¬ 
densation  of  electron-hole  pairs  in  closely  spaced  layers  of 
graphene.  It  falls  under  the  category  of  vertical  2-D  crystal 
heterostructures  [54].  Its  single-particle  counterpart,  a 
tunneling  transistor  that  takes  advantage  of  the  symmetry 
of  the  bandstructure  of  some  2-D  crystals,  is  called  the 
“SymFET”  [55].  These  tunneling  devices  that  are  rooted  in 
2-D  crystals  are  described  in  Section  V. 


able  [56],  [57].  Chemical  vapor  deposition  (CVD)-grown 
graphene  has  been  realized  on  metals,  and  transferred  to 
other  substrates  [58].  Nanoribbons  have  been  fabricated 
on  CVD-grown  graphene  [59].  CVD-grown  graphene  has 
shown  promise  for  larger  area  crystals  than  epitaxial 
graphene,  which  is  limited  to  the  size  of  the  starting  3-D 
crystal  substrate.  The  crystal  quality  is  not  perfect  yet,  but 
as  was  the  case  in  the  development  of  3-D  crystals,  there  is 
reason  to  believe  it  will  undergo  drastic  improvements  in 
the  near  future. 

Similarly,  BN  2-D  crystals  have  been  grown  by  CVD,  as 
have  electronic-grade  M0S2  and  WS2  layered  materials 
[60]-[62].  However,  it  is  also  important  to  realize  that 
most  forms  of  2-D  crystals  have  been  produced  in  large 
volumes  in  their  layered  forms  [63].  They  have  already 
found  industrial  applications  in  chemical  catalysis  (M0S2, 
graphite),  lithium-ion  batteries  (lithium  cobaltate  and 
layered  carbon),  lubricants  (M0S2),  neutron  moderation  in 
nuclear  reactors  (graphite),  and  thermally  and  mechani¬ 
cally  refractory  crucibles  used  in  much  of  electronic 
material  and  device  processing  (BN  and  graphite).  The 
development  of  electronic  grade  counterparts  thus  is 
expected  to  heavily  leverage  the  considerable  prior  existing 
knowledge  and  industrial  base  for  these  materials. 

A  major  immediate  challenge  is  to  develop  methods  of 
doping  and  controlling  the  Fermi  level  in  2-D  crystals. 
Possible  methods  with  TMD  2-D  crystals  include  chemical 
substitutional  doping,  and/or  modulation  doping  by  taking 
advantage  of  the  rich  intercalation  chemistry  of  such  la¬ 
yered  materials.  Since  doping  control  is  intimately  con¬ 
nected  to  the  ability  to  form  low-resistance  contacts,  this 
challenge  assumes  increased  importance. 

The  development  of  electronic  grade  2-D  crystals  is 
expected  to  be  rapid.  The  first  active  device  applications 
are  expected  to  be  in  traditional  FETs.  For  example,  TMD- 
based  transistors  offer  attractive  routes  to  large-area  thin- 
film  transistors  (TFTs)  by  virtue  of  low  SS  values  and 
respectable  mobilities  when  compared  to  organic  semi¬ 
conductors  and  3-D  oxide  materials  [64].  But  can  they 
offer  new  functionalities  for  high-performance  devices  be¬ 
yond  what  is  being  envisioned  with  3-D  crystal  semi¬ 
conductors?  To  address  that  question,  we  focus  the  rest  of 
the  paper  on  one  of  the  possible  candidates  for  high- 
performance  and  low-power  energy-efficient  logic  devices: 
the  tunnel  FET  (TFET). 

III.  TUNNELING  TRANSPORT  IN 
SEMICONDUCTORS 


D.  Maturity  of  2-D  Crystals  and  Material  Challenges  Following  the  motivation  provided  earlier,  we  start  with  a 
Since  the  field  of  2-D  crystal  semiconductors  is  rela-  short  introduction  to  tunneling  transport  and  its  incor- 
tively  young,  a  short  discussion  of  the  material  challenges  poration  into  the  heart  of  the  transistor  operation.  The 
is  necessary.  Since  the  initial  demonstrations  in  2004,  the  discussion  starts  with  an  evaluation  of  the  effect  of  dimen- 
large-area  growth  capability  of  single-layer  graphene  has  sionality  on  interband  Zener  tunneling  [65],  [66]. 
expanded  rapidly  [21].  At  this  time,  epitaxial  single-layer  Consider  the  p  —  i  —  n  junction  shown  in  Fig.  4.  We 
graphene  on  several-inch-diameter  SiC  wafers  are  avail-  make  some  simplifying  assumptions  that  allow  us  to  zone 
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Fig.  4.  Interband  tunneling  in  a  reverse-biased  p  -  /  -  n  junction 
diode.  Most  TFETs  use  the  reverse-bias  Zener  tunneling  as  the 
mechanism  of  current  conduction  in  their  on- states.  The  current  may 
be  calculated  by  integrating  over  the  k-states  at  the  injection  point  as 
outlined  in  the  text. 


into  the  relevant  physics  immediately.  Assume  the  doping 
in  the  p-  and  n-sides  are  just  enough  to  align  the  Fermi 
levels  at  the  respective  band  edges.  Then,  under  no  bias, 
Ey  =  Enc  and  no  net  current  flows  across  the  junction. 
Under  the  application  of  a  reverse  bias  voltage  V,  a  finite 
energy  window  is  created  for  electrons  since  Ey  —  E^  =  qV. 
Within  this  energy  window,  electrons  from  the  valence 
band  can  tunnel  into  the  conduction  band  on  the  other 
side,  as  indicated. 

The  current  is  calculated  by  summing  the  individual 
contributions  by  each  k-state  electron.  There  are  many 
approaches  to  evaluate  currents,  but  none  is  as  transparent 
as  the  formalism  in  the  k-space.  To  illustrate,  we  write  the 
tunneling  current  as 


Then,  we  remove  the  bandgap  criteria  to  allow  for  special 
cases  such  as  graphene. 

The  tunneling  probability  is  obtained  by  the  Wentzel- 
Kramers-Brillouin  (WKB)  approximation  [67].  For  elec¬ 
trons  in  the  valence  band  of  the  p-side  with  transverse 
kinetic  energy  E^  =  h2k2±/2 m*,  the  WKB  tunneling  proba¬ 
bility  is  given  by  [68] 


7wKB  =  exp 


«  T0  exp 


E± 

E 


(2) 


where  Tq  =exp[—4y/2m^E^2/3qHF],  E  =  qhF / 2  ^J2 m^Eg , 
F  is  the  (constant)  electric  field  in  the  junction,  and 
is  the  reduced  effective  mass  given  by  = 
mX  /  (m*  +  mv)-  This  expression  is  found  to  be  consistent 
with  experimental  results  [69].  Note  that  the  tunneling 
probability  of  electrons  is  lowered  exponentially  with  their 
transverse  kinetic  energy.  To  evaluate  the  tunneling  cur¬ 
rent,  we  attach  this  tunneling  probability  to  each  electro¬ 
nic  k-state,  and  sum  it  over  all  electrons  incident  on  the 
tunneling  barrier. 

Three-dimensional  semiconductors:  Consider  the  case 
when  the  p  —  i  —  n  junction  is  made  of  3-D  crystal  semi¬ 
conductors.  In  Fig.  4,  we  concentrate  on  a  particular  2-D 
plane  as  shown  by  the  dashed  line,  at  the  p  —  i  junction. 
Half  of  the  electrons  in  the  valence  band  in  that  plane 
move  to  the  right  in  the  +kz  direction,  as  indicated  in  the 
hemisphere  in  the  k-space.  Since  there  are  negligible 
electrons  in  the  conduction  band  in  that  plane,  the  current 
there  must  be  carried  by  electrons  in  the  valence  band.  But 
which  of  these  right-going  electrons  are  allowed  to  tunnel 
through  the  gap?  In  the  absence  of  phonon  scattering, 
tunneling  is  an  elastic  process.  This  enforces  the  energy 
requirement 


El  — 


2m\ 


(kxp  +  kyp  +  kzp)  ~K  + 


2m* 


Xn+Kn  +  K 


n)  (3) 


iT  =  q^rJ2Xk)(f»-fc)T  (T 

L  k 

where  gs  =  2  is  the  spin  degeneracy  and  gv  is  the  valley 
degeneracy.  L  is  the  macroscopic  length  along  the  electric 
field  (which  will  cancel  out),  vg(k)  =  H~1VE(k )  is  the 
group  velocity  of  carriers  in  the  band  E(k),  fv,fc  are  the 
Fermi-Dirac  occupation  factors  of  the  valence  and  con¬ 
duction  bands,  and  T  is  the  tunneling  probability.  The  sum 
is  over  k-states  for  electrons  that  are  allowed  to  tunnel.  We 
illustrate  the  clarity  of  this  approach  by  using  the  same 
expression  for  evaluating  Zener  tunneling  currents  for 
p  —  i  —  n  junctions  made  of  3-D,  2-D,  and  1-D  crystals.  We 
first  consider  semiconducting  crystals  that  have  a  bandgap. 


with  the  additional  requirement  that  the  lateral  momen¬ 
tum  be  conserved.  To  simplify  the  analytical  treatment, 
and  in  preparation  for  2-D  crystals,  we  further  assume  that 
the  bands  are  symmetric,  i.e.,  m*  ~  m*  =  2 m£.  The  energy 
and  momentum  conservation  requirements  thus  lead  to 
the  relation 


2k2±+k 


2 

zp 


4KqV  ,2 

h2 


(4) 


where  k2±  =  k2p  +  k2p.  Let  us  define  k2mllx  =  4 m*RqV/H2. 
Since  there  is  an  electric  field  in  the  z-direction,  mo¬ 
mentum  in  that  direction  will  not  be  conserved.  For  the 
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electron  to  emerge  on  the  right  (n-)side,  kzn  must  be  non¬ 
zero,  and  thus  k2n  >  0,  which  implies 


2k 2  +  k2  <  k2 

zr<.±  -f  r<zp  \  Kmax* 


(5) 


The  above  condition  defines  a  restricted  volume  flj  of 
the  fe-space  hemisphere  for  electron  states  that  are  allowed 
to  tunnel.  We  are  now  in  a  position  to  evaluate  the  tun¬ 
neling  current  for  3-D  semiconductor  p  —  i  —  n  junctions. 
In  the  expression  for  the  tunneling  current  [see  (1)],  the 
group  velocity  term  is  that  of  the  valence  band  fe-state 
vg(k )  =  Hkz/m*.  We  skip  the  p-  or  n-subscripts,  since  it  is 
clear  that  the  electrons  tunnel  from  the  valence  band.  The 
expression  for  the  tunneling  current  is  then 


J  gsgv  f  n 

lT=£i“r  U  -v(/v-/c)Toex p 

z  (kx,k,,kz)enT 


Kkl 


2m*  E 


(6) 


The  sum  over  fe-states  is  converted  into  an  integral  via 
the  recipe  ^fe(. .  .)  — >  LxLyLz/ (2n)3  x  f  dkxdkydkz(. . .). 
To  evaluate  the  tunneling  current  in  the  restricted  volume, 
we  use  spherical  coordinates  (kx,ky,kz)  =  (ksin#cos  0, 
ksin#sin</>,  kcosO)  to  obtain  the  restricted  fe-space  vol¬ 
ume  fe2  <  femax/ (!  +  sin2#).  This  relation  is  representa¬ 
tive  of  the  “filtering”  brought  about  by  the  requirements  of 
energy  and  momentum  conservation.  Electrons  incident 
normal  to  the  junction  have  no  transverse  momentum.  For 
them  6  =  0,  and  they  are  allowed  to  tunnel.  The  number  of 
electron  states  allowed  to  tunnel  reduces  as  their  trans¬ 
verse  directed  momentum  increases.  The  current  carried 
by  these  states  with  transverse  momentum  is  further 
damped  by  the  exp[— Ej_/E]  factor,  leading  to  further  fil¬ 
tering  and  momentum  collimation. 

To  evaluate  the  current,  the  integral  in  fe-space  should 
be  evaluated.  To  simplify  the  evaluation  in  3-D  without 
losing  much  accuracy,  we  assume  fv  —  fc  «  1  for  the  energy 
window  of  current-carrying  electrons.  This  relation  is  exact 
at  0  K,  and  remains  an  excellent  approximation  even  at 
room  temperature.  The  tunneling  current  density  is  then 
given  by 
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Fig.  5.  Calculated  interband  tunneling  current  densities  in  a  few  3-D 
and  2-D  semiconductor  crystal  p  -  n  junctions .  The  left  figure  shows 
the  calculated  tunneling  current  densities  in  reverse-biased  p-n 
homojunctions.  If  the  current  per  unit  area  is  assumed  constant  for 
a  layer  thickness  of  10  nm,  then  the  effective  current  per  unit  width 
is  shown  in  the  right  axis  of  the  left  plot.  This  estimation  neglects 
quantization.  The  right  figure  shows  the  calculated  tunneling 
current  per  unit  widths  of  some  2-D  crystals.  The  transition  metal 
dichalcogenides  have  low  current  densities  due  to  high  bandgaps, 
whereas  2-D  graphene  has  the  highest  current  density. 
Two-dimensional  tunneling  currents  for  two  small  bandgap 
and  effective  masses  are  also  shown. 


where  the  fe-space  integral  is  evaluated  over  the  restricted 
volume  Qt-  The  units  are  in  current  per  unit  area  (A/cm2), 
as  it  should  be.  The  integral  yields  an  analytical  result. 
Using  the  symmetric  band  approximation  m*  ~  m*  =  2 m£, 
we  get 


t3-D 

JT 


6gsgvy/2^F 

8t T2h.2^/Eg 


To 


,V-2l{l-exp(-g)} 


(8) 


where  the  symbols  have  been  defined  earlier.  For  ex¬ 
tremely  small  reverse  bias  voltages  qV  2 E,  the  tunneling 
current  varies  as  J^"D  ~  V2  to  leading  order.  For  larger 
voltages  when  qV  2 E,  J^"D  ~  V  and  this  is  the  condition 
used  in  most  TFETs.  The  expression  for  the  tunneling 
current  shows  the  dependences  on  various  band  structure 
and  junction  parameters  explicitly. 

The  calculated  interband  tunneling  current  densities 
for  3-D  semiconductors  are  shown  in  Fig.  5(left)  for  a 
reverse  bias  voltage  of  0.3  V.  As  is  evident,  the  smaller 
bandgaps  of  InSb  and  InAs  favor  high  tunneling  current 
densities  that  approach  ~106  A/cm2.  If  we  assume  that 
the  body  thickness  of  the  p  —  i  —  n  junction  is  10  nm,  the 
effective  current  per  unit  width  is  also  shown  in  the  right 
axis  of  Fig.  5 (left).  However,  this  value  of  the  current  does 
not  account  for  the  increase  in  the  bandgap  due  to 
quantization,  which  we  address  shortly.  We  now  apply  the 
same  technique  for  calculating  tunneling  currents  in  2-D 
crystal  semiconductors. 
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Two-dimensional  semiconductors:  The  same  recipe  is 
repeated  for  2-D  crystals.  If  the  transport  is  along  the 
x-direction,  the  transverse  momentum  component  consists 
of  one  component  ky,  and  the  restricted  fc-space  volume  is 
given  by  2k2  +  k2  <  fe^ax.  The  interband  tunneling 
current  per  unit  width  in  a  2-D  crystal  p  —  i  —  n  junction 
then  evaluates  to 


t2-D 

JT 


Rgs& 


:»v/2  m£E 


27 r2R2 


Tn  x 


(qV-E)^Erf  A 


qV  •  2F  exp  -  -4=  (9) 


where  Erf  [. . .]  stands  for  the  error  function,  and  E  = 
qhF/2Ay2m^Eg  as  before.  For  extremely  small  reverse  bias 
voltages  qV  2 E,  the  tunneling  current  varies  as  J2_D  ~ 
V3/2  to  leading  order.  For  larger  voltages  when  qV  2 E, 
Erf  [. . .]  — >  1,  and  we  get  a  linear  dependence  of  the  tun¬ 
neling  current  on  the  reverse-bias  voltage  J2_D  ~ 
(q2gsgv  \J 27T mJF / 2tt2 ti2 )  T0 V.  We  note  that  the  units  are 
in  current  per  unit  width  (mA//xm),  as  should  be  the  case 
for  2-D  crystals.  In  quasi-2-D  systems,  multiple  subbands 
may  be  involved  in  transport.  Then,  we  sum  the  current 
from  each  subband  with  the  respective  band  parameters. 

For  the  special  case  of  2-D  graphene,  the  band  struc¬ 
ture  is  conical,  and  the  bandgap  is  zero.  The  interband 
tunneling  probability  for  a  graphene  in-plane  p  —  n  junc¬ 
tion  is  given  by  T(E,  6)  =  exp  [—7 rF2  sin2  0/qKvpF],  where 
0  is  the  angle  between  the  incident  electron  momentum 
and  the  junction  electric  field  F,  and  E  is  the  electron 
energy  [70].  The  requirement  of  lateral  momentum  con¬ 
servation  effectively  opens  a  bandgap  proportional  to  the 
lateral  momentum  of  electrons.  The  doping  in  the  p-  and 
n-graphene  regions  are  such  that  the  Fermi  level  to  Dirac 
point  energies  are  EFp  and  EFn,  respectively,  and  the  junc¬ 
tion  “depletion  width”  is  Lpn.  The  reverse-bias  tunneling 
current  in  the  2-D  graphene  p  —  n  junction  is  then  given 
by  [71] 


jCr  _  [El  +  EFn  ~  E 

T  7 r2H  y  hvFLpn 


(10) 


Let  us  assume  that  the  applied  reverse  bias  voltage  is 
small  compared  to  the  degeneracy  energies,  and  approx¬ 
imate  the  junction  field  by  qF  ~  ( EFp  +  EFn) /Lpn.  Then,  we 
obtain  an  approximate  expression  for  the  interband 
reverse-bias  tunneling  current  per  unit  width  in  2-D 
graphene  p  —  n  junctions  to  be  Jjr  ~  (q2  ^/qF / n2h^/hv^)V . 

The  interband  tunneling  current  densities  of  various 
2-D  crystals  are  plotted  in  Fig.  5 (right).  The  material  con¬ 
stants  (bandgaps  and  effective  masses)  are  obtained  from 


[39].  The  values  of  tunneling  currents  for  transition-metal 
dichalcogenides  are  low  owing  to  their  large  bandgaps. 
For  example,  the  current  density  approaches  ^0.1  fiA/fim 
for  MoTe2  at  a  high  field  of  4  MV/cm.  The  tunneling 
current  density  of  2-D  graphene  is  the  highest  (^several 
mA/pm),  but  it  lacks  a  bandgap.  As  new  2-D  crystals  come 
to  the  fore,  it  is  desirable  to  have  smaller  bandgaps  for 
boosting  the  current,  as  indicated  by  the  two  curves 
corresponding  to  hypothetical  2-D  crystals  with  bandgaps 
of  0.5  and  1.0  eV,  respectively.  Such  small-bandgap 
materials  could  be  intrinsic  2-D  crystals,  or  derived  from 
interaction-induced  bandgap  of  Dirac-cone  surface  states 
in  thin  topological  insulator  materials  [29].  Another  pos¬ 
sibility  is  in  bilayer  graphene,  where  breaking  the  layer 
symmetry  by  vertical  electric  fields  opens  a  small  bandgap 
[72]-[74].  It  is  clear  that  at  this  stage  the  currently 
available  TMD  family  of  2-D  crystal  semiconductors  can 
enable  tunneling  transistors.  But  for  in-plane  tunneling 
geometries,  the  current  densities  will  be  low.  This  feature 
can  be  effectively  addressed  by  either  narrower  gap  2-D 
crystal  semiconductors,  or  by  interlayer  tunneling  device 
geometries.  We  address  interlayer  tunneling  devices  in 
Section  V,  after  discussing  the  treatment  of  tunneling  in 
1-D  semiconductors. 

One -dimensional  semiconductors:  For  1-D  tunneling,  we 
obtain  an  exact  analytical  result  even  when  we  include  the 
Fermi-Dirac  occupation  factors  in  the  source  and  the  grain 
sides  of  the  p  —  i  —  n  junction.  In  the  ideal  1-D  case, 
E±  =  0  since  electrons  cannot  have  transverse  momen¬ 
tum.  When  a  voltage  V  is  applied,  fv  =  1/(1+  exp[(F  — 
qV)/feT])  and  fc  =  1/(1  +  exp[F/feT])  are  the  occupation 
functions  of  the  source  and  drain  sides.  The  interband 
tunneling  current  is  evaluated  by  the  same  prescription 
followed  for  the  3-D  and  2-D  cases  to  be  [75] 
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Note  the  explicit  appearance  of  the  Landauer  conduc¬ 
tance  in  the  expression.  This  expression  for  tunneling 
current  holds  for  quasi-l-D  semiconductors  such  as  semi¬ 
conducting  nanowires,  carbon  nanotubes,  or  semiconduct¬ 
ing  graphene  nanoribbons  (GNRs).  The  appropriate  WKB 
tunneling  probability  should  be  used.  For  nanowires  made 
from  conventional  3-D  semiconductor  crystals,  the  prob¬ 
ability  is  Tq  =  exp [—4^ 2mFE2J2 / 3 qhF]  as  before.  For 
CNTs  and  GNRs,  the  unconventional  band  structure  is 
captured  in  a  modified  WKB  tunneling  probability,  which 
is  given  by  T0  =  exp[— nE2 /4qhvFF],  where  vF  is  the  Fermi 
velocity  [75].  If  there  are  multiple  subbands  involved  in 
the  transport,  we  add  the  currents  from  each  subband  with 
the  right  bandgap. 

Fig.  6  shows  the  effect  of  quantization  on  bandgaps 
of  3-D  crystals  on  the  left,  and  the  calculated  1-D 
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Quantization  length  (nm) 


Fig.  6.  The  effect  of  quantization  on  the  bandgap  of  some  3-D  crystals. 
The  plot  is  generated  assuming  a  particle-in-a-box  quantization , 
and  is  meant  to  illustrate  the  approximate  trends.  The  effect  of 
quantization  is  the  most  severe  for  narrow  bandgap  semicondctors. 
The  increase  in  bandgap  will  reduce  interband  tunneling  currents. 
The  right  figure  shows  the  1-D  tunneling  currents  for  GNRs,  and  InSb 
and  Ge.  Note  the  large  reduction  of  current  due  to  quantization  effects 
in  Ge  and  especially  in  InSb.  A  major  advantage  of  2-D  crystals  is 
their  inherently  thin  nature.  In  addition,  their  large  effective  masses 
make  them  robust  to  quantization  effects  when  rendered  1-D. 


semiconductor  tunneling  current  densities  on  the  right. 
The  effect  of  quantization  is  to  increase  the  bandgap,  in 
turn  reducing  the  interband  tunneling  current.  The  right 
figure  shows  the  current  densities  of  “1-D”  semiconductors 
such  as  graphene  nanoribbons  (GNRs)  and  Ge  and  InSb 
nanowires.  The  tunneling  current  densities  of  GNRs  are 
the  highest  of  all  materials  calculated  that  possess  band- 
gaps.  The  effect  of  quantization  on  Ge  and  InSb  nanowire 
1-D  p  —  n  junction  structures  is  evident  from  the  preci¬ 
pitous  drop  in  the  interband  tunneling  currents  in  them. 

The  increase  in  the  bandgap  for  both  2-D  and  1-D 
confinement  is  calculated  using  a  simple  particle-in-a-box 
model  with  the  band-edge  effective  masses  of  the  3-D 
semiconductors.  The  values  are  meant  to  be  representative 
of  the  trends;  more  accurate  electronic  structure  calcula¬ 
tions  should  be  used  for  direct  validation.  However,  it  is 
clear  that  as  3-D  crystals  are  scaled  in  thickness  (for  mak¬ 
ing  them  2-D)  or  in  diameter  (for  making  them  1-D),  the 
corresponding  increase  in  bandgap  is  rapid.  Large  bandgap 
semiconductors  are  more  robust  to  quantization  since 
they  possess  heavier  effective  masses.  This  is  a  dilemma 
for  the  scaling  of  tunneling  transistors.  As  shown  in  the 
shaded  region  in  the  left  of  Fig.  6,  2-D  crystals  are 
typically  of  ^nanometer  thicknesses  and  span  bandgaps 
from  0  eV  (graphene)  to  several  eVs  (BN).  This  regime 
remains  inaccessible  to  3-D  crystal  semiconductors  due  to 
quantization.  It  is  possible  to  access  this  regime  with  3-D 
semiconductors  only  if  the  band  structure  allows  for 
extreme  anisotropies  [76],  but  such  highly  desirable 
properties  are  yet  to  be  demonstrated  in  3-D  crystal 
semiconductors. 


IV.  TUNNELING  TRANSISTORS  WITH 
3-D  CRYSTALS 

The  unified  view  of  tunneling  transport  discussed  in  the 
last  section  provides  a  framework  for  comparative  studies 
of  the  effect  of  dimensionality  on  tunneling  transistors. 
Based  on  the  discussion  of  transport  in  two-terminal  tun¬ 
nel  junctions,  we  now  discuss  the  electrostatics  and  device 
embodiments  of  the  corresponding  three-terminal  TFETs. 

In  a  TFET,  a  gate  terminal  electrostatically  controls  the 
energy-band  alignment  of  the  p  —  n  junction,  as  indicated 
in  Fig.  7.  In  the  OFF-state  of  the  device,  electrons  in  the 
valence  band  of  the  source  are  energetically  forbidden  to 
tunnel  to  the  drain  since  the  channel  length  exponentially 
damps  the  direct  source-to-drain  tunneling  probability.  To 
turn  the  device  on,  the  gate  pushes  the  channel  bands  to 
align  the  conduction  band  edge  of  the  channel  region  with 
the  valence  band  of  the  source.  Electrons  can  now  tunnel 
through  the  tunneling  barrier,  which  is  much  smaller  than 
the  OFF-state.  The  goal  therefore  is  to  allow  a  large  current 
to  flow  in  the  ON-state,  while  cutting  the  current  off  as 
much  as  possible  in  the  OFF-state. 

The  performance  requirements  of  a  TFET  are  indicated 
schematically  in  Fig.  7.  Compared  to  a  MOSFET,  the 
steeper  SS  slope  of  a  TFET  enables  a  higher  ON-current  at  a 
smaller  gate  overdrive  voltage.  This  feature  is  expected  to 
enable  scaling  of  the  voltage  supply  Vdd  to  lower  values 
while  maintaining  a  substantial  on/off  ratio.  The  issues  of 
electrostatics  and  transport  have  been  discussed  at  length 
in  various  articles  [12],  [69].  We  refer  the  reader  to  these 
articles  for  detailed  historical  perspectives  and  further 
technical  details.  Here,  we  qualitatively  discuss  a  few  em¬ 
bodiments  and  issues  with  TFETs  realized  with  3-D  crystal 
semiconductors.  The  discussion  naturally  motivates  the 
case  for  2-D  crystal  realizations  of  the  device. 

The  electric  field  lines  emanating  from  the  gate  metal 
of  a  TFET  need  to  access  the  p  —  n  junction.  Therein  lies  a 
dilemma  for  TFETs  based  on  3-D  crystal  semiconductors. 
As  shown  in  Fig.  8,  if  the  tunneling  current  flows  in  the 
lateral  direction  and  the  gate  field  is  vertical,  the  channel 
needs  to  be  thinned  down  to  exercise  substantial  electro¬ 
static  control  over  the  entire  junction  thickness.  As  the 


Fig.  7.  TFET  operation  and  requirements.  The  left  figure  shows  the 
OFF-state  energy-band  diagram  of  the  TFET  along  the  tunneling 
direction.  The  right  figure  shows  the  ON-state.  The  channel  band  is 
controlled  by  the  gate.  TFETs  are  expected  to  lower  the  supply  voltage 
VDD  since  a  steeper  SS  swing  leads  to  a  higher  ON-current  at  a  smaller 
voltage,  as  shown  in  the  middle. 
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Fig.  8.  Schematic  representation  of  various  topologies  of  TFETs.  The  top  row  shows  TFETs  where  the  tunneling  current  flows  laterally  in  the 
p  n  -  n+  junction.  The  circles  are  indications  of  the  region  in  space  where  most  of  the  interband  tunneling  current  flows.  Since  the  gate  is 
on  the  top,  parts  of  the  junction  farther  away  from  it  are  not  effectively  gated  in  the  top  left  TFET.  The  top  middle  geometry  is  the  same  as  the  left, 
but  with  a  thinner  channel  for  more  uniform  electrostatic  gate  control  of  tunneling  current.  The  right  figure  on  the  top  row  is  the  2-D  crystal 
realization  of  the  lateral  TFET.  As  the  channel  thickness  is  reduced  in  3-D  semiconductors,  quantum  confinement  increases  the  bandgap  and 
reduces  the  tunneling  current.  This  is  avoided  in  2-D  crystals.  To  increase  the  net  current,  vertical  TFETs  are  being  considered.  The  bottom  row 
indicates  some  realizations  of  TFETs  in  which  the  tunneling  current  flows  vertically.  The  left  figure  shows  a  side-gate  geometry,  and  the  middle 
figure  is  a  geometry  in  which  the  current  flow  is  not  over  a  “line,”  but  an  “area,”  as  shown  by  the  shaded  ellipse.  The  right  figure  shows  the 
realization  of  a  vertical  double-gate  TFET  with  p+  and  n+  2-D  crystal  layers.  It  highlights  the  electrostatic  advantage  and  simplicity. 


channel  thickness  is  scaled  down,  quantum  confinement 
increases  the  bandgap,  and  thus  the  interband  tunneling 
current  reduces  (see  Fig.  6).  A  2-D  crystal  does  not  suffer 
from  such  a  problem,  and  thus  offers  a  way  to  fight  quan¬ 
tization  effects.  In  addition,  it  offers  a  solution  to  surface 
state  related  trap  states,  and  simpler  integration  of  double¬ 
gate  geometries,  as  indicated  in  Fig.  8. 

A  number  of  TFETs  with  subthreshold  slopes  less  than 
the  SS  limit  of  60  mV/decade  have  been  demonstrated, 
proving  the  feasibility  of  the  concept.  Such  devices  have 
been  made  with  3-D  crystal  semiconductors  (Si,  Ge,  etc.) 
as  well  as  with  carbon  nanotubes  [77]  — [80] .  However,  for 
most  realizations,  the  ON-state  current  falls  below  the 
mAI/im  range  necessary  for  high-performance  ope¬ 
ration.  Low  ON-current  TFETs  can  enable  various  new 
applications  where  performance  (speed)  requirements  are 
not  as  critical  as  the  requirement  of  low  power  consump¬ 
tion.  For  high-performance  TFETs,  various  approaches  are 
being  pursued  to  increase  the  ON-current.  These  ap¬ 
proaches  involve  using  heterojunctions  that  have  staggered 
or  broken-gap  band  alignments,  or  through  changes  in  the 
device  topology.  An  approach  based  on  the  device  topology 
is  indicated  in  Fig.  8. 

The  shaded  regions  in  Fig.  8  indicate  the  location  of 
current  flow.  To  increase  the  tunneling  current  per  unit 
width,  it  is  necessary  to  increase  the  net  area  of  tunneling 
current  flow.  The  vertical  geometries  shown  in  Fig.  8  allow 
this  change  [81],  [82].  The  gate  field  effect  is  in  the  same 
direction  as  the  tunneling  current  flow  in  such  devices. 
The  tunneling  current  follows  a  nonlinear  path  (shaped 
like  an  “S”)  laterally  from  the  source,  vertically  into  the 
drain,  and  then  out  laterally  into  the  drain.  The  device 
geometry  requires  careful  processing.  For  this  geometry, 
two  layers  of  2-D  crystals,  one  doped  p-type  and  the  other 


n-type,  promise  efficient  vertical  scaling  and  electrostatic 
control  as  shown  in  the  figure.  It  may  also  enable  a  sim¬ 
plification  of  the  processing  requirements. 

V.  TUNNELING  TRANSISTORS  WITH 
2-D  CRYSTALS 

The  2-D  crystal  realizations  of  TFETs  discussed  here  in¬ 
volve  in-plane  tunneling  for  the  lateral  device  and  interlayer 
tunneling  in  the  vertical  TFET.  We  discuss  them  in  greater 
detail  here.  Note  that  due  to  the  relatively  early  phase  of 
material  development,  we  estimate  and  project  the  perfor¬ 
mance  advantages  in  cases  where  experimental  results  are 
not  available  yet. 

A.  In-Plane  Tunneling:  2-D  Crystal  Semiconductors 

The  in-plane  interband  tunneling  currents  calculated 
in  Fig.  5  show  that  smaller  bandgap  semiconductor  2-D 
crystals  are  required  for  boosting  the  ON-state  current  of 
the  devices.  Low-power  TFETs  are  realizable  with  the 
transition-metal  dichalcogenide  semiconductors.  The  ef¬ 
fective  masses  of  the  conduction  and  valence  band  edges  of 
TMD  2-D  crystals  have  been  calculated  to  be  rather  sym¬ 
metric.  For  example,  the  electron  effective  mass  of  M0S2  is 
~0.5 7,  and  the  hole  effective  mass  is  ~0.66  [39].  The 
symmetry  in  the  band  structure  is  expected  to  lead  to 
symmetric  performance  of  nTFETs  and  pTFETs,  which 
would  be  essential  for  complementary  logic  circuits. 

It  has  been  found  that  multilayer  versions  of  TMD  2-D 
crystals  have  smaller  bandgaps  than  the  single-layer 
counterpart,  and  are  generally  of  indirect  bandgap  nature 
[40],  [41].  This  is  also  true  when  one  considers  single-layer 
graphene  (direct  bandgap)  and  graphite  (which  is  a  semi¬ 
metal).  Furthermore,  it  has  recently  been  reported  that 
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carrier  inversion  can  be  achieved  in  multilayer  TMD 
crystals  by  the  field  effect.  A  hole  channel  was  observed  in 
a  nominally  n-type  layered  semiconductor  [64].  Consider  a 
few-layer  stack  of  2-D  TMD  crystals.  By  using  two  gates,  it 
is  possible  to  create  an  electron  channel  at  one  interface 
and  a  hole  channel  in  the  other.  These  channels  can  be 
placed  several  nanometers  apart  by  controlling  the 
number  of  layers.  The  wave  function  overlap  between 
these  states  is  small  at  no  bias  owing  to  the  high  effective 
mass  for  carrier  motion  between  planes.  The  geometry 
then  allows  for  a  TFET  similar  to  the  vertical  structure 
shown  in  Fig.  8,  but  without  the  need  to  chemically  dope 
the  individual  layers.  A  major  challenge  in  such  structures 
is  in  the  formation  of  ohmic  contacts  to  the  individual 
layers.  Note  that  such  a  device  has  also  been  recently 
proposed  for  thin  layer  Si  [83].  The  realization  with 
multilayer  version  of  2-D  crystals  can  be  an  alternative 
approach  that  can  leverage  the  robustness  against  quan¬ 
tization  effects,  and  relative  insensitivity  to  surface  and 
interface  trap  effects. 

B.  In-Plane  Tunneling:  2-D  Graphene 

As  shown  in  Fig.  5,  the  ON-state  interband  tunneling 
current  density  in  2-D  graphene  is  the  highest  due  to  the 
absence  of  a  bandgap.  For  the  same  reason,  it  is  difficult  to 
obtain  the  OFF-state  condition  using  monolayer  2-D 
graphene.  Field-tunable  bandgaps  in  bilayer  graphene 
have  been  proposed  as  a  possible  approach  to  achieving 
on/off  ratios  in  TFETs  [84].  There  have  also  been  recent 
reports  of  the  observation  of  negative  differential  resis¬ 


tance  in  monolayer  2-D  graphene  FETs  [85].  The  proposed 
mechanism  responsible  for  such  behavior  relies  entirely  on 
gate  electrostatics  and  the  unique  band  structure  with  the 
zero-gap  nature  of  2-D  graphene.  More  experimental  work 
and  understanding  of  NDR  mechanisms  in  2-D  graphene 
can  lead  to  useful  device  applications  in  the  analog  arena  to 
complement  TFETs.  To  decrease  the  OFF-state  current  for 
in-plane  tunneling  devices,  it  is  necessary  to  create 
bandgaps  in  graphene.  One  approach  is  to  use  CNTs  or 
lithographically  patterned  GNRs,  which  is  discussed  next. 

C.  In-Plane  Tunneling:  CNTs  and  GNRs 

One  of  the  early  reports  of  sub-60-mV/decade  SS  slope 
TFET  behavior  was  observed  in  semiconducting  carbon 
nanotubes  at  room  temperature  [80].  Analysis  of  the  de¬ 
vice  performance  for  CNT  TFETs  [86]  and  GNR  TFETs 
[87]  shows  that  they  are  attractive  for  desirable  on- 
currents,  on/off  ratios,  and  sub-60-mV/decade  SS  slopes. 
CNTs  do  not  have  edge  states,  and  are  the  most  attractive 
from  a  performance  viewpoint.  Bandgap  control,  chemical 
doping,  and  patterned  assembly  on  large  wafers  still  re¬ 
main  challenging  for  CNTs,  though  rapid  progress  is  being 
made  [88]. 

Their  close  cousins,  GNRs  are  also  highly  attractive 
candidates  for  TFETs.  For  example,  Fig.  9  shows  the  device 
structure,  energy  band  diagrams,  and  the  projected  device 
characteristics  of  complementary  GNR  TFETs.  The  inclu¬ 
sion  of  parasitic  elements  to  the  intrinsic  model  still 
maintains  a  high  performance.  GNRs  can  be  integrated 
on  planar  surfaces,  and  can  be  made  lithographically  in 
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Fig.  9.  A  proposed  GNR  TFET  geometry,  energy  band  diagram,  and  the  calculated  transistor  transfer  curves.  The  device  structure  consists  of 
a  CNR  p  -  n  junction  that  is  gated  through  an  insulator  from  the  top  gate.  The  energy  band  diagrams  are  for  a  20-nm-long  channel  device 
with  a  5-nm-wide  CNR.  The  energy  band  diagrams  indicate  the  off-  and  ON-states  of  the  device,  where  the  channel  potential  is  moved  with 
the  gate  voltage.  The  resulting  transfer  curve  shows  a  high  ON-current,  a  low  OFF-current ,  and  a  lowSS  slope,  below  the  60-mV/decade  limit. 
Though  the  calculations  are  for  an  ideal  case,  they  represent  the  attractiveness  of  GNRs  as  possible  candidates  for  TFETs.  The  figure  has 
been  adapted  from  [87]. 
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parallel  arrays  to  boost  the  net  current  in  a  realistic  TFET 
device  geometry.  The  available  dangling  bonds  at  the  edges 
can  be  used  to  chemically  dope  them;  initial  reports  indi¬ 
cate  this  possibility  [89].  The  major  challenges  at  this  stage 
for  the  realization  of  GNR  TFETs  lie  in  the  narrowness  of 
the  GNR  widths  necessary  to  avail  high-performance 
levels.  The  energy  bandgap  of  a  semiconducting  GNR  of 
width  W  is  Eg  ~  1.4/W  eV,  where  W  is  in  nanometers. 
Based  on  theoretical  estimates,  GNRs  of  widths  <  10  nm 
are  necessary.  The  line-edge  roughness  that  might  result 
from  process  variations  for  the  thinnest  GNRs  can  degrade 
the  performance  of  GNR  TFETs,  as  has  been  analyzed  in 
[90].  On  the  other  hand,  advances  in  process  control  in  the 
fabrication  of  thin  films  in  Si  FinFETs  can  be  effectively 
leveraged  for  fabrication  of  wafer-scale  GNRs.  A  number  of 
variants  of  the  GNR  TFETs  have  also  been  proposed  to 
improve  the  device  performance  [91]-[95]. 

The  symmetry  of  the  band  structure  of  CNTs  and  GNRs 
is  a  major  advantage  that  allows  for  the  realization  of 
nTFETs  and  pTFETs  on  equal  footing.  Combined  with  the 
scaling  advantages  that  stem  from  their  atomically  thin 
body  nature,  they  are  highly  desirable  for  nanoscale 
TFETs.  The  approach  to  high-performance  TFETs  using 
3-D  crystal  semiconductors  is  taking  the  path  toward 
materials  with  successively  smaller  bandgaps  to  increase 
the  ON-current.  The  approach  with  graphene,  CNTs,  and 
GNRs  is  from  the  other  extreme,  where  we  start  from  zero 
bandgap  and  very  high  ON-currents,  and  now  need  to  open 
bandgaps  controllably  to  lower  the  OFF-current.  While  this 
is  an  attractive  and  complementary  approach,  2-D  crystals 
also  offer  the  possibility  of  interlayer  tunneling  transistors, 
which  we  discuss  now. 

D.  Interlayer  Tunneling  Devices,  BiSFETs, 
and  SymFETs 

Electron  tunneling  out  of  the  plane  of  a  2-D  crystal  is 
under  intense  scrutiny  at  this  time  [96].  The  electronic 
band  structure  of  the  2-D  crystal  is  defined  in  the  plane  but 
not  out  of  it.  The  conventional  approach  to  tunneling 
calculations  requires  the  knowledge  of  band  parameters 
such  as  the  effective  mass  of  the  evanescent  band  structure 
in  the  direction  of  the  tunneling.  Since  this  feature  is  not 
well  defined  for  2-D  crystals,  it  is  more  feasible  to  use 
scattering  rate  formalisms  for  quantitative  calculations  of 
interlayer  tunneling.  The  Bardeen  transfer-Hamiltonian 
approach,  used  in  scanning  tunneling  microscopy  [97], 
[98]  and  in  superconducting  Josephson  junctions  [99] 
allows  such  evaluation.  We  do  not  derive  the  quantitative 
results  here,  but  refer  the  reader  to  recent  articles  that 
approach  the  subject  of  interlayer  tunneling  using  the 
Bardeen  method. 

A  prototype  interlayer-tunneling  device  is  a  graphene- 
insulator-graphene  (GIG)  junction.  In  a  recent  work 
[100],  the  interlayer  tunneling  current  in  such  a  GIG 
junction  was  explicitly  evaluated  using  the  Bardeen  meth¬ 
od.  The  predicted  I-V  characteristics  are  rather  remark- 


Fig.  10.  Band  alignments  of  GIG  interlayer  tunnel  junctions  under 
various  bias  conditions  from  [lOO].  The  graphene  layers  are 
doped  to  form  ap-n  junction .  In  (a)  and  (b),  the  symmetry  of 
the  band  structure  restricts  electrons  at  only  one  energy  to  carry 
interlayer  current  due  to  the  requirement  of  transverse  momentum 
conservation.  A  special  case  occurs  when  the  Dirac  points  align: 
electrons  at  all  energies  are  now  allowed  to  tunnel ,  leading  to 
a  spike  in  the  current ,  as  shown  schematically  in  (d). 


able,  and  highlight  the  strong  role  of  the  symmetry  of  the 
band  structure  of  graphene. 

Fig.  10  shows  the  energy  band  alignments  and  pro¬ 
jected  device  performance  of  a  GIG  interlayer  tunnel 
junction  device.  The  two  graphene  layers  are  “indepen¬ 
dent”  in  the  sense  that  they  do  not  form  a  bilayer,  and  they 
are  doped  p-  and  n-type  as  captured  by  their  Fermi  level 
degeneracies.  Ohmic  contacts  are  made  to  the  two  layers 
independently.  A  voltage  is  applied  across  the  junction. 
When  the  Dirac  points  of  the  two  layers  are  misaligned,  a 
small  interlayer  tunneling  current  flows.  The  circles  indi¬ 
cated  on  the  Dirac  cones  in  Fig.  10(a)  and  (b)  show  the 
states  that  participate  in  the  interlayer  tunneling  process. 
Electrons  that  have  energy  halfway  between  the  Dirac 
points  carry  the  current.  This  is  because  transverse  mo¬ 
mentum  conservation  requires  the  radii  of  the  iso-energy 
circles  to  be  the  same  in  both  layers. 

However,  at  the  particular  voltage  when  the  Dirac 
points  align,  as  shown  in  Fig.  10(c),  electrons  at  all  ener¬ 
gies  are  now  allowed  to  tunnel,  leading  to  a  large  spike  in 
the  current.  This  is  schematically  shown  in  Fig.  10(d)  as  a 
Dirac-delta  function.  A  quantitative  evaluation  leads  to 
broadening,  but  with  a  very  large  NDR  effect.  Note  that 
the  peak  would  be  much  smaller  if  the  band  structure  was 
not  symmetric.  Since  then  the  requirement  of  transverse 
momentum  conservation  would  restrict  the  current  to  flow 
at  a  particular  energy,  and  a  collective  tunneling  condition 
as  in  Fig.  10(c)  cannot  be  achieved.  The  large  tunneling 
current  peak  is  a  direct  consequence  of  the  symmetric 
band  structure  of  2-D  graphene. 
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The  GIG  p  —  n  junction  structure  can  be  connected  to 
gates  to  realize  an  interlayer-tunneling  transistor.  Such  a 
device,  called  the  symmetric-FET  (SymFET)  has  been  re¬ 
cently  proposed  [55].  In  addition  to  performing  logic  ope¬ 
rations,  the  inherently  fast  tunneling  feature  and  large 
NDR  promises  to  also  enable  analog  applications  such  as 
high-harmonic  generation,  and  high-speed  oscillator  de¬ 
sign.  Note  that  the  SymFET  device  structure  is  similar  to 
gated  RTD  structures  [101]  realized  in  3-D  semiconductor 
heterostructures,  but  takes  advantage  of  the  band  structure 
symmetry  of  graphene  to  deliver  a  stronger  NDR  behavior. 
The  first  experimental  report  of  such  a  structure  did  not 
show  NDR,  but  exhibited  TFET-like  behavior  with  a  few 
orders  on/off  ratio  at  room  temperature.  The  structure 
used  consisted  of  graphene-BN-graphene  and  graphene- 
MoS2-graphene  heterostructures  [53]. 

The  SymFET  structure  is  based  on  single-particle 
tunneling.  Realistic  fabrication  of  the  device  calls  for 
rotational  alignment  of  the  graphene  layers.  By  adjusting 
the  interlayer  distance  and  the  carrier  densities,  the 
Coulombic  forces  between  the  electrons  and  the  holes  in 
the  two  graphene  layers  can  be  made  strong  enough  to  form 
excitonic  quasi-particles.  Under  suitable  bias  conditions, 
the  interlayer  current  flow  can  take  a  collective  many-body 
form  triggered  by  a  Bose-Einstein  condensation  of  the 
excitons.  The  condensate  can  boost  the  interlayer  current 
significantly.  The  proposed  device,  called  the  bilayer 
pseudospin  FET  (BiSFET)  is  insensitive  to  the  rotational 
alignment  of  the  two  graphene  layers.  It  has  been  shown 
that  if  the  BISFET  can  be  realized,  it  can  perform  digital 
logic  by  consuming  many  orders  of  magnitude  lower  energy 
than  conventional  MOSFETs  [54].  The  SymFET  and  the 
BiSFET  are  fundamentally  new  types  of  devices  with  no 
direct  analogs  to  conventional  semiconductors.  This  is 
because  of  their  unique  band  structures  and  their  2-D 


crystal  nature.  Their  discussion  is  an  ideal  point  to  end  this 
review  paper  and  to  wrap  up  with  a  few  concluding 
remarks. 

VI.  FUTURE  PERSPECTIVES  AND 
CONCLUSION 

The  emergence  of  2-D  crystal  materials  has  marked  a  new 
phase  for  the  development  of  semiconductor  devices.  It 
may  rank  at  the  same  level  as  the  origin  and  proliferation 
of  heterostructures  in  3-D  semiconductors.  The  materials 
and  the  resulting  devices  are  at  their  infancy,  as  are  many 
device  ideas  based  on  tunneling  that  are  at  proposal  stages. 
But  the  novelty  the  family  of  2-D  crystal  has  brought  to  the 
field  becomes  evident  by  the  string  of  new  device  concepts 
based  on  tunneling.  The  addition  of  graphene  with  its 
unique  band  structure,  BN  as  a  2-D  crystal  insulator,  and 
transition-metal  dichalcogenides  with  material  properties 
ranging  from  semiconducting  to  metallic  and  super¬ 
conducting  casts  a  much  wider  net  than  has  been  possible 
with  conventional  materials.  The  possibility  of  integration 
of  diverse  material  properties  in  2-D  crystal  heterostruc¬ 
tures  has  breathed  new  life  into  existing  paradigms  of 
electronic  device  technologies.  This  is  an  exciting  time 
when  creative  ideas  are  needed  to  exploit  the  power  of  this 
new  material  system.  Though  it  is  impossible  to  predict  the 
exact  path  forward,  we  can  be  sure  that  electronic  devices 
that  go  far  beyond  the  current  state  of  the  art  will  result 
from  the  new  material  family.  ■ 
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The  electron  transport  properties  of  atomically  thin  semiconductors  such  as  MoS2  have  attracted 
significant  recent  scrutiny  and  controversy.  In  this  work,  the  scattering  mechanisms  responsible  for  limiting 
the  mobility  of  single-layer  semiconductors  are  evaluated.  The  roles  of  individual  scattering  rates  are 
tracked  as  the  two-dimensional  electron  gas  density  is  varied  over  orders  of  magnitude  at  various 
temperatures.  From  a  comparative  study  of  the  individual  scattering  mechanisms,  we  conclude  that  ah 
current  reported  values  of  mobilities  in  atomically  thin  transition-metal  dichalcogenide  semiconductors  are 
limited  by  ionized  impurity  scattering.  When  the  charged  impurity  densities  are  reduced,  remote  optical 
phonon  scattering  will  determine  the  ceiling  of  the  highest  mobilities  attainable  in  these  ultrathin  materials 
at  room  temperature.  The  intrinsic  mobilities  will  be  accessible  only  in  clean  suspended  layers,  as  is  also  the 
case  for  graphene.  Based  on  the  study,  we  identify  the  best  choices  for  surrounding  dielectrics  that  will  help 
attain  the  highest  mobilities. 
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Two-dimensional  (2D)  layered  crystals  such  as  single 
layers  of  transition-metal  dichalcogenides  represent  the 
thinnest  possible  manifestations  of  semiconductor  materi¬ 
als  that  exhibit  an  energy  band  gap.  For  example,  a  single¬ 
layer  (SL)  MoS2  is  around  ~0.6  nm  thick  and  exhibits  an 
energy  band  gap  of  around  ~1.8  eV  [1].  Such  semicon¬ 
ductor  layers  differ  fundamentally  from  ultrathin  hetero¬ 
structure  quantum  wells  or  thin  membranes  carved  out  of 
three-dimensional  (3D)  semiconductor  materials  because 
there  are,  in  principle,  no  broken  bonds,  and  no  roughness 
over  the  2D  plane.  In  heterostructure  quantum  wells,  the 
electron  mobility  suffers  from  variations  in  the  quantum- 
well  thickness.  A  classic  “sixth-power  law”  from  Sakaki 
et  al.  [2]  shows  that  since  the  quantum-mechanical  energy 
eigenvalues  in  a  heterostructure  quantum  well  of  thickness 
L  go  as  e~l/L2,  variations  in  thickness  A L  lead  to 
perturbations  of  the  energy  Ae  ~  —2AL/L3.  Since  the 
scattering  rate  depends  on  the  square  of  Ae ,  the  rough¬ 
ness-limited  mobility  degrades  as  jiR  ~  L6.  When  L 
reduces  from  about  ~7  to  ~5  nm  for  example,  pR  reduces 
from  about  104  to  103  cm2/ Vs  in  GaAs/AlAs  quantum 
wells  at  4.2  K  [2].  Though  low-temperature  mobilities 
exceeding  106  cm2/  Vs  have  been  achieved  in  such  hetero¬ 
structures  by  scrupulous  cleanliness  and  design  to  reduce 
roughness  scattering,  the  statistical  variations  in  the  quan¬ 
tum-well  thickness  during  the  epitaxial  growth  process 
pose  a  fundamental  limit  to  electron  mobility. 
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Because  of  the  absence  of  intrinsic  roughness  in  atomi¬ 
cally  thin  semiconductors,  the  expectation  is  that  higher 
mobilities  should,  in  principle,  be  attainable.  However, 
recent  measurements  in  MoS2  and  similar  semiconductors 
[3-5]  exhibit  rather  low  mobilities  in  single  layers,  which 
are,  in  fact,  lower  than  in  their  multilayer  counterparts. 
Many-particle  transport  effects  can  appear  in  transition- 
metal  dichalcogenides  under  special  conditions  because 
of  the  contribution  of  highly  localized  d-orbitals  to  the 
conduction  and  valence-band-edge  eigenstates.  Collective 
effects  have  been  observed  in  multilayer  structures,  such  as 
charge-density  waves  [6,7]  and  the  appearance  of  super¬ 
conductivity  at  extremely  high  metallic  carrier  densities  [8] 
under  extreme  conditions.  We  do  not  discuss  such  collec¬ 
tive  phenomena  here.  Instead,  we  focus  on  single-particle 
transport  in  single-layer  MoS2 ;  the  only  many-particle  effect 
included  is  free-carrier  screening.  In  this  work,  we  perform 
a  comprehensive  study  of  the  scattering  mechanisms  that 
limit  electron  mobility  in  atomically  thin  semiconductors. 
The  mobility  is  calculated  in  the  relaxation-time  approxi¬ 
mation  (RTA)  of  the  Boltzmann  transport  equation.  The 
results  shed  light  on  the  experimentally  achievable  electron 
mobility  by  designing  the  surrounding  dielectrics  and  low¬ 
ering  the  impurity  density.  The  findings  thus  offer  useful 
guidelines  for  future  experiments. 

With  the  advent  of  graphene,  it  was  realized  that  for 
ultrathin  semiconductors,  the  dielectric  environment  plays  a 
crucial  role  in  electron  transport.  It  has  now  been  demon¬ 
strated  that  the  dielectric  mismatch  significantly  modifies 
the  Coulomb  potentials  inside  a  semiconductor  thin  layer 
[9-12].  Electrons  in  the  semiconductor  can  also  remotely 
excite  polar-optical-phonon  modes  in  the  dielectrics 
[13-19].  Such  long-range  interactions  become  stronger  as 
the  thickness  of  the  semiconductor  layer  decreases.  Thus, 
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FIG.  1.  Coulomb  potential  contours  due  to  an  on-center  point  charge  for  three  different  dielectric  environments:  ee  —  1,  7.6  (=£*),  100. 


one  can  expect  the  dielectric  environment  to  significantly 
affect  electron  transport  properties  in  SL  gapped  semi¬ 
conductors.  In  this  work,  we  take  SL  MoS2  as  a  case  study 
to  investigate  such  effects.  The  results  and  conclusions  can 
be  extended  to  other  SL  gapped  semiconductors. 

We  first  study  the  effect  of  the  dielectric  environment 
on  Coulomb  scattering  of  carriers  from  charged  impurities 
located  inside  the  MoS2  single  layer.  Figure  1(a)  shows  a 
point  charge  located  at  the  center  (z0  =  0)  of  a  SL  MoS2  of 
thickness  a.  Assuming  the  surrounding  dielectric  provides 
a  large  energy  barrier  for  confining  electrons  in  the  MoS2 
membrane,  we  consider  scattering  of  electrons  within  the 
conduction  band  minima  at  the  K  point,  i.e.,  in  the  ground 
state.  The  envelope  function  of  mobile  electrons  is  then 

ip-ip,  z)  =  x(z)elk'p /VS,  where j(z)  =  Jlfaco^inzl a), 

k 

S  is  the  2D  area,  k  is  the  in-plane  2D  wave  vector,  and  p  is 
the  in-plane  location  vector  of  the  electron  from  the  point 
charge.  The  dielectric  mismatch  between  the  MoS2  (rela¬ 
tive  dielectric  constant  es )  layer  and  its  environment  ( ee ) 
creates  an  infinite  array  of  image  charges  at  points  zn  =  na, 
where  n  =  ±1,  ±2. . .  [9,10,20].  The  nth  point  charge  has  a 
magnitude  of  ey W,  where  y  =  (es  —  ee)/(es  +  ee).  These 
image  charges  contribute  to  the  net  electric  potential  seen 
by  the  electron,  which  is  given  by 


ci  cy\n\ 

UnSc('/?’  ^  „=-oo4tf£o£svV  +  k  -  z«|2 


a) 


where  e  is  the  elementary  charge,  and  £0  is  the  vacuum 
permittivity.  Figure  1  shows  the  net  unscreened  Coulomb 
potential  contours  in  the  dielectric/MoS2 /dielectric  system 
with  three  different  ee.  The  Coulomb  interaction  is  strongly 
enhanced  for  a  low -k  dielectric  environment  and  is  damped 
for  the  high-/c  case. 

When  a  point  charge  is  located  inside  a  3D  semicond¬ 
uctor,  its  Coulomb  potential  is  lowered  by  the  dielectric 


constant  of  the  semiconductor  host  alone.  For  thin  semi¬ 
conductor  layers,  the  Coulomb  potential  is  determined 
by  the  dielectric  constants  of  both  the  semiconductor  itself 
and  the  surrounding  dielectrics.  When  a  high  density  of 
mobile  carriers  is  present  in  the  semiconductor,  the 
Coulomb  potential  is  further  screened.  For  atomically  thin 
semiconductors,  understanding  the  dielectric  mismatch 
effect  on  the  free-carrier  screening  of  scattering  potentials 
is  necessary.  At  zero  temperature,  static  screening  by  the 
2D  electron  gas  is  captured  by  the  Lindhard  function  [21]: 

£2d(q,  «J  -►  o)  =  1  +- - n (q,co  0)(4>!  +  $2)>  (2) 

2e0esq 


where  q  is  the  2D  scattering  wave  vector,  and  II  is  the 
polarizability  function  at  zero  temperature  [22], 


a{q-°’  - o) = ss~  { 1  - e[?  -  UF]^f~Wf\ 

(3) 

where  gs ,  gv  are  the  spin  and  valley  degeneracy  factors, 
respectively,  m*  is  the  electron  mass,  kF  is  the  Fermi  wave 
vector,  and  ©[...]  is  the  Heaviside  unit-step  function.  The 
function  is  the  form  factor,  and  <F2  is  the  dielectric 
mismatch  factor,  which  are  defined  by  the  equations  [23] 

$1  =y  X2{z)dz  j X2(z')  exp(-q\z  -  z'\)dz',  (4) 

$  _  2x+x_  exp (~qa)(ee  -  e,)2  -  (jg.  +  J+)(g^  ~  g?) 
exp(qa)(ee  +  es )2  -  exp {-qa)(ee  -  es)2 

(5) 

where  x±  =  /  <fcexp(zbgz)/2(z).  The  free-carrier  screen¬ 
ing  is  taken  into  account  by  dividing  the  unscreened 
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scattering  matrix  elements  by  e2( j.  Equation  (2)  can  be 
recast  as  the  Thomas-Fermi  formula:  e2d  —  1  +  qFF/q,  in 
analogy  to  the  case  in  the  absence  of  a  dielectric  mismatch. 
Here,  q corresponds  to  the  Thomas-Fermi  screening 
wave  vector  qFF  without  a  dielectric  mismatch.  Figure  2(a) 
shows  the  ratio  qFF /  qFF  that  captures  the  effect  of  the 
dielectric  mismatch  on  screening  at  zero  temperature.  The 
2D  electron  density  is  ns  ~  1012  cm-2  in  this  figure.  As 
can  be  seen,  the  free-carrier  screening  is  weakened  by  a 
high-/c  dielectric,  and  it  is  enhanced  in  the  low-/c  case.  This 
dependence  is  opposite  to  the  effect  of  the  dielectric 
environment  on  the  net  unscreened  Coulomb  interaction. 

The  momentum  relaxation  rate  (rm)_1  due  to  elastic 
scattering  mechanisms  is  evaluated  using  Fermi’s  golden 
rule  in  the  form 


1 


T 


m 


2 n  [  d2k'  \Mkkf\2 
~hj  (2n)2  e2d 


(1-cos  0)S(Ek^Ek,), 


(6) 


where  Mkk>  is  the  matrix  element  for  scattering  from  state  k 
to  k ',  6  is  the  scattering  angle,  and  Ek  and  Ek>  are  the 
electron  energies  for  states  k  and  k ',  respectively.  For  the 
charged  impurity  scattering  momentum  relaxation  rate 
(t^)_1,  the  scattering  matrix  element  is  evaluated  as 


e2  1  f  y  4;r2sinh(^) 

2e0ssS  q  \exp(ga)  —  y  47i2(qa)  +  ( qa )3 
2[l-exp(-f)]^2  +  M2| 

4  n2  (qa)  +  (qa)3  J 


(7) 


Figure  2(b)  shows  (t^)_1  with  the  impurity  density  of 
Nj  ~  1012  cm-2.  £e  and  ns  are  varied  over  2  orders  of 
magnitude  to  map  out  the  parameter  space.  Evidently, 
( Tcm)~l  still  reduces  monotonically  with  increasing  ee 


FIG.  2.  Effect  of  dielectric  mismatch  on  the  (a)  free-carrier 
screening  and  (b)  Coulomb  momentum  relaxation  rate  at  zero 
temperature.  The  inset  of  (a)  shows  schematically  the  scattering 
angle  for  different  electron  densities. 


because  the  weakening  of  the  unscreened  Coulomb  poten¬ 
tial  is  stronger. 

The  reduction  of  ( Tcm)~l  for  a  high-/c  environment  is 
much  enhanced  for  high  ns> ,  as  indicated  in  Fig.  2(b). 
When  ee  varies  from  1  to  100,  ( Tcm)~l  decreases  about 
~1.4  times  for  ns  ~  1011  cm-2,  and  about  ~2.6  times  for 
ns  ~  1013  cm-2.  From  the  perspective  of  screening,  notice 
from  Fig.  2(a)  that  in  a  low-/c  environment,  qFF  is  higher  for 
small-angle  scattering  events.  This  means  the  smaller  the 
scattering  angle,  the  stronger  is  the  screening.  Thus  screen¬ 
ing  favors  randomizing  the  electron  momentum.  A  high-/c 
environment  reverses  this  process:  small  angle  scattering 
events  are  weakly  screened,  and  thus  such  scattering  events 
are  favored.  Thus,  as  ee  increases,  the  electron  transport 
become  more  directional.  Though  qFF  decreases,  the  net 
screening  efficiency  increases.  These  tendencies  are 
enhanced  as  ns  increases.  From  the  scattering  potential 
point  of  view,  a  higher  ns  leads  to  a  larger  Fermi  wave 
vector  kF.  As  shown  schematically  in  the  inset  of  Fig.  2(a), 
the  same  q—  |kz-  —  ky|  with  high  ns  corresponds  to  a 
smaller  scattering  angle  than  a  lower- /?s,  case,  leading  to  a 
reduced  ( Tcm)~l .  This  effect  on  the  Coulomb  scattering 
matrix  element  is  multiplied  by  the  dielectric  mismatch 
factor;  thus,  a  high -ns  system  shows  stronger  ee  depend¬ 
ence  at  zero  temperature. 

For  finite  temperatures,  following  Maldague  [22,24,25], 
the  static  polarizability  function  is 


n  (q,T,EF)  =  J 
0 


n (q,  co  ^  0) 

4kBTcosh2[(EF  -  E)/2kBT } 


dE, 


(8) 


where  EF  is  the  Fermi  energy  and  kB  is  the  Boltzmann 
constant.  Figure  3(a)  shows  the  calculated  temperature- 
dependent  polarizability  normalized  to  the  zero- 
temperature  value  at  different  ns.  The  electron  gas  is  less 
polarizable  at  higher  temperatures  and  lower  ns. 
Polarizability  is  caused  by  the  spatial  redistribution  of 
the  electron  gas  induced  by  the  Coulomb  potential;  thus,  it 
is  proportional  to  ns.  As  temperature  increases,  the  thermal 
energy  randomizes  the  electron  momenta,  accelerating  the 
transition  of  the  electron  system  back  into  an  equilibrium 
distribution,  consequently  weakening  the  polarization. 
The  decrease  of  polarizability  reduces  the  free-carrier 
screening.  Figure  3(b)  shows  the  temperature-dependent 
Coulomb-scattering-limited  mobility  (//imp)  at  two  different 
ns.  The  dielectric  mismatch  effect  is  more  significant  for 
low  ns  because  of  the  fast  decrease  of  the  polarizability 
with  increasing  temperature.  For  high  ns ,  on  the  other  hand, 
the  dielectric  mismatch  effect  is  not  as  drastic.  The  shape  of 
the  temperature-dependent  /qmp  curve  is  highly  dependent 
on  the  polarizability  and  ns.  Consequently,  if  the  electron 
transport  is  dominated  by  impurity  scattering,  one  can  infer 
ns  from  the  shape  of  the  temperature  dependence  of  the 
electron  mobility. 
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Temperature  (K) 


FIG.  3.  (a)  The  normalized  polarizability  and  (b)  impurity- 

limited  mobility  at  different  electron  densities  as  a  function  of 
temperature. 


“+”  and  ”  are  associated  with  phonon  emission  and 
absorption,  respectively.  For  optical  deformation  potentials 
(ODP)  [26], 

1  _Dlm*(Nq+\±$*n  l-(k?/k)cos0)d6 

^0— odp  4  nh2psm  J  4/ 

—n 

1  D\m*(Nq  +  \  ±  |)  r  4(1  —  (k'/k)  cos6)d9 
4-odp  Anh2ps(o  J  4/ 

—71 

(12) 

where  D  is  the  optical  deformation  potential,  Nq  = 
l/[exp(hco/kBT)  —  1]  is  the  Bose-Einstein  distribution 
for  optical  phonons  of  energy  hco ,  and  the  subscripts  0 
and  1  denote  the  zero-  and  first-order  ODP,  respectively. 

The  scattering  rate  by  polar-optical  (LO)  phonons  is 
given  by  the  Frohlich  interaction  [28], 


Much  interest  exists  in  using  atomically  thin  semi¬ 
conductors  as  possible  channel  materials  for  electronic 
devices,  in  which  such  layers  are  in  close  proximity  to 
dielectrics.  To  that  end,  we  investigate  both  the  intrinsic  and 
extrinsic  phonon  scattering  in  SL  MoS2.  Kaasbjerg  et  al. 
[26]  have  predicted  the  theoretical  intrinsic  phonon-limited 
mobility  (/q_ ph)  of  SL  MoS2  from  first  principles  using  a 
density-functional-based  approach.  They  estimated  a  room- 
temperature  upper  limit  for  the  experimentally  achievable 
mobility  of  about  410  cm2/Vs,  which  weakly  depended  on 
ns.  Their  estimate  did  not  include  the  effects  of  free-carrier 
screening  and  dielectric  mismatch.  In  light  of  the  strong 
effect  of  these  factors  on  the  Coulomb  scattering,  we 
evaluate  //;-ph  in  MoS2  in  the  Boltzmann  transport  formal¬ 
ism  with  the  modified  free-carrier  screening.  The  material 
parameters  for  SL  MoS2  were  obtained  from  Ref.  [27].  The 
momentum  relaxation  rate  due  to  quasielastic  scattering  by 
an  acoustic  phonon  is  given  by 

J_  =  ?  (1  —  cos  6)d0 

4C  2 nh3psv2s  J  e\d 

—n 


where  ps  is  the  areal  mass  density  of  SL  MoS2,  vs  is  the 
sound  velocity,  and  Eac  is  the  acoustic  deformation  poten¬ 
tial.  For  inelastic  electron-optical  phonon  interactions,  the 
momentum  relaxation  rate  in  the  RTA  is  obtained  by 
summing  the  emission  and  absorption  processes, 


1  &[Ek  —  flG>op\  1 

op  +  '  — 

^ m  ^"op  ^"op 


(10) 


where  afov  is  the  frequency  of  the  i/th  optical-phonon 
mode.  The  momentum  relaxation  rates  with  superscripts 


1  e2conf  1/1  1  \  /  1  i  1\ 

t/q  8 nh2  6*0  Uoo  e  J  \  q  +  2  2 ) 

x  (13) 

J  <1  £2d 

—Jl 

where  e ^  is  the  high-frequency  relative  dielectric  constant, 
and  Tq  is  the  form  factor  defined  by  Eq.  (4). 

Figure  4(a)  shows  the  iq-dependent  screened  /q  ph  at 
room  temperature.  For  comparison,  the  unscreened  /q_ph  is 


FIG.  4.  Electron  mobility  in  MoS2  due  to  intrinsic  phonon 
scattering  at  room  temperature  with  the  electron-phonon  interaction 
(a)  fully  screened  and  (b)  partially  screened.  The  dashed  lines  show 
mobilities  limited  by  unscreened  phonon  modes,  and  the  solid  lines 
show  the  mobilities  limited  by  fully  screened  modes. 
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also  shown  as  a  reference  (blue  line).  The  unscreened 
values  remain  effectively  constant  (about  380  cm2 /Vs)  over 
the  range  of  ns  of  interest  (10n-1013  cm-2).  This  is  in 
agreement  with  the  previous  predictions  (32CM-10  cm2/Vs) 
[26,29].  However,  the  screened  increases  sharply  with 
increasing  ns.  As  can  be  seen  in  Fig.  4(a),  introducing  a 
high-/c  dielectric  leads  to  a  reduction  of  ^/  ph;  the  highest 
values  of  ^/  ph  reduce  from  3100  to  1500  cm2/Vs  as  ee 
increases  from  about  ~7.6  to  about  ~20.  The  strong 
dependence  of  fiiph  on  the  dielectric  environment  is  entirely 
due  to  the  dielectric-mismatch  effect  on  free-carrier  screen¬ 
ing  since  the  unscreened  phonon-scattering  matrix  element 
is  not  affected  by  ee.  Over  the  entire  range  of  ns ,  longitudinal 
optical  phonon  scattering  is  dominant.  This  finding  is 
different  from  previous  works  on  multilayer  MoS2  transport 
where  the  room-temperature  /i/  ph  was  determined  by 
homopolar  phonon  scattering  [30-32]. 

We  have  used  the  static  dielectric  function  for  calculat¬ 
ing  the  screened  interactions  due  to  different  modes  of 
phonons  in  the  limit  go  0.  Scattering  mechanisms  via 
long-range  Coulomb  interactions,  such  as  charged  impu¬ 
rities,  polar-optical  phonons,  and  piezoelectric  acoustic 
phonons,  can  be  effectively  screened  by  free  carriers. 
However,  free  carriers  may  not  respond  to  rapidly  chang¬ 
ing  scattering  potentials  originating  from  short-range 
interactions.  There  are  arguments  about  to  what  extent 
the  short-range  deformation  potentials  induced  by  acoustic 
(ADP)  and  optical  phonons  (ODP)  are  screened  by  free 
carriers.  Boguslawski  and  Mycielski  [33]  argue  that  in  a 
single- valley  conduction  band,  the  deformation  potentials 
(both  ADP  and  ODP)  are  screened  in  the  same  way  as  the 
macroscopic  (long-range)  phonon  potentials.  But  for 
multivalley  semiconductors  (Ge),  only  the  longitudinal 
acoustic  (LA)  mode  of  the  ADP  can  be  effectively 
screened  by  free  carriers.  The  free-carrier  screening  of 
the  transverse  acoustic  (TA)  mode  ADP  and  ODP  can,  to  a 
good  approximation,  be  neglected.  [34].  In  SL  MoS2, 
Kaasbjerg  et  al.  [27]  have  argued  that  the  LA  mode  of  the 
ADP  can  be  treated  as  screened  by  the  long-wavelength 
dielectric  function,  while  the  screening  of  the  TA  mode 
ADP  by  free  carriers  can  be  neglected. 

Figure  4(b)  highlights  the  effect  of  the  partially  screened 
electron-phonon  interaction  compared  to  the  fully  screened 
version  in  Fig.  4(a).  For  the  plot  in  Fig.  4(b),  we  have 
screened  the  polar-optical  and  LA  phonon  scattering  as 
in  Fig.  4(a),  and  we  leave  the  TA  and  ODP  interactions 
unscreened.  The  highest  ///  ph  reached  by  free-carrier 
screening  effects  is  reduced  to  about  750  cm2/Vs  by  not 
screening  the  DP  modes.  The  mobility  is  dominated  by 
the  polar-optical  phonon  interaction  at  low  carrier  density 
and  by  TA  and  ODP  at  moderate  and  high  densities.  The 
scattering  of  electrons  due  to  piezoelectric  phonons  is  not 
considered  because  it  is  relevant  only  at  very  low  temper¬ 
atures  and  because  there  are  still  uncertainties  in  the 
piezoelectric  coefficients  of  SL  MoS2  [27,35]. 


In  both  cases,  the  calculated  room-temperature  ///  ph  are 
much  higher  than  reported  experimental  values,  implying 
that  there  is  still  much  room  for  improvement  of  mobilities 
in  atomically  thin  semiconductors.  For  the  rest  of  this  work, 
we  use  the  fully  screened  intrinsic  phonon  scattering,  as 
shown  in  Fig.  4(a).  To  pinpoint  the  most  severe  scattering 
mechanisms  limiting  the  mobility  in  current  samples,  we 
discuss  an  extrinsic  phonon-scattering  mechanism  at  play 
in  these  materials,  again  motivated  by  similar  processes  in 
graphene. 

Electrons  in  semiconductor  nanoscale  membranes  can 
excite  phonons  in  the  surrounding  dielectrics  via  long- 
range  Coulomb  interactions,  if  the  dielectrics  support  polar 
vibrational  modes.  Such  “remote  phonon”  or  “surface- 
optical”  (SO)  phonon  scattering  has  been  investigated 
recently  for  graphene  and  found  to  be  far  from  negligible 
[15-17].  SO  phonon  scattering  can  severely  degrade 
electron  mobility;  however,  this  process  has  not  been 
studied  systematically  in  atomically  thin  semiconductors. 
The  electron-SO  phonon  interaction  Hamiltonian  is 
[15,17,18] 


H 


e-SO 


-qz 


—  (eWap  +  e-Wa*) 
.  v  Q 


(14) 


where  aq+(aq)  represents  the  creation  (annihilation)  oper¬ 
ator  for  the  zdh  SO  phonon  mode.  Neglecting  the  dielectric 
response  of  the  atomically  thin  MoS2  layer  in  lieu  of  the 
surrounding  media,  the  electron-SO  phonon  coupling 
parameter  Fv  is 


^so  f _ 

256*0  \£™x  + 


(15) 


where  £™x  (6*qX)  is  the  high-  (low-)  frequency  dielectric 
constant  of  the  dielectric  hosting  the  SO  phonon,  and  6*^,  is 
the  high-frequency  dielectric  constant  from  the  dielectric 
on  the  other  side  of  the  membrane.  The  frequency  of  the  SO 
phonon  govso  is  [17,36] 


^so  =  ** 


TO 


+  6*° 


+  e° 


1/2 


(16) 


where  G0jO  is  the  vth  bulk  transverse  optical-phonon 
frequency  in  the  dielectric.  The  scattering  rate  due  to  the 
SO  phonon  is  then  given  by 


1  32n3e2F2vm*S  /  1  ,  1\ 

ib~  (^+2±2j 

j  1  sinh2(f)  (l -(k'/k) cos O)d0 

XJ  q{Ajq  +  aW?  4T 
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Table  I  summarizes  the  parameters  for  some  commonly 
used  dielectrics. 

Figure  5  shows  the  room-temperature  electron  mobility 
for  various  dielectric  environments  for  two  representative 
temperatures,  100  K  and  300  K.  Nj  and  ns  are  both  about 
1013  cm-2.  The  solid  lines  show  the  net  mobility  by 
combining  the  scattering  from  charged  impurities,  and 
intrinsic  and  SO  phonons,  whereas  the  dashed  lines  show 
the  cases  neglecting  the  SO  phonons.  When  SO  phonon 
scattering  is  absent,  the  electron  mobility  is  limited  almost 
entirely  by  //imp,  which  increases  with  ee  because  of  the 
reduction  of  Coulomb  scattering  by  dielectric  screening. 
The  addition  of  the  SO  phonon  scattering  does  not  change 
things  much  at  100  K,  except  for  the  highest  ee  case 
(Hf02/Zr02).  But  it  drastically  reduces  the  electron  mobil¬ 
ity  at  room  temperature,  as  is  evident  in  Fig.  5.  For  instance, 
neglecting  SO  phonon  scattering,  one  may  expect  that 
by  using  Hf02/Zr02  as  the  dielectrics  instead  of  Si02/air, 
the  RT  mobility  //imp  should  improve  from  about  ~45  to 
80  cm2/Vs.  However,  when  the  SO  phonon  scattering  is 
in  action,  the  mobility  in  the  Hf02/MoS2/Zr02  structure  is 
actually  degraded  to  around  25  cm2/Vs,  even  lower  than 
the  Si02/air  case.  Thus,  SL  MoS2  layers  suffer  from 
enhanced  SO  phonon  scattering  if  they  are  in  close 
proximity  to  high-/c  dielectrics  that  allow  low-energy  polar 
vibrational  modes. 

To  calibrate  our  calculations,  we  study  the  temperature- 
dependent  electron  mobility  for  SL  MoS2  embedded 
between  Si02  and  Hf02  and  compare  the  calculations 
with  reported  experimental  results.  This  structure  is  often 
used  in  top-gated  MoS2  field  effect  transistors  (FETs);  thus, 
understanding  the  transport  in  it  provides  a  pathway  to 
understanding  the  device  characteristics.  In  Fig.  6(a),  the 
blue  curves  indicate  calculated  values  of  //imp  with  different 
A/,  and  the  red  line  shows  the  SO  phonon-scattering 
limited  mobility  (//so)’  with  ns  ~  1013  cm-2.  The  temper¬ 
ature-dependent  juso  of  each  SO  phonon  mode  follows 
the  Arrhenius  rule:  juso  cx  exp (hco0/kBT),  and  the  net  //so 
is  dominated  by  the  softest  phonon  mode  with  the  lowest 
energy.  The  black  curves  indicate  the  net  mobilities 
considering  all  scattering  mechanisms  discussed  in  this 
work.  The  open  squares  are  the  experimental  results 
measured  by  the  Hall  effect  on  SL  MoS2  FETs  from 
Ref.  [4].  The  Nj  and  ns  necessary  to  fit  the  data  are 
indicated  in  Fig.  6(a).  At  low  temperatures,  the 


TABLE  I.  SO  phonon  modes  for  different  dielectrics. 


Si02a 

AlNa 

BNb 

Al203a 

Hf02 

Zr02‘ 
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cox 
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24 
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3.2 

5.03 

4 

^SO 

55.6 

81.4 

93.07 

48.18 

12.4 

16.67 

^SO 

138.1 

88.5 

179.1 

71.41 

48.35 

57.7 

aRef.  [15] 
bRef.  [37] 


FIG.  5.  Electron  mobility  as  a  function  of  an  environment 
dielectric  constant.  Dashed  lines  show  the  mobility  without 
considering  the  SO  phonons. 

experimental  electron  mobility  in  SL  MoS2  is  entirely 
limited  by  //imp.  This  is  really  not  unexpected;  it  took 
several  decades  of  careful  epitaxial  growth  and  ultraclean 
control  to  achieve  the  high  mobilities  in  III-V  semicon¬ 
ductors  at  low  temperatures.  Based  on  this  study,  we  predict 
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FIG.  6.  (a)  Temperature-dependent  electron  mobility  (black 

lines)  in  Si02/MoS2/Hf02  structure.  The  blue  lines  indicate  //imp 
and  the  red  lines  show  //so.  Open  squares  show  experimental 
results  from  single-layer  MoS2  FETs  from  Ref.  [4].  (b)  Room- 
temperature  phonon-determined  electron  mobilities//ph  and 
(c)  the  critical  impurity  densities  Ncr  corresponding  to  //imp  = 
//ph  in  SL  MoS2  surrounded  by  different  dielectrics.  Dashed  lines 
show  the  fitted  //ph  and  Ncr. 
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that  the  low-temperature  mobilities  in  atomically  thin 
semiconductors  can  be  significantly  improved  by  lowering 
the  impurity  density.  The  room-temperature  mobility  in 
III-V  semiconductors  is  limited  by  intrinsic  polar-optical 
phonon  scattering.  For  comparison,  we  find  that  for  SL 
MoS2,  the  room-temperature  mobility  is  considerably 
degraded  by  SO  phonon  scattering,  even  with  Nj  as  high 
as  6  x  1012  cm-2,  as  shown  in  Fig.  6.  When  SO  phonon 
scattering  is  absent,  the  room-temperature  mobility  is 
expected  to  be  about  130  cm2/Vs  with  Nj  =  6x  1012cm-2, 
but  the  measured  values  are  typically  lower  (about 
50  cm2/Vs).  Consequently,  using  Hf02  as  gate  dielectrics 
can  modestly  improve  //imp.  However,  the  strong  SO 
phonon  scattering  that  comes  with  Hf02  can  severely 
decrease  the  high-temperature  electron  mobility  in  clean 
MoS2  with  low  charged  impurity  densities. 

An  important  question  then  is,  which  dielectric  can  help 
improve  the  room-temperature  electron  mobility  in  SL 
MoS2?  To  answer  that  question,  in  Fig.  6(b),  we  plot  the 
room-temperature  (intrinsic  +  SO)  phonon-limited  electron 
mobility  (//ph)  in  SL  MoS2  surrounded  by  different  dielec¬ 
trics.  From  the  overall  trend,  //ph  decreases  with  increasing 
se,  and  suspended  SL  MoS2  shows  the  highest  potential 
electron  mobility  (over  10, 000  cm2/Vs).  It  is  worth  noting 
that  if  the  scattering  of  electrons  by  intrinsic  phonons  is 
only  partially  screened,  as  shown  in  Fig.  4(b),  the  highest 
achievable  mobility  in  SL  MoS2  will  be  an  order  lower 
(around  1000  cm2/Vs).  However,  these  high  values  are 
attainable  in  suspended  SL  MoS2.  Because  //ph  for  MoS2 
surrounded  by  high-/c  materials  is  dominated  by  SO  phonon 
scattering,  the  values  do  not  vary  much.  The  critical 
impurity  densities  (Acr)  corresponding  to  //imp  =  //ph  are 
shown  in  Fig.  6(c).  As  long  as  Nj  >  Ncr ,  //imp  completely 
masks  juph.  When  Nj  <  Ncr,  the  electron  mobility  becomes 
dominated  by  phonons  and  moves  towards  the  upper  limit. 
High  fiph  indicates  a  greater  potential  for  attaining  higher 
electron  mobilities.  However,  we  also  need  the  sample  to  be 
highly  pure.  In  high-/c  environments  that  support  low- 
energy  polar  vibrational  modes,  there  is  not  as  much  room 
for  improving  the  electron  mobility  as  in  low-/c  structures. 
A  compromise  is  seen  for  Aluminum  Nitride  (AIN)-  and 
Boron  Nitride  (BN)-based  dielectrics,  which  by  virtue  of 
the  light  atom  N,  allows  high-energy  optical  modes  in 
spite  of  their  polar  nature.  From  Figs.  6(b)  and  6(c),  one 
can  obtain  two  useful  relationships  for  single-layer  MoS2: 
//ph  ~  35000/£2-2  cm2/Vs  and  Ncr  ~  1010£2-5  cm-2,  with 
ns  set  at  a  typical  on-state  carrier  density  of  1013  cm-2,  as 
shown  by  dashed  lines.  These  empirical  relations  should 
guide  the  proper  choice  of  dielectrics  and  the  maximum 
allowed  impurity  densities. 

To  further  illustrate  the  relative  importance  of  SO 
phonon  and  charged  impurity  scattering  in  SL  MoS2,  we 
vary  Nj  and  ns  in  different  dielectric  environments  and 
check  the  changing  trends  of  electron  mobilities  at  room 
temperature.  Figure  7(a)  shows  the  net  electron  mobilities 


in  SL  MoS2  as  a  function  of  Nj  with  ns  =  1013  cm-2. 
Figures  7(b)  and  7(c)  show  the  electron  mobility  as  a 
function  of  ns  for  Nj  =  1011  and  1013  cm-2.  The  electron 
mobility  is  weakly  dependent  on  the  dielectric  environment 
at  high  Nj  (>  1013  cm-2),  as  shown  in  the  dashed  box  in  the 
bottom  right  corner  of  Fig.  7(a).  Within  this  window,  high-/c 
dielectrics  can  improve  the  mobility,  but  only  very  nomi¬ 
nally  because  the  unscreened  mobilities  are  already  quite 
low.  When  Nj  is  lowered  below  about  1012  cm-2,  a  low-/c 
environment  shows  higher  electron  mobility.  For  most  of 
the  dielectric  environments,  when  Nj  >  1012  cm-2,  the 
mobility  fits  the  following  empirical  impurity-scattering- 
dominated  relationship:  // «4200/[A7/10n  cm-2]cm2/Vs, 
as  shown  by  the  dashed  line  in  Fig.  7(a).  Using  this 
expression,  one  can  estimate  Nj  from  measured  electron 
mobility  for  high  ns.  As  ns  decreases,  electron  mobilities 
in  different  dielectric  environments  start  to  separate 
from  each  other,  as  shown  in  Fig.  7(c).  In  this  case,  the 
electron  mobilities  can  fit  the  following  relationship: 

[VeJ  +  (lofe)L2]cm2/Vs  for  ns< 
1013  cm-2,  shown  as  dashed  lines  in  Fig.  7(c).  A(ee)  is 
a  fitting  constant  depending  on  se,  and  some  values  are 
listed  in  the  inset  table  of  Fig.  7(c).  High-/c  dielectrics  with 
low-energy  phonons  (Hf02,  Zr02)  severely  degrade  the 
electron  mobility  over  the  entire  Nj  range  because  of 
the  dominant  effect  of  SO  phonon  scattering.  Note  that  the 
dielectric  mismatch  effect  can  be  slightly  overestimated 
here  since  we  have  assumed  the  thickness  of  the  dielectric 
to  be  infinite  [25].  In  top-gated  FETs,  the  top  dielectric 
could  be  very  thin.  Thus,  the  capability  of  improving 
electron  mobility  by  high-/c  dielectrics  can  be  even  less 
significant.  Since  most  applications  require  high  mobilities, 


FIG.  7.  The  room-temperature  net  electron  mobilities  in  SL 
MoS2,  considering  all  kinds  of  scattering  mechanisms  as  a 
function  of  (a)  Nj  with  fixed  ns  at  1013  cm-2;  (b)  and  (c)  ns 
with  Nj  fixing  at  1011  and  1013  cm-2,  respectively.  The  numbers 
on  the  curves  show  the  average  dielectric  constant  of  the 
surrounding  dielectrics.  Dashed  lines  show  the  fitted  electron 
mobilities. 
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high  ns ,  and  high  ee  to  be  present  simultaneously  in  the 
same  structure  for  achieving  the  highest  conductivities, 
AIN/A1203  or  BN/BN  encapsulation  emerges  as  the  best 
compromise  among  the  dielectric  choices  considered  here. 
One  can  also  conceive  of  dielectric  heterostructures,  with  a 
few  BN  layers  closest  to  MoS2  to  damp  out  the  SO  phonon 
scattering,  followed  by  higher-/c  dielectrics  to  enhance 
the  gate  capacitance  for  achieving  high  carrier  densities. 
All  this,  however,  requires  ultraclean  MoS2  to  start  with, 
with  Nj  well  below  1012  cm-2  to  attain  the  high  room- 
temperature  mobilities,  about  1000  cm1 2 3/Vs.  The  presence 
of  high  impurity  densities  will  always  mask  the  intrinsic 
potential  of  the  materials,  and  this  is  the  most  important 
challenge  moving  forward. 

In  conclusion,  carrier  transport  properties  in  atomically 
thin  semiconductors  are  found  to  be  highly  dependent  on 
the  dielectric  environment  and  on  the  impurity  density.  For 
current  2D  crystal  materials,  electron  mobilities  are  mostly 
dominated  by  charged  impurity  scattering.  Remote  pho¬ 
nons  play  a  secondary  role  at  high  temperature  depending 
on  the  surrounding  dielectrics.  The  major  point  is  that 
the  mobilities  achieved  to  date  are  far  below  the  intrinsic 
potential  in  these  materials.  High-/c  gate  dielectrics  can 
increase  the  electron  mobility  only  for  samples  infected 
with  very  high  impurity  densities.  Clean  samples  with 
low-/c  dielectrics  show  much  higher  electron  mobilities. 
AIN-  and  BN-based  dielectrics  offer  the  best  compromise 
if  a  high-mobility  and  high-gate  capacitance  are  simulta¬ 
neously  desired,  as  is  the  case  in  field-effect  transistors.  The 
truly  intrinsic  mobility  limited  by  the  atomically  thin 
semiconductor  itself  can  only  be  achieved  in  ultraclean 
suspended  samples,  as  is  the  case  for  graphene. 
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ABSTRACT  Layered  2-D  crystals  embrace  unique  features  of  atomically  thin  bodies,  dangling  bond  free 
interfaces,  and  step-like  2-D  density  of  states.  To  exploit  these  features  for  the  design  of  a  steep  slope 
transistor,  we  propose  a  Two-dimensional  heterojunction  interlayer  tunneling  field  effect  transistor  (Thin- 
TFET),  where  a  steep  subthreshold  swing  (SS)  of  ~14  mV/dec  and  a  high  on-current  of  ~300  /xA//rm  are 
estimated  theoretically.  The  SS  is  ultimately  limited  by  the  density  of  states  broadening  at  the  band  edges 
and  the  on-current  density  is  estimated  based  on  the  interlayer  charge  transfer  time  measured  in  recent 
experimental  studies.  To  minimize  supply  voltage  V do  while  simultaneously  maximizing  on  currents, 
Thin-TFETs  are  best  realized  in  heterostructures  with  near  broken  gap  energy  band  alignment.  Using 
the  WSe2/SnSe2  stacked-monolayer  heterostructure,  a  model  material  system  with  desired  properties  for 
Thin-TFETs,  the  performance  of  both  /2-type  and  p-type  Thin-TFETs  is  theoretically  evaluated.  Nonideal 
effects  such  as  a  nonuniform  van  der  Waals  gap  thickness  between  the  two  2-D  semiconductors  and  finite 
total  access  resistance  are  also  studied.  Finally,  we  present  a  benchmark  study  for  digital  applications, 
showing  the  Thin-TFETs  may  outperform  CMOS  and  III-V  TFETs  in  term  of  both  switching  speed  and 
energy  consumption  at  low-supply  voltages. 

INDEX  TERMS  Tunnel  FET,  2-D  crystals,  transport  model,  steep  slope,  subthreshold  swing  (SS),  layered 
materials,  benchmarking. 


I.  INTRODUCTION 

Tunnel  Field  Effect  Transistors  (FETs)  are  perceived  as 
promising  electronic  switches  that  may  enable  scaling  the 
supply  voltage  Vdd  down  to  0.5  V  or  lower  by  reduc¬ 
ing  the  subthreshold  swing  (SS)  below  60  mV/dec  at  room 
temperature. 

To  date,  numerous  Tunnel  FETs  have  been  demonstrated, 
among  which  heterostructures  with  near  broken  gap  band 
alignment  are  favored  in  order  to  achieve  sub-60  mV/dec 
SS  and  high  on  currents  simultaneously  [1].  Tunnel  FETs 
also  require  a  very  strong  gate  control  over  the  channel 
region  to  obtain  sub-60  mV/dec  SS  values;  this  in  turn 
demands  ultra-thin  body  or  nanowire  structures,  where  size 
induced  quantization  enlarges  the  bandgap  and  impedes  the 


realization  of  near  broken  gap  alignment  [2]-[4].  Layered 
2D  crystals,  such  as  monolayers  of  transition  metal  dichalco- 
genides  (TMD)  MX2  (e.g.,  M  =  Mo,  W;  X  =  S,  Se,  Te)  and 
other  metal  chalcogenides  MX*  (e.g.,  M  =  Ga,  Sn;  X  =  O, 
S,  Se)  offer  a  native  thickness  of  about  0.6  nm  with  a  vari¬ 
ety  of  bandgaps  and  band- alignments  [4],  [5].  Furthermore, 
2D  crystals  possess  a  sharp  turn  on  of  density  of  states  at 
the  band  edges  and  have  no  surface  dangling  bonds  thus 
potentially  enabling  a  low  interfacial  density  of  state,  which 
are  highly  desired  for  achieving  a  sharp  SS  [6].  Recent 
experimental  results  show  that  the  band  alignment  in  stacked- 
monolayer  2D  crystal  heterostructures  can  be  tuned  by  an 
external  electric  field  perpendicular  to  the  heterojunction 
plane  [7]  and  the  charge  transfer  in  stacked-monolayer  2D 
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crystal  heterojunctions  is  reasonably  fast  [8].  In  such  a 
context,  we  propose  the  Two-dimensional  Heterojunction 
Interlayer  Tunneling  FET  (Thin-TFET)  based  on  a  verti¬ 
cal  arrangement  of  2D  layered  materials.  In  particular,  we 
discuss  both  ft-type  and  p- type  Thin-TFETs  employing  a 
promising  material  system  of  2H-WSe2  and  lT-SnSe2.  Our 
simulations  suggest  that  very  competitive  SS  values  and  a 
high  on-current  can  be  achieved  in  the  Thin-TFETs.  Along 
with  the  low  intrinsic  gate-to-drain  and  gate- to- source  capac¬ 
itances  in  comparison  to  CMOS  and  p-i-n  III-V  TFETs 
benchmarked  in  Section  III-D,  the  Thin-TFETs  enable  fast 
switching  and  low  energy  consumption.  The  effect  of  a 
non-uniform  van  der  Waals  gap  thickness  and  the  external 
source  and  drain  total  access  resistance  are  also  discussed. 
At  the  end  of  the  paper,  we  will  also  share  some  insights 
on  the  experimental  realization  of  Thin-TFETs  derived  from 
the  ongoing  investigations  in  our  laboratory. 
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FIGURE  1.  Schematic  device  cross  section  of  a  Thin-TFET. 


II.  DEVICE  STRUCTURE  AND  MODELING  APPROACH 

The  Thin-TFET  device  structure  is  shown  in  Fig.  1,  where  the 
bottom  and  top  2D  semiconductors  act  as  the  source  and  the 
drain  respectively.  A  van  der  Waals  gap  separates  the  top 
and  bottom  2D  semiconductors  and  the  thickness  of  the  van 
der  Waals  gap  is  defined  as  the  distance  from  the  center 
of  the  chalcogenide  atom  in  the  top  2D  layer  to  the  center 
of  the  nearest  chalcogenide  atom  in  the  bottom  2D  layer 
(see  Fig.  1).  The  device  working  principle  can  be  explained 
as  follows:  take  the  /7-type  Thin-TFET  as  the  example,  when 
the  conduction  band  edge  of  the  bottom  2D  semiconductor 
E cb  is  higher  than  the  valence  band  edge  of  the  top  2D 
semiconductor  Eyy  (see  Fig.  2),  tunneling  from  the  bottom 
layer  is  inhibited  and  the  device  is  nominally  off.  When  a 
negative  top  gate  voltage  pulls  E yy  above  E cb  (see  Fig.  3(a)), 
a  tunneling  window  is  opened  thus  current  can  flow. 

To  calculate  the  band  alignment  between  Ecb  and  Eyj 
along  the  direction  perpendicular  to  the  2D  semiconductors 
we  first  use  Gauss’s  law  and  write  [9] 

CtoxVtox  —  CvdwVvdw  =  e(pr  —  nr  +  Ay) 
CboxVbox  +  CvdwVvdw  =  e(ps  -  m  +  Nb)  (1) 

where  e  is  the  magnitude  of  an  electron  charge,  Ct(b)OX 
is  the  capacitance  per  unit  area  of  top  (back)  oxide,  and 
Cvdw  is  the  capacitance  per  unit  area  of  the  van  der  Waals 
gap.  Vt(b)OX  and  Vvdw  are  the  corresponding  potential  drops. 
n( p)t{B )  is  the  electron  (hole)  density  in  the  top  (bottom) 


2D  semiconductor  layer,  and  Ay,  Nb  are  the  net  chemical 
doping  concentrations  (donor  minus  acceptor)  in  the  layers, 
which  are  set  to  zero  in  this  work.  The  potential  drops  can 
be  written  in  terms  of  the  top  gate  VyG,  back  gate  \bg>  and 
drain- source  voltage  Yds  (which  sets  the  split  of  the  quasi- 
Fermi  levels  in  the  top  and  bottom  semiconductor  layers), 
and  of  the  material  properties  as 


eVvdw  =  eVos  ~  e(/)p,B  —  ecf)nj  +  Eqb  +  X2d,b  —  X2d,t 
eVrox  =  eVrG  +  e<t>n,T  —  eVos  +  X2d,t  —  e$>M,T 
eVsox  =  eVsG  —  +  Eqb  +  X2d,b  +  e®M,B  (2) 


where  we  define  ecj)nj(B)=EcT{B)—EFT(B )  and 
e(/)p,T(B)=EFT(B)-Evr(B ),  EGB  is  the  energy  gap  in 
the  bottom  2D  semiconductor  and  Eft(B)  is  the  Fermi  level 
in  the  top  and  bottom  layers,  X2Dj{B)  is  the  electron  affinity 
of  the  top  (bottom)  2D  semiconductor,  and  <&m,t(B)  is  the 
metal  workfunction  of  the  top  (back)  gate  (see  Fig.  2). 

Using  the  effective  mass  approximation  and  assuming  that 
the  majority  carriers  of  the  two  2D  semiconductors  are  at 
thermodynamic  equilibrium  with  their  Fermi  levels  [10],  the 
carrier  densities  can  be  written  as 


gvm*  (m*)  kBT 
n(p)  =  - - m 


7T 


h2 


,  <#/yr(0y,5)  .  1 

exp| — 1+1 


(3) 


where  gv  is  the  valley  degeneracy  and  m*(m*)  is  the  con¬ 
duction  (valence)  band  effective  mass,  and  the  rest  of  the 
parameters  assume  their  common  meanings. 

By  inserting  Eqs.  2  and  3  in  Eq.  1  we  obtain  two  equations 
determining  0n?y,  0^#  and  thus  the  band  alignment. 

We  calculate  the  tunneling  current  by  using  the  transfer- 
Hamiltonian  method  [11],  which  was  also  recently  revisited 
for  resonant  tunneling  graphene  transistors  [12],  [13].  We 
here  summarize  the  basic  equations;  a  more  thorough  dis¬ 
cussion  can  be  found  in  our  earlier  work  [9].  The  tunneling 
current  density,  Jy,  is  expressed  as  [9]: 


_  gyg  \Mbo\2A  2KTvdw 

4: Ph 

x  J  J dkj  dkB  Sp(q )  Se(Eb  —  Ej)  (fB  —fr)  (4) 

k  t  k  b 


where  k  is  the  decay  constant  of  the  wave-function  in  the 
van  der  Waals  gap  [12],  [13],  Tvdw  is  the  thickness  of  the 
van  der  Waals  gap,  k t(b),  Et{b)  and  fj{B)  are  the  wave- 
vector,  the  energy  and  Fermi  occupation  function  in  the  top 
(bottom)  2D  semiconductor  and  Mbo  is  the  tunneling  matrix 
element  [9],  which  is  a  property  of  the  material  system  and 
is  further  discussed  in  Section  III.  Equation  4  assumes  that 
in  the  tunneling  process  electrons  interact  with  a  random 
scattering  potential,  whose  spectrum  is  taken  as  S y(g)  = 
ttL2c/(\  +q2L^/2)3/2,  where  g=|ky— ky|  and  LG  is  the  cor¬ 
relation  length.  The  scattering  relaxes  the  momentum  conser¬ 
vation,  i.e.,  allowing  tunneling  for  ky^ky.  A  similar  S y(g) 
has  been  used  to  analyze  the  resonance  linewidth  in  graphene 
tunneling  transistors  [13].  The  S'y(q)  may  be  representa¬ 
tive  of  different  scattering  mechanisms  that  are  discussed 
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in  [9]  and  [13].  The  energy  broadening  in  the  2D  semi¬ 
conductors  is  described  by  5'^(JE)=exp(— E1  /o1) /(^/ncr1), 
where  a  is  the  energy  broadening  parameter  [9]. 

Finally,  after  discussing  the  intrinsic  device  performance, 
the  contact  resistance  is  included  in  our  model  by  self- 
consistently  calculating  the  tunnel  current  density  and  the 
voltage  drop  on  the  total  access  resistance.  The  effect  of  the 
lateral  resistance  in  the  intrinsic  Thin-TFET  has  been  dis¬ 
cussed  in  our  prior  work  [14].  The  key  finding  is  that:  when 
the  tunnel  current  is  sufficiently  low  (~  1  iiAJ/im  in  the  sub¬ 
threshold  region),  the  tunnel  junction  resistance  associated 
with  the  vertical  current  flow  is  much  higher  than  the  lateral 
resistance  of  the  2D  semiconductor  source  and  drain  layers; 
as  a  result,  the  current  distribution  across  the  junction  is 
rather  uniform  laterally  in  the  sub-threshold  region. 

III.  SIMULATION  RESULTS  AND  DISCUSSIONS 
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FIGURE  2.  An  example  to  realize  both  n-type  and  p-type  Thin-TFETs  using 
one  pair  of  2-D  semiconductors  (2H-WSe2  and  lT-SnSe2)  with  near  broken 
gap  band  alignment.  For  the  n-type  Thin-TFET,  SnSe2  is  the  top  (i.e.,  drain) 
2-D  layer  and  WSe2  is  the  bottom  (i.e.,  source)  2-D  layer,  along  with  the 
top  and  back  gate  labeled  as  n-type  in  blue.  While  for  the  p- type  Thin-TFET, 
WSe2  is  the  top  (i.e.,  drain)  2-D  layer  and  SnSe2  is  the  bottom  (i.e.,  source) 
2-D  layer,  along  with  the  top  and  back  gate  labeled  as  p- type  in  red;  band 
gaps,  electron  affinities,  effective  masses  are  shown  for  WSe2  and  SnSe2. 
The  n-type  and  p- type  metal  work  functions  are  tuned  to  give  symmetric 
threshold  voltages  for  the  n-type  and  p- type  Thin-TFETs. 


A.  MATERIAL  SYSTEM  AND  N-TYPE  &  P-TYPE  THIN-TFETS 


(e) 


-0.2  0.0  0.2 
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FIGURE  3.  For  the  n-type  and  p- type  Thin-TFETs  shown  in  Fig.  2.  (a)  Band 
alignment  versus  VTC.  (b)  Current  density  versus  VTC,  the  average  SS  is 
calculated  from  10-3  fiA/fim  to  10  (c)  Current  density  versus  VDS 

at  various  VTC.  (d)  Transconductance  versus  VTC.  (e)  Carrier  concentration 
in  the  top  and  bottom  2-D  layers  versus  VTC  at  various  VDS.  (f)  Quantum 
capacitances  of  the  top  and  bottom  2-D  layers  versus  VTC  at  various  VDS. 


Out  of  various  2D  semiconductors  studied  by  density  func¬ 
tion  theory  calculations  [5]  and  experimental  efforts,  we 
chose  the  trigonal  prismatic  coordination  monolayer  (2H) 
WSe2  and  the  octahedral  coordination  (Cdl2  crystal  struc¬ 
ture)  monolayer  (IT)  SnSe2  (see  Fig.  2).  WSe2/SnSe2 
stacked-monolayer  heterojunction  can  potentially  form  a  near 
broken  band  alignment,  which  reduces  the  voltage  drop  in  the 
van  der  Waals  gap  in  the  on-state  condition  [1].  Since  there 
is  no  experimental  band  alignment  reported  for  monolayer 
WSe2  and  SnSe2,  the  band  alignment  of  the  WSe2/SnSe2 
system  used  in  this  work  are  based  on  the  existing  exper¬ 
imental  results  of  multilayer  WSe2  and  SnSe2  [15]— [17], 
while  their  approximated  effective  masses  are  based  on  the 
DFT  results  of  monolayer  WSe2  and  SnSe2  [5]  (see  Fig.  2). 


Following  the  complex  band  method  [18],  we  assume 
the  effective  barrier  height  E#  of  the  van  der  Waals  gap 
is  1  eV  and  the  electron  mass  in  the  van  der  Waals  gap 
is  the  free  electron  mass  mo,  thus  the  decay  constant  is 
k  =  +j2moEB/h  =  5.12  nm-1.  In  our  model,  we  set  the 
scattering  correlation  length  Lc  in  Sp  (q)  to  Lc= 10  nm,  which 
is  also  consistent  with  the  value  employed  in  [13];  the  energy 
broadening  o  is  set  to  be  10  meV.  Mbo  in  Eq.  4  is  directly 
related  to  the  interlayer  charge  transfer  time  r  across  the 
van  der  Waals  gap,  which  can  be  written  as  [19] 

T-1  =  p  \Mbo  1 2e~2lcTvdw  Sp  (q)  (5) 

h 
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where  p=gvm*/nh2  is  the  density  of  states  (DOS).  As  can  be 
seen  from  Eq.  5  and  the  expression  of  the  scattering  potential 
spectrum  SHq)  (given  after  Eq.  4),  due  to  scattering  in  our 
model,  t  increases  with  increasing  q,  which  is  the  magnitude 
of  the  wave-vector  difference  across  the  van  der  Waals  gap 
defined  as  q=  |k^— k#  | .  In  a  recent  experiment,  a  charge  trans¬ 
fer  time  of  25  fs  has  been  observed  across  the  van  der  Waals 
gap  between  a  stacked-monolayer  M0S2/WS2  heterostructure, 
which,  according  to  Eq.  5,  gives  us  Mpo  ~0.02  eV  when  q= 0. 
We  recognize  that  the  charge  transfer  time  might  be  different 
for  different  2D  heterojunctions,  nevertheless,  this  experimen¬ 
tally  determined  charge  transfer  time  is  a  reasonable  value 
to  use  for  the  first  pass  estimate.  Thus,  we  choose  M#o=0.02 
eV  in  all  following  simulations. 

Throughout  this  work,  the  gate  length  is  set  to  be  15 
nm,  the  back  gate  and  source  are  grounded.  An  effective 
oxide  thickness  (EOT)  of  1  nm  is  used  for  both  the  top 
and  back  oxide,  which  gives  a  top  (back)  oxide  capacitance 
Ctg  (Cbg)  of  0.518  fF/fim.  The  thickness  of  the  van  der 
Waals  gap  is  set  to  3.5  A,  unless  specified  otherwise.  We 
assume  the  relative  dielectric  constant  of  the  van  der  Waals 
gap  is  1.0,  therefore  the  van  der  Waals  gap  capacitance 
C vdw  is  0.38  fF//zm.  The  external  total  access  resistances  are 
considered  after  the  intrinsic  device  performance  is  discussed 
first  (Figs.  3  and  4). 

The  example  material  systems  for  ft-type  and  p- type  Thin- 
TFETs  based  on  the  stacked-monolayer  WSe2  and  SnSe2  are 
shown  in  Fig.  2.  The  metal  work  functions  are  tuned  to  obtain 
a  symmetric  threshold  voltage  for  the  n- type  and  the  /7-type 
Thin-TFET.  Fig.  3(a)  shows  the  band  alignment  versus  Vtg- 
Vtg  can  effectively  control  the  vertical  band  alignment  in 
the  device  by  controlling  primarily  the  band  edge  of  the  top 
(i.e.,  drain)  layer  while  having  a  weak  effect  on  the  band 
edge  of  the  bottom  (i.e.,  source)  layer,  so  that  a  tunneling 
window  is  modulated.  Fig.  3(b)  shows  Id  versus  Vtg  trans¬ 
fer  curves  with  very  compelling  average  SS  of  ~14  mV/dec 
averaged  from  10-3  /xA//zm  to  10  fiA/pm.  The  I d  versus 
Yds  family  curves  are  shown  in  Fig.  3(c).  I d  saturates 
for  Yds  when  Vds>~ 0.2  V.  The  superlinear  onset  is  also 
observed  and  the  so  called  Yds  threshold  voltage  increases 
at  lower  Yjg  [20].  A  peak  transconductances  of  ~4  mS//zm 
is  observed  around  Ytg= 0.12  V  (Fig.  3(d)),  which  are 
much  larger  than  ~0.8  mS//zm  reported  peak  transcon¬ 
ductances  of  10  nm  Fin-FET  [21].  In  Fig.  3(e),  the  top 
gate  changes  the  carrier  concentrations  of  the  top  2D  semi¬ 
conductor  much  faster  than  of  the  bottom  2D  semiconductor 
under  different  Yds-  The  ability  to  efficiently  change  a  hole 
(electron)  concentration  in  the  top  2D  semiconductor  while 
keeping  a  high  electron  (hole)  concentration  in  the  bot¬ 
tom  2D  semiconductor  is  vital  to  achieve  good  electrostatics 
control  of  these  Thin-TFETs.  The  quantum  capacitance  asso¬ 
ciated  with  the  top  and  bottom  semiconductor  layers  can 
be  expressed  as  Eq.  6: 


Cq,t(B) 


edpr(B )  +  edriT(B)  ^ 

_d(pp,T(B)  d<t>nJ{B)  _ 


(6) 


The  quantum  capacitances  are  plotted  in  Fig.  3(f)  under 
various  bias  conditions. 


FIGURE  4.  Effect  of  van  der  Waals  gap  thickness  variation  on  a  p- type 
Thin-TFET.  (a)  Tunnel  current  density  versus  VTC  for  different  van  der  Waals 
gap  thicknesses  Tvdw.  (b)  Differential  SS  versus  current  density  assuming 
an  evenly  distributed  van  der  Waals  gap  thickness  Tvdw  in  the  specified 
range. 


B.  EFFECTS  OF  NONUNIFORM  VAN  DER  WAALS  CAP 
THICKNESS  AND  ACCESS  RESISTANCE 

Due  to  the  nature  of  van  der  Waals  bonds,  the  van  der  Waals 
gap  thickness  is  subject  to  intercalation  of  atoms/ions, 
interlayer  rotational  misalignment  between  2D  layers  etc. 
For  instance,  in  bilayer  mechanically  stacked  Molybdenum 
Disulfide  (M0S2)  with  an  interlayer  twist,  a  maximum  vari¬ 
ation  of  0.59  A  [22]  was  experimentally  verified  in  the  van 
der  Waals  gap  thickness  [22].  Surface  roughening  due  to 
ripples  in  2D  crystals  or  roughness  of  the  underlying  sub¬ 
strates  can  also  introduce  van  der  Waals  gap  variations  [23]. 
Meanwhile,  tunneling  probability  is  very  sensitive  to  the 
tunneling  distance,  namely  the  van  der  Waals  gap  thickness 
in  a  Thin-TFET,  which  makes  it  important  to  investigate 
effects  of  a  non-uniform  van  der  Waals  thickness.  First,  the 
Thin-TFET  I-V  curves  are  calculated  by  varying  the  van  der 
Waals  gap  thickness  T vdw  from  3.0  A  to  6.0  A  and  a  step  of 
0.5  A  (which  is  roughly  half  of  the  Se  covalent  radius  [24]). 
The  results  are  shown  in  Fig.  4(a)  for  a  /7-type  Thin-TFET: 
the  on  current  density  decreases  and  the  threshold  voltage 
moves  towards  0  when  increasing  the  T vdw-  We  note  that, 
as  long  as  the  T vdw  is  uniform,  the  SS  remains  as  steep 
as  ~14  mV/dec.  However,  for  a  non-uniform  T vciw,  SS  will 
degrade.  To  estimate  its  impact,  an  evenly  distributed  T vdw 
over  several  ranges  is  used  in  the  calculated  differential  SS 
shown  in  Fig.  4(b).  For  example,  for  a  2D  heterojunction 
with  an  evenly  distributed  T vdw  from  3.0  A  to  5.0  A  and 
a  step  of  0.5  A,  we  take  the  corresponding  \d~Ytg  curve 
for  each  T vdw  (i.e.,  3.0  A,  3.5  A,  4.0  A,  4.5  A,  and  5.0  A) 
shown  in  Fig.  4(a)  and  average  them  over  the  T vdw  range 
to  obtain  the  overall  I d~Ytg  curve  for  the  calculation  of  SS. 
Fig.  4(b)  shows  that  up  to  1  A  variation  in  T vdw  is  tolerable, 
resulting  in  a  sub-60  mV/dec  SS  over  a  decent  current  win¬ 
dow  (up  to  50  /xA//xm).  Depending  on  how  Thin-TFETs  are 
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fabricated,  the  T vdw  non-uniformity  may  have  different  dis¬ 
tributions.  Our  first  look  at  its  impact  in  this  work  highlights 
the  importance  to  precisely  control  T vdw- 


FIGURE  5.  Effect  of  total  access  resistance  on  a  p-type  Thin-TFET. 

(a)  Iq  versus  VTC.  (b)  lD  versus  VDS  with  various  total  access  resistance  Rc 
values. 


A  finite  total  access  resistance  has  a  critical  impact  on 
ultrascaled  transistors.  To  date,  how  to  minimize  the  total 
access  resistance  in  2D  crystal  based  device  still  remains  an 
open  question.  In  Fig.  5,  we  show  its  effects  on  Thin-TFET 
by  assuming  several  values  for  the  total  access  resistance 
Rc-  At  a  sufficiently  high  |  V ds  I  of  0.4  V,  maximum  Id  is 
almost  the  same  for  a  Rc  of  up  to  320  Q/im;  a  higher 
Rc  decreases  maximum  Id  appreciably.  Understandably,  a 
lower  Rc  is  necessary  for  a  lower  Ydd •  In  an  ideal  2D 
conductor,  the  quantum  limit  of  the  total  access  resis¬ 
tance  is  inversely  proportional  to  the  square  root  of  the 
carrier  concentration;  e.g.,  ~52  Q/im  for  a  carrier  concen¬ 
tration  of  10 13  cm-2  [25].  Thus  the  access  region  of  2D 
semiconductors  can  be  degenerately  doped  to  minimize  Rc. 


C.  CAPACITANCE  EVALUATION 

The  gate-to-drain  and  gate-to-source  capacitances  (i.e.,  Cgd , 
Cgs )  can  be  readily  calculated  from  the  capacitance  network 
shown  in  Fig.  6. 

The  quantum  capacitances  Cqjj)  of  the  top  (bottom)  2D 
semiconductor  are  defined  in  Eq.  6  and  indicated  as  the  red 
non-linear  capacitances  in  Fig.  6.  First  we  define  C s  as: 


1  /Cs  =  1  /Cvdw  +  1  /(Cq,b  +  CBg )  (7) 


Then,  Cgd  and  Cgs  can  be  written  as  Eqs.  8: 


Cgs  = 


Cgd  = 


CtgCs 

Ctg  +  Cqj  +  Cs 
CtgCqj 
Ctg  +  Cqj  +  Cs 


(8) 


Due  to  the  symmetry  in  these  /7-type  and  n- type  Thin- 
TFETs  as  well  as  the  similar  hole  and  electron  effective  mass 
in  these  2D  crystals,  we  expect  similar  C-V  characteristics 
for  the  p-type  and  n- type  Thin-TFETs.  In  Fig.  7  we  plot  the 
calculated  C-V  curves  for  the  p-type  Thin-TFETs  shown  in 


FIGURE  6.  Capacitance  network  model  of  the  Thin-TFET. 


Fig.  2.  In  the  linear  region  of  the  I d~Y  ds  family  of  curves, 
Cgd  is  significant,  where  the  drain  is  coupled  with  the  top 
gate  to  modulate  the  tunnel  current.  From  the  linear  region 
to  the  saturation  region,  Cgd  drops  to  be  near  zero  while 
Cgs  increases  to  its  maximum.  What  is  worthy  noting  is  that 
the  magnitude  of  a  Thin-TFET  capacitance  is  smaller  than 
CMOS  and  III-V  TFET  benchmarked  in  Section  III-D  for  a 
given  gate  oxide  EOT  thus  capacitances,  which  stem  from 
the  serially  connected  capacitance  components  as  shown  in 
Fig.  6.  The  capacitance  model  is  useful  for  implementing 
the  Thin-TFET  into  circuit  simulations. 


FIGURE  7.  For  the  p-type  Thin-TFET.  (a)  CGD  and  Ccs  versus  VDS  at 
Vrc  =  —0.2,  -0.3,  -0.4  V.  (b)  CCD  and  Ccs  versus  Vrc  at  VDS  =  -0.2,  -0.3, 
-0.4  V. 

D.  BENCHMARKING 

The  Semiconductor  Research  Corporation  (SRC) 
Nanoelectronic  Research  Initiative  (NRI)  has  sup¬ 
ported  research  on  beyond  CMOS  devices  as  reported 
by  Bernstein  et  al.  [26]  As  part  of  the  initiative,  the 
projected  performance  of  the  beyond-CMOS  devices  and 
the  CMOS  of  the  same  technology  node  was  compared,  i.e., 
benchmarked.  The  benchmarking  activity  has  continued  by 
Nikonov  and  Young  [27],  [28].  Thin-TFET  being  proposed 
by  us  primarily  under  the  support  of  SRC  STARnet,  we 
participated  in  the  recent  benchmarking  using  the  Nikonov 
and  Young  (N&Y)  methodology. 

The  N&Y  methodology  uses  basic  device  performance 
parameters  such  as  operating  voltage  (Ydd  =  |VdsI),  satu¬ 
ration  current  (I Dsat),  and  average  gate  capacitance  (Cc,avg)i 
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to  project  logic  switching  energy  and  delay.  The  change 
of  the  net  charge  under  the  gate  (AQ=qAn?)  when  Ytg 
switches  from  0  to  Ydd  is  the  sum  of  the  change  of  the 
net  charge  in  the  top  2D  semiconductor  and  the  bottom 
2D  semiconductor.  The  average  gate  capacitance  ( Cc,avg )  is 
defined  as  AQ/Ydd-  Here  we  take  the  p-type  Thin-TFET 
as  an  example,  I Dsat  and  C c,avg  are  provided  in  Table  1  for 
a  few  Ydd  values  of  0.2,  0.3,  and  0.4  V  and  a  few  total 
access  resistance  R c  values  of  52  and  320  Qfi m.  The  device 
parameters  for  High  Performance  (HP)  CMOS,  Low  Power 
(LP)  CMOS,  InAs  Homojunction  TFET  (HomJTFET)  and 
InAs/GaSb  Heterojunction  TFET  (HetJTFET)  are  taken  from 
Ref.  [28]  and  we  use  the  same  geometrical  parameters  for 
all  the  devices  as  shown  in  Table  1,  while  neglecting  the 
contact  capacitance. 

The  intrinsic  switching  delay  t \nt  and  the  intrinsic  switch¬ 
ing  energy  E int  are  calculated  by  [28]: 

+  _  CG,avgVDD 

Hnt  —  T 

Wsat 

Eint  =  CG,avgWV2DD  (9) 


l102 
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FIGURE  8.  Intrinsic  switching  energy  and  delay  for  HP  CMOS,  LP  CMOS, 
HetJTFET,  HomJTFET,  and  Thin-TFETs  with  VDD  =  0.2,  0.3,  0.4  V,  and  Rc  =  52, 
320  S2/£m. 


Therefore,  given  that  C c,avg  stays  roughly  the  same  (increas¬ 
ing  slightly  with  decreasing  Ydd),  the  intrinsic  switching 
time  t int  slightly  decreases  when  decreasing  Ydd • 


In  Fig.  8,  we  plot  the  projected  values  of  t int  and  E int  of 
the  devices  listed  in  Table  1. 


TABLE  1.  Benchmarking  parameters. 


Parameters  for  Thin-TFETs  with  various  V  dd  and  Rc 

Vdd  (V) 

0.2 

0.3 

0.4 

Rc  (O/xm) 

52 

320 

52 

320 

52 

320 

iDsat  (jlA/lim) 

263 

233 

325 

317 

349 

348 

AQ  (fC//zm2) 

2.34 

2.80 

3.33 

3.72 

4.30 

4.47 

Ansxl0i2(/cm-2) 

1.46 

1.75 

2.08 

2.32 

2.69 

2.79 

Cc,avg  (fF//xm) 

0.175 

0.210 

0.167 

0.186 

0.161 

0.168 

Parameters  for  HP/LP  CMOS  and  HetJ/HomJ  TFET  [28] 

Vdd  (V) 

IDsat  (flA/fim) 

Cc,avg  (fF//im) 

HP  CMOS 

0.73 

1805 

1.29 

LP  CMOS 

0.3 

2 

1.29 

HetJTFET 

0.4 

500 

1.04 

HomJTFET 

0.2 

25 

1.04 

Geometrical  Parameters  for  Benchmarking 

Half-pitch 
(F)  (nm) 

EOT  (nm) 
(nm) 

Gate  Length 
(L)  (nm) 

Gate  Width 
(W)  (nm) 

15 

1 

15 

60 

As  far  as  the  intrinsic  switching  energy-delay  product  is 
concerned,  the  Thin-TFET  shows  distinct  energy  consump¬ 
tion  and  performance  advantages.  For  instance,  Thin-TFET 
operation  at  a  Ydd  as  low  as  0.2  V  is  fast  because  its  current 
is  still  significantly  high.  The  most  distinguishing  feature  of 
a  Thin-TFET  is  its  low  intrinsic  capacitance  in  comparison 
to  the  other  devices.  This  advantage  will  be  less  signifi¬ 
cant  when  device  parasitics  become  dominant  in  completed 
circuits. 

It  is  observed  that  the  Thin-TFET  intrinsic  switching 
energy-delay  product  moves  toward  the  desired  corner  when 
decreasing  Ydd  from  0.4  V  to  0.2  V.  This  is  an  unusual  but 
favorable  behavior  for  ultrascaled  switches.  In  the  case  of 
15  nm  CMOS,  Id  is  roughly  proportional  to  Ydd •  While  in 
the  ON  state  of  Thin-TFET,  I d  has  much  weaker  dependence 
on  Ytg  (see  Fig.  5(a))  than  CMOS,  thus  Ydd  to  I d  ratio  actu¬ 
ally  decreases  when  scaling  down  Ydd  from  0.4  V  to  0.2  V. 


£.  EXPERIMENTAL  INSIGHTS 

Since  our  proposal  of  Thin-TFET  in  2012  [29]  that  is 
derived  from  our  III-V  TFET  design  [1],  several  key  chal¬ 
lenges  have  been  identified  along  our  pursuit  in  experimental 
demonstration  of  Thin-TFETs  [30].  The  foremost  is  the 
scarcity  of  electronic-grade  layered  materials  and  knowl¬ 
edge  of  their  properties,  in  particular,  the  semiconductor 
heterojunctions  with  near  broken  gap  alignment.  The  reason¬ 
ably  well-characterized  material  properties  in  the  literature 
are  largely  based  on  bulk  layered  materials.  An  exponen¬ 
tially  growing  number  of  publications  in  the  recent  years 
on  monolayer  and  few-layer  materials  are  mainly  theoretical 
calculations  or  based  on  exfoliation  of  naturally  occurring 
crystals  or  synthesized  by  chemical  vapor  transport,  which 
typically  contains  a  few  atomic  percent  of  defects  (impu¬ 
rities,  vacancies  etc).  Both  chemical  vapor  deposition  and 
molecular  beam  epitaxy  [31]  are  actively  pursued  by  the 
community  to  grow  electronic  grade  layered  materials. 

Besides  lack  of  high  quality  layered  materials  and  het- 
eroj unctions,  the  fabrication  development  of  Thin-TFET  is 
also  challenging.  It  inherits  all  the  fundamental  fabrication 
challenges  of  a  TFET  including  doping  profile,  alignment 
especially  gate  registry,  gate  dielectrics,  ohmic  contacts. 
Atomic  layer  deposition  has  been  improved  over  years  to 
achieve  good  quality  gate  dielectrics  on  2D  crystals  [32]. 
Using  2D  dielectrics  such  as  hexagonal  boron  nitride  as 
the  gate  dielectrics  has  also  been  pursued  [33].  Third,  low 
resistance  ohmic  contacts  to  2D  crystal  are  vital  to  device 
performance.  Various  techniques  such  as  external  chemical 
doping  [34],  internal  chemical  doping  [35],  electrostatic  dop¬ 
ing  such  as  ion  doping  [36]  and  phase-engineering  from 
the  semiconductor  phase  to  the  metallic  phase  of  a  2D 
crystal  [37],  have  been  implemented  to  reduce  the  contact 
resistances.  Furthermore,  Thin-TFETs  demand  true  preci¬ 
sion  layer  number  control  since  the  properties  of  nearly  all 
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layered  materials  critically  depend  on  the  layer  number  when 
the  layer  number  is  in  the  range  of  1-3  nm. 

IV.  CONCLUSION 

A  new  tunnel  transistor,  Thin-TFET,  has  been  proposed  and 
a  model  material  system  identified.  Simulations  based  on 
the  transfer  Hamiltonian  method  suggest  that  Thin-TFETs 
can  achieve  desired  sub-threshold  swing  (SS)  and  high  on- 
current.  A  uniform  van  der  Waals  gap  thickness  and  low 
total  access  resistance  are  vital  to  optimize  the  Thin-TFET 
performance.  The  benchmark  study  shows  Thin-TFETs  may 
have  distinct  advantages  over  CMOS  and  III-V  TFETs  in 
term  of  both  performance  and  energy  consumption  at  low 
supply  voltages. 
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We  report  the  realization  of  field-effect  transistors  (FETs)  made  with  chemically  synthesized 
multilayer  2D  crystal  semiconductor  MoS2.  Electrical  properties  such  as  the  FET  mobility, 
subthreshold  swing,  on/off  ratio,  and  contact  resistance  of  chemically  synthesized  (s-)  MoS2  are 
indistinguishable  from  that  of  mechanically  exfoliated  (x-)  MoS2,  however,  flat-band  voltages  are 
different,  possibly  due  to  polar  chemical  residues  originating  in  the  transfer  process.  Electron 
diffraction  studies  and  Raman  spectroscopy  show  the  structural  similarity  of  s-MoS2  to  x-MoS2. 

This  initial  report  on  the  behavior  and  properties  of  s-MoS2  illustrates  the  feasibility  of  electronic 
devices  using  synthetic  layered  2D  crystal  semiconductors.  ©  2013  American  Institute  of  Physics. 
[http://dx.doi.org/10.1063/L4789975] 


Two-dimensional  (2D)  crystal  materials  are  receiving 
increased  attention  for  future  electronic  devices.  Short- 
channel  effects  in  modern  transistors  originate  from  the 
3-dimensional  (3D)  nature  of  the  gate  control;  introduction 
of  2D  materials  in  the  channel  significantly  improves  the 
gate  electrostatics.  Graphene  is  a  true  2D  material,1  but  its 
lack  of  a  bandgap  results  in  high  leakage  in  the  off-state  in 
conventional  transistor  geometries.  Alternative  switching 
mechanisms2-4  could  enable  electronic  switching  with  high 
on-off  current  ratios.  2D  transition-metal  dichalcogenide 
(TMD)  materials  such  as  MoS2,5’6  WSe2,7’8  and  WS2 
(Refs.  9  and  10)  have  drawn  considerable  attention  due  to 
the  presence  of  a  bandgap,  in  contrast  to  graphene.  Prior 
studies  of  TMD  multilayered  materials  have  been  conducted 
using  exfoliated  layers,5-9  similar  to  the  initial  work  with 
graphene.1  The  size  and  quality  of  naturally  occurring  exfoli¬ 
ated  layered  semiconductor  materials  are  limited  and  uncon¬ 
trollable,  so  it  is  important  to  develop  synthetic  techniques. 
Multilayer  TMDs,  like  graphite,  are  excellent  solid  lubricat¬ 
ing  materials  and  have  been  chemically  synthesized  in  large 
volumes,  but  have  not  been  investigated  intensively  for  tran¬ 
sistors11-13  except  in  very  recent  work.14  Here,  we  report  the 
fabrication  and  demonstration  of  chemically  synthesized 
multilayer  (s-)  MoS2  field-effect  transistors  (FETs)  and  com¬ 
pare  the  properties  with  exfoliated  (x-)  MoS2  FETs. 

s-MoS2  flakes  were  grown  by  an  iodine-transport 
method  from  previously  synthesized  MoS2  (0.6  g)  at  1060  K 
in  an  evacuated  silica  ampoule  at  a  pressure  of  10-3  Pa,  and 
with  temperature  gradient  of  6.8  K/cm.  The  volume  concen¬ 
tration  of  iodine  was  11  mg/cm3.  After  21  days  of  growth, 
the  silica  ampoule  was  slowly  cooled  to  room  temperature  at 
a  rate  of  30°C/h.  The  MoS2  flakes  were  then  dispersed  by 


a) Electronic  addresses:  whwang@us.ibm.com  and  whwangl@nd.edu. 
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sonication  in  isopropyl  alcohol  (IP A)  and  transferred  onto  a 
30  nm  thick  atomic -layer-deposited  (ALD)  A1203  dielectric 
on  a  p  +  Si  substrate  held  at  100  °C  until  dry.  For  compari¬ 
son,  x-MoS2  flakes  were  released  and  transferred  from  bulk 
MoS2  using  scotch  tape.  The  height  of  the  MoS2  flakes  are  in 
the  range  of  2CM0  nm.  Source  and  drain  contacts  were 
defined  by  electron  beam  lithography  (EBL)  using  Ti/Au 
(5/100 nm)  contacts.  The  devices  were  annealed  at  300  °C 
for  3  h  under  Ar/H2  flow  to  decrease  the  contact  resistance. 
A  schematic  cross-sectional  image  of  the  back-gated  (BG) 
MoS2  device  is  shown  in  the  inset  of  Fig.  1(a). 

Figure  1(a)  shows  the  measured  drain  current,  ID ,  per 
unit  gate  width,  versus  the  back-gate-to-source  voltage,  VBg> 
at  room  temperature  for  a  multilayer  s-MoS2  channel  at  three 
drain  biases.  The  gate  modulation  is  ~105  and  the  gate 
leakage  current  is  much  lower  (less  than  1  pA/^um)  than  the 
drain  current.  This  large  gate  modulation  relative  to  gra¬ 
phene  is  attributed  to  the  presence  of  a  bandgap.  The  device 
shows  clear  ft-type  behavior  indicating  accumulation  of 
electrons  (n- type  conductivity)  for  positive  back-gate  bias. 
The  comparative  transfer  curves  of  x-MoS2  FETs  are  also 
shown  in  Fig  1(a).  s-MoS2  and  x-MoS2  FETs  show  highly 
similar  transfer  characteristics.  The  extracted  field-effect 
mobilities  of  both  x-MoS2  and  s-MoS2  FETs  are  ^30  cm2/V  s 
at  room  temperature.  The  subthreshold  swing  (SS)  of  the 
x-MoS2  FET  is  200mV/dec.  and  that  of  the  s-MoS2  FET  is 
190mV/dec.  The  average  values  of  the  FET  mobility  of  each 
type  of  both  x-MoS2  and  s-MoS2  are  ^15cm2/V  s  and 
the  SS  of  those  are  ^170mV/dec.  Another  s-MoS2  (W/L 
=  1/2  pm)  FET  is  compared  with  x-MoS2  and  electrical 
properties  still  work  out  to  be  similar. 

The  subthreshold  swing  is  similar  for  the  s-MoS2  and 
x-MoS2  FETs,  but  higher  than  the  ideal  Boltzmann  limit  of 
60mV/dec.  The  similar  SS  suggests  that  the  interface  charge 
leading  to  the  higher  subthreshold  swing  likely  arises  from 
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FIG.  1.  Transport  properties  and  Raman  spectroscopy  of  FETs  of  s-MoS2  with  W/L  =  1/2  /an  and  x-MoS2  with  W/L  =  5/1.5  /an.  (a)  Drain  current,  IDr  per  unit 
gate  width  vs.  back-gate  voltage,  VBG,  of  s-MoS2  at  various  drain  voltages,  VDS.  The  transfer  curve  of  an  x-MoS2  is  also  shown  for  comparison,  dotted  line, 
(b)  Common  source  transistor  characteristics  comparing  s-MoS2  (solid  lines)  and  x-MoS2  (dashed  lines)  FETs.  (c)  Raman  spectra  ( Xexc  —  488  nm)  of  both 
s-MoS2  and  x-MoS2  materials  with  a  laser  power  of  1.5  mW  and  a  spot  size  of  0.5-1  /an.  The  inset  sketch  shows  the  two  primary  vibrational  modes  in  MoS2 
leading  to  the  two  peaks  in  the  Raman  spectrum. 


traps  in  the  A1203  dielectric,  which  are  identical  for  the 
s-MoS2  and  x-MoS2  FETs.  The  major  difference  between 
s-MoS2  and  x-MoS2  FETs  is  in  the  threshold  voltage.  The 
flat-band  voltage  of  the  x-MoS2  FETs  is  higher  than  that  of 
s-MoS2  by  approximately  1  V.  The  possible  reasons  for  this 
shift  could  be  (a)  different  unintentional  doping  densities  in 
s-MoS2  and  x-MoS2,  or  (b)  scotch  tape  residue-induced- 
charges  possibly  leading  to  a  higher  flat-band  voltage  for  the 
x-MoS2  compared  to  the  s-MoS2  FETs,  which  do  not  experi¬ 
ence  the  tape  exfoliation  procedure.  Identifying  the  precise 
reason  for  this  threshold  shift  requires  further  work. 

The  family  of  Id~Vds  curves  at  various  VGS  in  Fig.  1(b) 
shows  typical  transistor  behavior  including  a  linear  increase 
of  current  at  low  VDS  and  current  saturation  at  high  VDS.  The 
behavior  shows  desirable  transistors  attributes  such  as  ohmic 


contacts,  current  saturation,  and  good  gate  electrostatic  con¬ 
trol.  However,  the  current  levels  are  in  the  /iA//im  regime  for 
micron  long  gate-lengths,  which  is  low.  The  reason  for  the 
low  current  is  a  high  contact  resistance.  The  contact  resistan¬ 
ces  of  both  s-MoS2  and  x-MoS2  FETs  in  this  work  were 
extracted  to  be  ~80  Q  mm  at  low  VDS •  This  value  is  compa¬ 
rable  to  the  ^69  reported  for  exfoliated  MoS2  FETs.6  This  is 
an  extremely  high  value,  and  currently  holds  back  the  per¬ 
formance  of  the  FETs.  The  observation  here  is  that  s-MoS2 
and  x-MoS2  FETs  have  similar  contact  resistances,  and  are 
both  high.  A  significant  increase  in  the  current  drive  of  the 
FETs  is  expected  if  the  contact  resistance  can  be  lowered  to 
1  Qmm  regime  or  lower,  as  is  the  case  in  Si  and  IIT-V  FETs. 
The  contact  resistance  can  be  lowered  using  a  low  work 
function  metal  like  Sc.15 
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FIG.  2.  (a)  TEM  of  the  s-MoS2  multi¬ 
layers.  TEM  electron  diffraction  patterns 
from  (b)  region  I  revealing  the  single 
crystal  layer,  (c)  region  II  revealing  the 
superposition  of  two  single  crystal  layers 
with  rotation,  and  (d)  region  III  reveal¬ 
ing  the  superposition  of  four  single  crys¬ 
tal  layers  with  different  rotations.  The 
electron  diffraction  patterns  reveal  that 
the  lattice  parameter  of  s-MoS2  is 
0.32  nm.  (e)  An  atomic-scale  Moire  pat¬ 
tern  of  s-MoS2  multilayers,  (f)  FFT  from 
the  image  of  (e),  indicating  three  differ¬ 
ent  crystal  layers  with  different  rota¬ 
tional  angles. 
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FIG.  3.  (a)  TEM  of  the  x-MoS2  multi¬ 
layers.  (b)  Electron  diffraction  pattern 
from  region  I  revealing  the  single  crystal 
layer,  (c)  HR  TEM  image  of  region  I.  A 
lattice  constant  of  0.32  nm  can  be 
obtained  from  (b)  and  (c). 


The  measured  Raman  spectra  shown  in  Fig.  1(c)  at  an 
excitation  wavelength,  Xexc,  of  488  nm  are  highly  similar 
between  s-MoS2  and  x-MoS2  flakes.  The  spectrum  exhibits 
two  peaks:  one  in  the  E2g1  range,  corresponding  to  in-plane 
vibrations  at  ^385  cm-1,  and  the  other  in  the  Alg  range  cor¬ 
responding  to  out-of-plane  vibrations  at  ^410  cm-1.  Since 
the  Raman  spectrum  and  transport  properties  are  almost 
identical  for  s-MoS2  and  x-MoS2,  this  confirms  the  overall 
similarity  of  chemically- synthesized  and  exfoliated  MoS2.  In 
order  to  further  investigate  the  atomic  properties  of  s-MoS2, 
transmission  electron  microscopy  (TEM)  was  performed. 
The  information  obtained  from  s-MoS2  is  shown  in  Fig.  2 
and  compared  with  that  of  x-MoS2  in  Fig.  3. 

Figure  2(a)  shows  the  morphology  of  the  s-MoS2  on  a 
TEM  grid  where  a  single  crystal  layer  region  (I),  a  superposi¬ 
tion  of  two  single  crystal  layer  with  rotation  (II),  and  super¬ 
position  of  four  single  crystal  layers  with  rotation  (III)  are 
clearly  resolved.  The  corresponding  TEM  electron  diffrac¬ 
tion  patterns  are  shown  in  Figs.  2(b)-2(d),  respectively, 
showing  the  number  of  spots  corresponding  to  the  number  of 
layers  and  the  angular  separation  superimposed  on  the  hex¬ 
agonal  lattice  pattern.  The  electron  diffraction  patterns  in 
Figs.  2(b)-2(d)  reveal  that  s-MoS2  flakes  retain  the  crystal 
symmetry  and  lattice  constant,  ^0.32  nm,  when  compared 
with  that  of  x-MoS2  in  Fig.  3.  A  high-resolution  (HR)  TEM 
image  is  shown  in  Fig.  2(e)  and  its  fast  Fourier  transform 
(FFT)  is  shown  in  Fig.  2(f).  The  FFT  shows  that  the  Moire 
pattern  in  the  image  is  due  to  three  different  crystal  layers 
superimposed  with  different  angles. 

For  comparison,  the  morphology  of  x-MoS2  on  the  TEM 
grid  and  the  electron  diffraction  pattern  with  a  high- 
resolution  TEM  image  are  shown  in  Figs.  3(a)-3(c),  respec¬ 
tively.  The  measurement  confirms  again  a  lattice  parameter 
of  x-MoS2  of  ^0.32  nm  which  is  identical  to  that  of  s-MoS2 
and  that  of  bulk  MoS2.16  The  single-layer  region  is  clearly 
resolved  in  the  electron  diffraction  pattern,  and  the  high- 
resolution  atomic  image  clearly  resolves  the  hexagonal  crys¬ 
tal  structure  of  single-layer  MoS2. 

In  summary,  chemically  synthesized  MoS2  transistors 
were  fabricated  and  characterized,  and  compared  with  exfoli¬ 
ated  MoS2.  The  electronic  and  structural  properties  of  chemi¬ 
cally  synthesized  layers  were  found  to  be  highly  similar  to 
exfoliated  MoS2.  In  particular,  the  transistor  characteristics 
of  s-MoS2  were  found  to  be  almost  identical  to  that  of  the 
x-MoS2  in  terms  of  FET  mobility,  subthreshold  swing,  on/off 
ratio,  and  contact  resistance.  TEM  electron  diffraction 
patterns  and  Raman  measurements  prove  that  the  crystal 
symmetries  and  structural  properties  of  s-MoS2  are  also 


identical  to  that  of  x-MoS2.  Though  a  number  of  issues  need 
to  be  resolved  before  TMD  crystals  deliver  high-performance 
transistors,  this  initial  report  shows  that  synthetic  procedures 
can  also  realize  high-quality  channel  material. 
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Interband  quantum  tunneling  of  electrons  in  semiconductors  is  of  intense  recent  interest  as  the 
underlying  transport  mechanism  in  tunneling  field-effect  transistors.  Such  transistors  can 
potentially  perform  electronic  switching  with  lower  energy  than  their  conventional  counterparts. 
The  recent  emergence  of  two-dimensional  (2D)  semiconducting  crystals  provides  an  attractive 
material  platform  for  realizing  such  devices.  In  this  work,  we  derive  an  analytical  expression  for 
understanding  tunneling  current  flow  in  single-layer  2D  crystal  semiconductors  in  the  k- space.  We 
apply  the  results  to  a  range  of  2D  crystal  semiconductors,  and  compare  it  with  tunneling  currents  in 
three-dimensional  semiconductors.  We  also  discuss  the  implications  for  tunneling  devices.  ©  2013 
American  Institute  of  Physics.  [http://dx.doi.org/10.1063/L4799498] 


Two-dimensional  (2D)  crystals  of  graphene  were  first 
isolated  in  2004. 1  The  unique  electronic  and  optical  proper¬ 
ties  of  graphene  have  been  extensively  studied  since  then.2 
Electronic  transistors  have  been  proposed  with  graphene  for 
ultra-low  power  switching.3-5  These  proposed  devices  either 
exploit  the  symmetric  zero-bandgap  Dirac-cone  bandstructure 
of  graphene  or  require  the  opening  of  energy  bandgaps  by 
quantum  confinement.  Soon  after  the  isolation  of  graphene, 
semiconducting  2D  crystals  were  reported.6-8  2D  crystal 
semiconductors  have  been  found  in  the  transition  metal 
dichalcogenide  (TMD)  material  family,  and  the  list  is 
expected  to  expand  in  the  future.9  Taking  advantage  of  an 
energy  bandgap,  ultrathin  channels,  and  absence  of  broken 
bonds,  conventional  field-effect  transistors  (FETs)  using  2D 
crystal  semiconductors  have  shown  high  promise  in  initial 

.  .  .  7  in 

investigations.  ’ 

A  number  of  electronic  switching  devices  have  been  pro¬ 
posed  recently  to  address  the  power-dissipation  problems  of 
FETs.11  Among  the  proposed  devices,  the  tunneling  FET  (or 
TFET)  has  emerged  as  an  attractive  candidate.  These  devices 
take  advantage  of  the  potential  of  interband  Zener-tunneling 
of  electrons  to  beat  the  Boltzmann  thermal  limit  of  switching 
of  60  mV/decade.  Initial  experimental  demonstrations  show 
much  promise.12-14  The  availability  of  2D  semiconducting 
crystals  with  bandgaps  begs  the  question  whether  TFETs 
with  attractive  properties  can  be  realized  with  them.  Such  2D 
crystal  TFETs,  if  realized,  can  take  advantage  of  the  ultrathin 
nature  of  the  layers,  and  the  absence  of  broken  bonds  to  ena¬ 
ble  scaling  of  such  devices  to  much  smaller  dimensions  than 
three-dimensional  (3D)  crystal  semiconductors.  A  critical 
metric  for  TFETs  is  the  on- state  current,  which  is  limited  by 
interband  tunneling  of  electrons.  To  date,  interband  tunneling 
in  purely  2D  semiconducting  crystal  junctions  has  not 
received  sufficient  attention,  certainly  not  to  the  extent  it  has 
for  3D  semiconductor  p-n  junctions  since  Zener’s15  and 
Esaki’s  works.16  In  this  work,  we  solve  this  problem.  We 
derive  an  analytical  expression  for  the  tunneling  current  in 
2D  crystal  semiconductors.  The  expression  highlights  the  de¬ 
pendence  of  the  tunneling  current  on  the  material  parameters 
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of  the  2D  crystal  semiconductor,  such  as  its  bandgap  and 
effective  masses.  We  apply  the  results  to  a  range  of  2D  crys¬ 
tals,  and  discuss  the  implications  for  device  applications. 

Consider  the  2D  crystal  p-i-n  junction  shown  schemati¬ 
cally  in  Fig.  1(a)  with  ohmic  contacts  to  the  p-  and  zz-doped 
regions.  The  contacts  would  form  the  source  and  drain  con¬ 
tacts  of  the  corresponding  TFET.  We  do  not  address  the  ex¬ 
perimental  challenges  of  doping  and  electrostatic  gating  in 
this  work  and  focus  exclusively  on  evaluating  the  two- 
terminal  tunneling  current.  Assume  the  doping  in  the  p-  and 
n- sides  aligns  the  Fermi  levels  to  the  respective  band-edges. 
Then,  under  no  applied  bias,  Epv  =  Ef  and  no  net  current 
flows  across  the  junction.  Here,  Ep  is  the  valence-band  edge 
on  the  p-side  and  Enc  is  the  conduction-band  edge  on  the  n- 
side.  Under  the  application  of  a  reverse  bias  voltage  V,  a  finite 
energy  window  is  created  for  electrons  since  Ep  —  Enc  =  qV. 
Within  this  energy  window,  electrons  from  the  valence  band 
can  tunnel  into  the  conduction  band  on  the  other  side,  as  indi¬ 
cated  in  Fig.  1(b). 

The  electric  current  is  obtained  by  summing  the  individ¬ 
ual  quantum-mechanical  probability  current  contributions  by 
each  k- state  electron,  and  multiplying  it  by  q ,  the  electron 
charge.  The  tunneling  current  is  thus  given  by 

h  ~fc)Twkb,  (1) 

k 

where  gs  =  2  is  the  spin  degeneracy  and  gv  is  the  valley 
degeneracy  for  2D  crystal  single  layers.  Single-layer  TMD  2D 
crystals  have  been  found  to  have  a  direct  bandgap  with  two 
valleys  at  the  K  and  K'  points  of  the  Brillouin  zone  similar  to 
2D  graphene,17  so  we  use  gv  =  2.  Lx  is  the  macroscopic 
length  along  the  electric  field  (which  will  cancel  out),  vg(k)  = 
H~l\7E(k)  is  the  group  velocity  of  carriers  in  the  band 
£(/:), /v,  fc  are  the  occupation  functions  of  the  valence  and 
conduction  bands,  respectively,  and  Twkb  is  the  interband  tun¬ 
neling  probability  given  by  the  Wentzel-Kramers-Brillouin 
(WKB)  approximation.  The  sum  is  over  all  state  electrons 
that  are  allowed  to  tunnel.  We  note  here  that  the  tunneling 
probability  may  be  obtained  by  various  means  such  as  the 
Landau-Zener  approach,  Bardeen’s  transfer  Hamiltonian,  or 
direct  numerical  evaluation  by  integration  over  evanescent 
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FIG.  1.  Schematic  depiction  of  a  2D  crystal  p-i-n  junction  (a),  the  energy 
band-diagram  (b),  and  the  k-space  distribution  of  current  densities  (c).  The 
k-states  contributing  to  interband  tunneling  current  and  the  group  velocity 
are  indicated. 


line  move  to  the  right  in  the  +kx  direction,  as  indicated  in 
the  semi-circle  in  the  &-space  in  Fig.  1(c).  Since  there  are 
negligible  electrons  in  the  conduction  band  in  that  line,  the 
current  there  must  be  carried  by  electrons  in  the  valence 
band.  Which  of  these  right-going  electrons  are  allowed  to 
tunnel  through  the  gap?  In  the  absence  of  phonon  scattering, 
tunneling  is  an  elastic  process.  This  enforces  the  energy  con¬ 
servation  requirement 


+  n) 


(3) 


with  the  additional  requirement  that  the  lateral  momentum 
be  conserved,  i.e.,  kyp  =  kyn  =  ky.  The  energy  and  momen¬ 
tum  conservation  requirements  thus  lead  to  the  relation 


a  |  <  7.2  =  2mUv  K  i  i 

xp  rrfR  y  H2  m*  xn‘ 


(4) 


Let  us  define  £^ax  =  2 m*vqV /ft2  and  r\2  =  m*v/mR.  Note 
that  kmSLX ,  the  radius  of  the  semi-circle  in  the  k- space  shown 
in  Fig.  1(c),  is  controlled  by  the  applied  voltage.  Since  there 
is  an  electric  field  in  the  x-direction,  the  momentum  in  that 
direction  will  not  be  conserved.  For  the  electron  to  emerge 
on  the  right  (ft-)side,  kxn  must  be  non-zero,  and  thus  k2n  >  0, 
which  implies 


k2  +  n2k2  <  k2 

^xp  '  '  —  ^ma 


(5) 


k-states  in  the  bandgap.  Since  the  results  from  the  analytical 
WKB  approach  are  not  vastly  different  from  the  alternate 
approaches,  we  use  this  approach. 

The  tunneling  probability  T wkb  is  obtained  by  the  WKB 
approximation  in  the  following  manner.  For  2D  crystals, 
electrons  in  the  valence  band  of  the  /7-side  have  a  transverse 
kinetic  energy  Ey  =  h2k2/2m*v ,  where  hky  is  the  transverse 
quasi-momentum  and  m*  is  the  valence  band  effective  mass. 
The  WKB  tunneling  probability  is  then  given  by18 


The  above  condition  defines  a  restricted  elliptical  area  AT  of 
the  k- space  semi-circle  for  electron  states  that  are  allowed  to 
tunnel,  as  shown  in  Fig.  1(c).  We  can  now  evaluate  the  tunnel¬ 
ing  current  for  2D  semiconductor  p-i-n  junctions.  In  the 
expression  for  the  tunneling  current  (Eq.  (1)),  the  group  veloc¬ 
ity  term  is  that  of  the  valence  band  k- state  vg(k)  =  hkx/ml. 
We  skip  the  p- subscripts,  since  it  is  clear  that  the  electrons 
tunnel  from  the  valence  band  of  the  /7-side.  The  expression  for 
the  tunneling  current  is  then 


Twkb  =  exp 


T0  exp 


j  gsgv  hkx 


(fv  -/c)r0ex p 


n24 ' 

2m*E 


(6) 


where  T0  =  exp  [—4  yj 2 mRE3J2  /  3 qhF]  is  the  tunneling  proba¬ 
bility  of  perpendicularly  incident  electrons,  E  =  qhF / 2 
yj2 mREg,F  is  the  (constant)  electric  field  in  the  junction, 
and  mR  is  the  reduced  effective  mass  given  by  mR  = 

/(m*  +m*).  are  the  effective  masses  of  electrons  of 

the  valence  and  conduction  bands,  respectively.  The  above 
expression  is  found  to  be  consistent  with  experimental 
results.18  Note  that  the  tunneling  probability  of  electrons  is 
lowered  exponentially  with  their  transverse  kinetic  energy  as 
a  consequence  of  lateral  momentum  conservation  in  the  tun¬ 
neling  process.  We  neglect  phonon  emission  or  absorption 
processes  here.  To  evaluate  the  tunneling  current,  we  attach 
this  tunneling  probability  to  each  electronic  Estate  and  sum 
it  over  all  electrons  incident  on  the  tunneling  barrier. 

In  Fig.  1,  we  concentrate  on  a  particular  ID  line  as 
shown  by  the  dashed  line,  at  the  p-i  junction,  which  is  the 
source  side.  Half  of  the  electrons  in  the  valence  band  in  that 
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The  sum  over  £- states  is  converted  into  an  integral  via  the 
standard  recipe  ^LxLy/(2n)2  x  J dkxdky (...).  Due 

to  the  "filtering"  brought  about  by  the  requirements  of  energy 
and  momentum  conservation  (Eq.  (5)),  the  &-space  integral  is 
evaluated  over  the  restricted  area  Aj.  The  tunneling  current 
per  unit  width  or  the  current  density  is  then  given  by 


j'id jjL  VgsgyhTo 

T  Ly  (27 i)2m* 


T&max/  Y] 


dkv 


x  exp 


n24 
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Electrons  incident  normal  to  the  junction  have  no  transverse 
momentum,  and  carry  most  of  the  tunneling  current.  The 
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number  of  electron  states  allowed  to  tunnel  reduces  as  their 
transverse  directed  momentum  increases,  as  shown  sche¬ 
matically  in  Fig.  1(c).  The  current  carried  by  these  states 
with  transverse  momentum  is  further  damped  by  the 
[—Ey/E]  factor,  leading  to  further  filtering  and  momentum 
collimation.  The  maximum  tunneling  current  is  carried  by 
states  closest  to  (kx,ky)  =  (km ax,0)  as  indicated  by  the  shad¬ 
ing  in  the  figure. 

To  evaluate  the  current,  the  integral  in  k- space  should  be 
evaluated.  At  T  —>  OK ,/v  —fc~  1,  for  the  energy  window 
of  current-carrying  electrons.  This  relation  remains  an  excel¬ 
lent  approximation  at  room  temperature.  The  interband  tun¬ 
neling  current  per  unit  width  (pA/pm)  in  a  2D  crystal  p-i-n 
junction  then  evaluates  to 


r2 D  _<f( gsgyTp  \  fl m*E 
h  \  2n  )  y  fi2 


X 


where  Erf[...]  stands  for  the  error  function,  and  we  have 
defined  Vo  =  r\2E/q.  Equation  (8)  is  the  central  result  of  this 
work.  The  expression  shows  the  dependences  on  various 
bandstructure  and  junction  parameters  explicitly.  For  small 
reverse  bias  voltages  V  Vo,  the  tunneling  current  varies  as 
JjD  ~  y3/2  to  leading  order.  This  is  consistent  with  a  recent 
report  investigating  dimensionality  effects  on  tunneling.19 
For  larger  voltages  when  V  Vo,  Erf  [. . .]  — >  1,  and  we  get  a 
linear  dependence  of  the  tunneling  current  on  the  voltage 


t2D 


<72  fgsgy) 

j2nml 

qHF 

x  Tn 

h  \2nJ 

V  * 

^8m*REg  \ 
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where  the  Landauer  quantum  of  conductance  is  split  off. 
The  entire  square  root  term  is  an  effective  wavevector  with 
units  of  inverse  length,  leading  to  units  of  current  per  unit 
width.  The  corresponding  current  density  for  tunneling  cur¬ 
rent  in  3 -dimensional  semiconductors  is  given  by  JjD 


<r  (Mi) 

h  V  2n  ) 


2 mREg  qF 
i2  Et 


x  To 


V  with  units  of  current  per 


unit  area.  The  WKB  term  is  similar  for  2D  and  3D  crystals. 


The  prefactor  for  3D  semiconductors  goes  as  rsj  Fy/m* /Eg, 
whereas  for  2D  crystals,  it  goes  as  its  square  root 

F yjnf  /Eg.  In  quasi-2D  systems,  multiple  subbands 
may  be  involved  in  transport.  Then,  we  sum  the  current 
from  each  subband  with  the  respective  band  parameters. 
When  the  temperature  T  is  high,  the  assumption  fv  —fc~  1 
may  no  longer  be  suitable.  In  that  case,  the  2D  tunneling 
current  becomes 
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(10) 


where  f}  =  1/ksT  and  kB  is  the  Boltzmann  constant.  We 
have  not  simplified  this  expression  analytically,  but  the  nu¬ 
merical  evaluation  is  discussed.  The  interband  tunneling  cur¬ 
rent  densities  of  various  2D  crystals  at  T  =  4  K  and  300  K  are 
plotted  as  solid  and  dashed  lines  in  Fig.  2(a),  respectively. 
As  is  evident,  the  temperature  dependence  is  rather  weak. 
The  material  constants  (bandgaps  and  effective  masses)  are 
obtained  from  Refs.  17  and  20.  We  note  that  the  field  F  is 
determined  by  the  choice  of  doping,  bandgap,  external 
dielectrics,  and  in  the  case  of  a  3 -terminal  TFET  geometry, 
the  gate  voltage. 

The  tunneling  current  densities  for  MoS2  and  the  family 
of  TMDs  are  found  to  be  low  owing  to  their  large  bandgaps. 
For  example,  the  current  density  approaches  ~0.1  pA/pm  for 
MoTe2  at  a  high  field  of  4  MV/cm.  The  tunneling  current 
density  of  2D  graphene  can  be  higher  (~few  mA/pm),  but  it 
lacks  a  bandgap.  For  TFET  applications,  2D  crystals  with 
smaller  bandgaps  are  necessary  for  boosting  the  current.  For 
example,  tunneling  currents  for  2D  crystals  with  bandgaps  of 
0.5  eV  and  1.0  eV  with  corresponding  lower  effective  masses 


F  (MV/cm) 


Voltage  (V) 


FIG.  2.  (a)  Interband  tunneling  current  den¬ 
sity  for  various  2D  crystal  semiconductors 
at  a  reverse  bias  of  V  =  0.3V.  The  solid 
lines  are  at  T  =  4  K,  and  the  dashed  lines  at 
T  =  300  K,  and  the  temperature  dependence 
is  weak,  (b)  Current-voltage  curves  at  vari¬ 
ous  temperatures  at  a  junction  field  F  =  4 
MV/cm  for  a  2D  crystal  semiconductor  with 
band  parameters  indicated,  (c)  Same  as  (b), 
but  for  the  2D  crystal  MoTe2. 
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FIG.  3.  (a)  Interband  tunneling  currents  for  2D 
crystal  semiconductors  as  a  function  of  the 
energy  bandgap  for  various  effective  mass  pa¬ 
rameters.  (b)  The  high  current  part  of  (a) 
zoomed  in  for  more  details. 


are  plotted  in  Fig.  2(a).  The  currents  for  such  crystals  exceed 
^100  jiA/jim  at  the  highest  junction  fields,  and  thus  can  be 
attractive  for  high-performance  TFET  applications.  Such 
small-bandgap  materials  could  be  intrinsic  2D  crystals  or 
derived  from  interaction-induced  bandgap  of  Dirac-cone  sur¬ 
face  states  in  thin  topological  insulator  materials.21  Another 
possibility  is  in  bilayer  graphene,  where  breaking  the  layer 
symmetry  by  vertical  electric  fields  opens  a  small  bandgap. 

The  tunneling  current  of  MoTe2  and  a  2D  semiconductor 
crystal  with  Eg=  1.0  eV  and  m*  =  0.1rao  (m0  is  free  electron 
mass)  as  a  function  of  the  voltage  at  different  temperatures  is 
shown  in  Figs.  2(b)  and  2(c).  Note  the  JjD  ~  F3/2  dependence 
at  low  voltages,  the  approximately  linear  relation  JjD  ~  F  at 
high  voltages,  and  the  rather  weak  temperature  dependence. 
Note  that  the  current  for  the  small  effective  mass  2D  crystal  is 
orders  of  magnitude  higher  than  MoTe2  in  Figs.  2(b)  and  2(c), 
even  though  their  bandgaps  are  similar.  Is  this  always  tme? 
Comparing  the  material  parameters  of  Figures  2(b)  and  2(c)), 
a  natural  question  is  the  relative  importance  of  effective 
masses  and  bandgaps.  For  III-V  3D  semiconductors,  the  effec¬ 
tive  masses  are  proportional  to  the  bandgaps,  as  would  be 
expected  from  interband  repulsive  interaction  from  basic  per¬ 
turbation  theory.22  The  equivalent  picture  is  not  clear  yet  for 
2D  semiconducting  crystals.  Therefore,  we  discuss  all  possi¬ 
bilities  by  treating  the  effective  mass  and  bandgap  as  inde¬ 
pendent  material  parameters.  Fig.  3  shows  the  interband 
tunneling  currents  in  2D  crystals  at  a  high  junction  field  for 
various  bandgaps,  plotted  for  a  range  of  effective  masses. 

As  is  evident  from  Fig.  3(a),  there  is  a  tradeoff  in  the 
choice  of  effective  mass  and  bandgap  for  maximizing  the 
tunneling  current.  Fig.  3(b)  zooms  in  to  highlight  this  cross¬ 
over.  For  high-performance  TFETs  for  digital  switching 
applications,  currents  exceeding  100  fiA/fim  are  highly  desir¬ 
able.  For  2D  crystals  semiconductors  with  bandgaps  smaller 
than  ^0.3— 0.4  eV,  a  choice  of  a  higher  effective  mass  will 
maximize  the  interband  tunneling  current,  far  exceeding  typ¬ 
ical  transistor  on-currents  for  high-performance  switching. 
But  for  larger  bandgaps,  a  lower  effective  mass  is  more  de¬ 
sirable.  It  is  essential  that  TFET  devices  switch  off \  which 
may  be  problematic  for  2D  crystals  with  very  small 
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bandgaps.  For  high-performance  TFETs,  2D  crystals  with 
bandgaps  in  the  ~0. 6-0.7  eV  range  and  effective  masses  of 
0.1 -0.5 m0  can  thus  be  potentially  very  attractive. 

On  the  other  hand,  there  are  various  low-power  applica¬ 
tions  of  wider  bandgap  2D  crystals  that  do  not  require  high 
on-currents,  as  long  as  they  can  beat  the  Boltzmann  limit  of 
60  mV/decade.  Such  applications  may  include  low-power 
sensors.23  Other  applications  are  expected  to  emerge  as  such 
devices  become  available.  For  such  applications,  the  low 
tunneling  currents  of  TMD  2D  crystals  might  be  actually 
turned  into  an  advantage.  We  summarize  this  work  by  hop¬ 
ing  that  the  analytical  evaluation  of  tunneling  currents  in  2D 
crystals  will  be  found  useful  for  guiding  the  choice  of  the 
right  materials  for  the  specific  applications.  The  analytic 
expression  of  the  current  vs.  voltage  can  also  form  the  back¬ 
bone  for  compact  modeling  and  design  of  2D  crystal  TFETs 
when  combined  with  the  device  electrostatics. 
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devices.  Furthermore,  separation  and 
recombination  of  photoinduced  electron- 
hole  pairs  at  a  heteroj  unction  interface 
are  primary  mechanisms  governing  the 
operation  of  solar  cells,  photodetectors 
and  light- emitting  devices.  Both  Gong  and 
collaborators  and  Duan  and  colleagues 
managed  to  electrically  address  the  lateral 
heterojunctions  and  perform  transport 
measurements  to  investigate  their  electrical 
behaviour.  Even  though  the  contacted 
structures  are  not  predefined  in  shape 
and  thus  are  far  from  optimized  devices, 
rectifying  behaviour  across  the  junctions 
was  clearly  observed,  showing  a  concrete 
proof  of  principle.  Both  groups  also  showed 
clear  evidence  of  the  photoresponse  of 
the  heterojunctions  in  their  transport 
measurements.  Similarly,  charge 
recombination  at  the  lateral  heterojunctions 
has  been  characterized  in  the  three 
studies  by  means  of  photoluminescence 
imaging,  which  revealed  an  enhanced 
emission  response  at  the  interline  between 
the  materials  in  all  cases.  Mastering  the 


in-plane  growth  of  different  TMDs  down 
to  the  single -atom  level  could  enable  the 
synthesis  of  laterally  quantum -confined 
systems,  with  promise  of  exciting  physics. 
Also,  this  could  open  up  opportunities 
for  a  wide  range  of  devices,  including 
complementary  logic  circuits,  high- 
frequency  devices  and  photodetectors. 

Whereas  the  lateral  junctions  show  scope 
for  miniaturization,  the  development  of 
larger  scale  vertical  heterojunctions  will 
prove  crucial  for  commercial  applications 
in  light  emission  and  light  harvesting. 
Large-scale  photodiodes  using  vapour- 
phase-grown  MoS2  transferred  to  silicon 
wafers  have  been  realized  already12.  They 
have  shown  impressive  performance  and 
have  also  revealed  that  using  vapour-phase- 
grown  materials  allows  tuning  of  the  optical 
response  of  the  device.  However,  the  in  situ 
growth  of  a  bilayer  heterojunction  could 
largely  improve  on  this,  in  particular  by 
engineering  pristine  interfaces  that  are  not 
achievable  by  transfer  techniques.  Further 
optimization  of  the  CVD  process  will  lead 


to  better  control  over  the  grown  structures 
and  over  the  interface  quality.  Future  work 
will  extend  these  approaches  towards  other 
material  sets  with  arresting  discoveries 
yet  to  come.  □ 
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2D  CRYSTAL  SEMICONDUCTORS 

Intimate  contacts 

High  electrical  contact  resistance  had  stalled  the  promised  performance  of  two-dimensional  layered  devices. 
Low-resistance  metal-semiconductor  contacts  are  now  obtained  by  interfacing  semiconducting  MoS2  layers  with 
the  metallic  phase  of  this  material. 

Debdeep  Jena,  Kaustav  Banerjee  and  Grace  Huili  Xing 


Transistors  and  lasers  made  of 

semiconductor  materials  power  the 
information  age  by  providing  the 
building  blocks  for  electronic  switching, 
amplification  and  photonic  communication. 
The  electronic  and  photonic  properties 
of  the  semiconductor  play  a  primary  role 
in  determining  the  performance  of  such 
devices,  yet  comparable  importance  resides 
in  the  way  this  material  interfaces  with  the 
external  metallic  circuits.  A  high  electrical 
resistance  due  to  a  high  energy  barrier  — 
the  Schottky  barrier  —  encountered  by  the 
charge  carriers  moving  through  the  metal- 
semiconductor  contacts  saps  the  energy 
efficiency,  and  significantly  degrades  the 
performance  of  the  device. 

Semiconductor  researchers  are 
excited  about  the  potential  of  transition 
metal  dichalcogenide  (TMD)  layered 
semiconductors.  For  device  applications, 
the  excitement  is  fuelled  by  the  dream 
of  improving  device  performance  by 
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using  a  semiconductor  channel  that  is 
potentially  one  monolayer  thick1.  The  first 
generation  of  field- effect  transistors  made 
from  MoS2  and  related  TMDs  have  shown 
promise,  but  trail  the  silicon  and  III-V 
semiconductor  analogues  significantly  in 
performance.  A  high  contact  resistance  is 
the  root  cause.  Writing  in  Nature  Materials , 
Rajesh  Kappera  and  colleagues  now  report 
an  unconventional  method  to  address 
this  problem  head-on2,  based  on  the  local 
conversion  of  MoS2  semiconducting  layers 
into  a  metal  that  lowers  the  barrier  to  the 
charges  flowing  in  the  external  circuits. 

The  traditional  method  to  lower  the 
contact  resistance  between  a  metal  and  a 
semiconductor  relies  on  an  intermediate 
layer,  usually  a  heavily  doped  semiconductor 
that  helps  charge  carriers  from  the  metal 
to  be  injected  in  the  conduction  and 
valence  bands  of  the  semiconductor.  The 
metal  is  chosen  based  on  its  workfunction, 
so  that  it  creates  a  low  Schottky  barrier 


height  with  the  semiconductor  contact 
region,  whereas  heavy  doping  reduces 
the  thickness  of  the  barrier,  enabling  a 
significant  fraction  of  electrons  to  quantum 
mechanically  tunnel  through  it.  However, 
chemical  doping  has  proven  challenging 
with  TMD  semiconductors.  In  their  work, 
Kappera  et  al.  exploit  the  fact  that  MoS2 
exists  in  two  crystalline  phases:  the  2H 
and  IT  types,  which  differ  in  their  stacking 
geometry.  The  IT  phase  is  obtained  by 
twisting  one  set  of  the  2H  Mo-S  tetrahedron 
by  a  60°  rotation.  The  stable  2H  structural 
phase  has  a  semiconducting  behaviour  and 
is  desired  as  the  channel  of  the  transistor. 

The  IT  phase  is  structurally  metastable,  but 
can  be  stabilized  if  electrons  are  pumped 
into  it  —  it  then  becomes  metallic  because 
of  a  half-filled  d-orbital  band.  This  electronic 
stabilization  of  the  metallic  IT  phase  can 
be  achieved  by  converting  the  2H  phase  by 
various  routes,  for  example  by  irradiation 
with  an  electron  beam3.  Kappera  and 
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colleagues  use  a  chemical  technique,  treating 
the  MoS2  with  an  organometallic  solution 
containing  n -butyl  lithium.  Lithium  donates 
electrons  to  the  2H  MoS2  converting  it  into 
the  IT  metallic  phase.  The  chemical  process 
achieves  the  phase  transition  selectively: 
covered  areas  are  left  semiconducting, 
and  exposed  areas  become  metallic,  as 
confirmed  by  extensive  chemical,  structural 
and  optical  analytical  characterization. 

When  the  solution  is  washed  away,  the 
IT  region  is  likely  to  attract  immobile 
positive  charges  and  remain  stable.  Such  a 
semiconductor-to-metal  phase-transition 
process  can  be  amenable  in  a  device 
fabrication  environment. 

The  team  then  used  this  process  to  make 
batches  of  transistors  —  the  test  structures 
have  IT  metallic  MoS2  source-drain 
contact  areas  interfacing  gold  pads,  and 
the  control  structures  have  direct  gold-2H 
semiconductor  MoS2  contacts.  In  the  test 
devices,  the  contact  resistance  was  found  to 
drop  from  -1-10  kO  pm  to  -0.2-0. 3  kD  pm. 
As  a  result,  the  performance  of  the  IT 
contact  transistors  was  significantly  superior 
across  the  board  of  metrics:  the  drive 
current,  the  sharpness  of  switching  and  the 
gain  all  improved.  The  researchers  also  find 
that  the  IT  contact  transistors  have  a  much 
higher  reproducibility  and  yield  compared 
with  their  2H  contact  counterparts.  The 
transistor  characteristics  are  also  less 
sensitive  to  the  workfunction  of  the  metals, 
suggesting  that  the  contacts’  performance 
now  mainly  depends  on  the  1T/2H  interface. 
Other  TMD  semiconductors  have  similar 
metallic  counterparts,  implying  the  same 
principle  potentially  applies  to  them. 

The  phase-engineering  approach  to 
making  low- resistance  ohmic  contacts 
to  TMD  semiconductor  materials  is  thus 
an  exciting  advance  that  addresses  a 
critical  problem  holding  back  potential 
applications.  As  with  any  new  study,  a  list 
of  unknowns  remains  to  be  worked  out. 

By  itself,  the  IT  metallic  phase  of  MoS2  is 
negatively  charged  —  meaning  the  identity 
of  the  neutralizing  positive  charges  that 
are  presumably  immobile  remains  to  be 
determined  and  controlled.  The  nature  of 
the  1T-2H  metal-semiconductor  junction, 
the  band  alignments,  and  a  potential 
way  to  contact  the  valence  band  for  hole 
conduction  need  to  be  developed.  The 
barrier  between  the  metal  and  IT  phase 
due  to  weak  van  der  Waals  bonding  needs 
to  be  investigated.  The  lateral  diffusion  of 
the  organometallic  lithium  solution  under 
the  covered  areas  can  convert  part  of  the 
desired  2H  MoS2  channel  into  the  IT  phase, 
whereas  the  chemical  conversion  of  the 
exposed  regions  is  not  perfect;  alternative 
ways  to  seal  the  channel  during  the  chemical 


n2D  (1013  cm'2) 


Figure  1 1  Contact  resistances  for  various  semiconductor  materials  against  the  quantum  limits  for 
crystalline  materials.  Using  a  IT  metallic  phase  to  interface  MoS2  with  metals  shifts  the  performance  of 
TMD-based  transistors  closer  to  the  quantum  limit  predicted  by  Landauer  and  Sharvin.  The  inset  shows  a 
typical  transistor  configuration. 


treatment  and  achieve  complete  phase 
transformation  in  the  contact  areas  need 
to  be  investigated.  And  as  the  researchers 
state,  the  stability  of  the  IT  contacts  under 
high-performance  operation  —  for  example 
when  large  currents  are  driven  through  the 
transistor  —  remains  to  be  elucidated. 

But  how  low  a  contact  resistance  can 
one  obtain?  This  problem  in  various  forms 
has  been  studied  for  more  than  a  century. 
James  Clerk  Maxwell  calculated4  the 
classical  contact’  resistance  between  two 
regions  of  conductivity  a  separated  by  an 
insulator  and  connected  by  a  conducting 
circular  constriction  of  diameter  D  to 
be  Rc  -  1  /(aD).  In  the  1950s  and  1960s, 
Landauer5  and  Sharvin6  gave  the  problem 
a  quantum  facelift.  Their  work  predicted 
a  minimum  contact  resistance  value  for 
Rc  -  h/(2e2M),  where  h  is  Planck’s  constant, 
e  is  the  electron  charge  and  M  is  the  number 
of  electron  modes  whose  wavelength  fit  the 
narrow  conductor.  In  other  words,  even 
for  a  perfect  conducting  channel  with  no 
scattering,  only  those  electron  modes  that 
fit  are  allowed  access  into  the  channel,  the 
rest  are  reflected.  As  the  sheet  density  of 
electrons  n2D  =  kf/2n  in  a  two-dimensional 
(2D)  channel  increases,  the  wavelength 
A  =  2n/kP  of  the  energetic  electrons  riding 
the  Fermi  surface  shortens,  and  more  modes 
M  ~/cF  W  fit,  where  W  is  the  width  of  the 
channel  and  kv  is  the  Fermi  wavevector. 

The  minimum  contact  resistance  is  then 
RCW  ~h/2e2k¥  -0.02 6/fnfD  kO  pm,  which 


depends  strongly  on  the  electron  sheet 
density  (in  units  of  10 13  cm'2)  in  the 
semiconductor  channel,  and  weakly  on 
some  aspects  of  the  bandstructure7.  This  is 
the  quantum  limit  of  the  contact  resistance 
for  crystalline  semiconductors,  shown  by  the 
dashed  line  in  Fig.  1. 

This  limit  has  been  experimentally 
verified  in  atomic  break-junctions  and  in 
split-gate  quantum  point  contacts,  in  which 
a  quantized  conductance  was  observed8. 

The  highest-performance  semiconductor 
transistors  are  also  grazing  this  lower  limit9, 
as  shown  in  Fig.  1.  The  latest  achievements 
for  TMDs  represent  a  major  leap,  yet  there 
is  still  room  for  improvement.  Recently,  a 
joint  academic  and  industry  research  team 
has  reported  a  method  to  chemically  dope 
MoS2  and  have  achieved  low-resistance 
contacts  and  high-performance  transistors10. 
Their  technique  is  similar  to  the  traditional 
method  used  for  semiconductors,  and  the 
contact  resistance  values  are  similar  to  the 
IT  contacts  discussed  here.  Both  these  old 
and  new  approaches  mark  important  steps 
towards  harnessing  the  innate  potential 
of  2D  crystal  semiconductors.  More 
importantly,  metal- semiconductor  junctions 
can  enable  a  host  of  unanticipated  physical 
phenomena  exploiting  the  d- orbital  pedigree 
of  conduction  electrons  in  TMDs.  □ 
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BIOMINERALIZATION 

Crystals  competing  for  space 

Analysis  of  the  growth  patterns  of  calcitic  prisms  within  the  shell  of  the  fan  mollusc  Pinnanobilis  shows  that  growth 
can  be  predicted  using  grain  theory  and  that  the  organic  casings  of  the  prisms  set  the  thermodynamic  boundaries. 

Nico  A.J.M.  Sommerdijk  and  Maggie  Cusack 


Biominerals,  such  as  bone,  teeth  and 

shells,  are  organic-inorganic  composite 
materials  with  often  amazingly 
complex  shapes  and  structures.  Biology 
uses  a  variety  of  minerals  to  form  these 
structures  —  most  commonly,  calcium 
carbonate  (CaC03)  and  calcium  phosphate, 
which  are  used  by  marine  invertebrates 
and  vertebrates,  respectively,  as  well  as  iron 
oxides  and  silica1.  Growth  of  biominerals  is 
thought  to  be  controlled  by  their  interaction 
with  a  complex  organic  macromolecular 
matrix  that  consists  of  a  specialized 
dynamic  assembly  of  (glyco)proteins  and 
carbohydrate  polymers.  As  a  consequence 
of  these  interactions,  a  high  level  of 
control  is  achieved  over  the  composition, 


structure,  size  and  morphology  of  the 
resulting  material2.  Moreover,  through  these 
interactions  the  structure  and  properties 
of  the  biological  hybrids  are  precisely 
tuned  towards  their  specific  functions:  for 
example,  navigation,  mechanical  support, 
photonics  and  protection  against  predation, 
often  producing  physical  structures  that 
surpass  those  of  synthetic  analogues. 
However,  the  extent  to  which  the  growth 
of  biominerals  is  controlled  by  their 
interaction  with  these  organic  matrices, 
and  whether  there  is  also  control  through 
classical  thermodynamic  parameters, 
remains  unclear. 

Writing  in  Nature  Materials,  Zlotnikov 
and  colleagues  now  shed  light  on  the  shape 


evolution  of  calcite  prisms  in  a  mollusc, 
and  show  that  the  prismatic  structure 
and  morphology  can  be  predicted  by 
classical  thermodynamic  theories3.  Using 
synchrotron-based  microtomography,  they 
analyse  the  mesostructure  of  the  prisms 
forming  the  calcitic  outer  layer  of  the  giant 
Mediterranean  fan  mollusc  Pinna  nobilis 
and  conclude  that  only  a  minimal  amount 
of  biological  control  is  used  to  create  the 
well-organized  prismatic  crystals.  Indeed, 
other  than  setting  the  boundary  conditions 
of  the  thermodynamics  involved  in  the 
process,  the  findings  of  Zlotnikov  and 
colleagues  suggest  that  the  biological 
organism  has  little  involvement  in  the 
formation  of  the  finer  structural  detail  of  the 
outer  layer  of  its  shell. 

Biomineralization  has  intrigued 
scientists  for  many  decades  and  serves  as 
a  constant  source  of  inspiration  for  the 
development  of  new  materials  with  highly 
controllable  and  specialized  properties4. 
Also,  because  biological  materials  are 
normally  synthesized  in  aqueous  media 
and  at  ambient  temperatures  —  conditions 
that  are  prerequisites  for  the  synthesis  of 
green  materials  —  a  richer  understanding 
of  the  biomineralization  process  may  open 
new  sustainable  pathways  to  materials 
with  advanced  functional  and  structural 
properties.  Moreover,  the  findings  reported 
by  Zlotnikov  and  colleagues  point  towards 
the  simpler  nature  of  the  factors  that 
control  this  particular  biomineralization 
process,  and  if  similar  strategies  could 
be  exploited  by  synthetic  chemists  the 
accessibility  of  green  synthesis  routes  could 
be  greatly  improved. 

As  a  result  of  their  abundance  in 
nature,  biological  materials  composed  of 
CaC03,  for  example  mollusc  shells,  are  the 


a  b 


Figure  1 1  Secondary  electron  images  of  the  calcite  prisms  and  organic  casings  of  the  shell  of  the  fan 
mollusc,  Atrina  vexillum.  a,  Fracture  section  of  the  shell  of  A.  vexillum  showing  the  oriented  calcite  crystals 
in  the  prismatic  layer  of  the  shell.  Scale  bar,  250  pm.  b,  The  organic  casings  in  which  the  prisms  grow. 
The  prisms  were  dissolved  by  incubation  in  acetic  acid  (1  M)  for  24  h.  Scale  bar,  50  pm.  Images  courtesy 
of  Peter  Chung  (a,b)  and  Maggie  Cusack  (b),  University  of  Glasgow. 
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